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FOREWORD 


This  report  covers  work  accomplished  under  ARO  Contract  DAAG29-79-C-0008 
during  the  period  February  15,  1979  to  September  30,  1982.  The  original  contract 
included  provision  for  studying  the  high  pressure,  high  temperature  CARS  spectra 
of  nitrogen  and  carbon  dioxide  for  the  purpose  of  accurate  thermometry  in  combus¬ 
tion  media.  Subsequent  events  mandated  an  extension  to  the  contract  with  two  aims 
in  mind.  First,  the  study  of  picosecond  pulse  techniques  in  CARS  spectroscopy 
was  added  for  the  purpose  of  resonant  signal  enhancement  over  the  nonresonant  back¬ 
ground  as  a  means  to  improve  detectivity.  Second,  modifications  to  the  UTRC  high 
pressure  facility  were  necessary  to  allow  the  acquisition  of  high  pressure  and 
high  temperature  CARS  spectra. 

The  personnel  involved  with  this  contract  have  had  responsibilities  for 
different  phases  of  the  work.  James  F.  Verdieck  did  the  high  pressure  spectros¬ 
copy  at  room  temperature  for  N2  and  C02  and  also  the  work  in  N2  at  elevated  tem¬ 
perature  up  to  30  atmospheres.  John  H.  Stufflebeam  continued  the  high  pressure 
CARS  work  under  the  contract  extension  and  obtained  the  data  on  N2  for  the  condi¬ 
tions  of  both  high  pressure  and  high  temperature.  Robert  J.  Hall  was  responsible 
for  the  theoretical  aspects  of  high  pressure  CARS  spectroscopy  including  the  com¬ 
puter  analyses  employed  to  model  the  experimental  data.  John  A.  Shirley  carried 
out  the  picosecond  CARS  studies.  Alan  C.  Eckbreth,  Manager,  Chemical  Physics, 
provided  the  overall  direction  of  the  experimental  and  theoretical  aspects  of  this 
work. 
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CARS  Diagnostics  of  High  Pressure  Combustion 


SUMMARY 


Under  Contract  DAAG29-79-C-0008  sponsored  by  the  Army  Research  Office,  the 
United  Technologies  Research  Center  (UTRC)  has  conducted  basic  research  investi¬ 
gations  into  coherent  anti-Stokes  Raman  spectroscopy  (CARS).  CARS  is  a  remote 
laser  diagnostic  technique  for  temperature  and  species  measurements  in  hostile 
combustion  environments.  As  such  it  possesses  considerable  relevance  to  the  Army 
in  the  general  areas  of  ballistics,  propulsion  and  internal  combustion  engines. 

This  final  report  describes  the  results  of  the  investigations  which  were  conducted 
in  two  specific  areas.  The  first  area  concerns  the  effects  of  high  pressure, 
specifically  the  phenomenon  of  collisional  narrowing,  on  CARS  spectra  from  which 
temperature  and  density  information  derive.  Experimental  studies  of  CARS  spectra 
were  conducted  in  N2  and  C02  in  a  heated,  high  pressure  cell.  The  experimental 
spectra  were  in  excellent  agreement  with  the  theoretical  model  developed  to  des¬ 
cribe  high  pressure  CARS  spectroscopy.  The  second  area  of  investigation  focused 
upon  the  use  of  picosecond  pulses  as  a  means  of  suppressing  the  sometimes  undes¬ 
irable  nonresonant  background  which  can  limit  CARS  species  detectivity.  The  main 
purpose  of  the  study  was  to  ascertain  whether  the  background  could  be  suppressed 
with  picosecond  pump-probe  techniques  with  less  resonant  signal  loss  than  accom¬ 
panies  the  more  widely  utilized  polarization  approach.  Using  broadband  CARS  genera¬ 
tion  and  detection,  as  normally  used  for  time-resolved  diagnostics,  pump-probe 
approaches  appear  to  possess  little  advantage  for  gas-phase  work  relative  to  the 
simpler  poiirization  techniques.  For  narrowband  CARS  work,  the  results  would  lead 
one  to  infer  less  signal  loss  with  pump-probe  techniques  vis-a-vis  the  polarization 
methods . 
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INTRODUCTION 


High  pressure  combustion  is  extremely  important  in  a  variety  of  practical 
applications  of  Army  relevance  such  as  propulsion,  airbreariiing  and  otherwise, 
ballistics,  and  diesel  engines.  Experimental  diagnostics  of  high  pressure  com¬ 
bustion  phenomena  are  important  to  gain  the  understanding  required  to  improve  and 
control  these  processes.  Very  little  is  known  about  high  pressure  combustion 
processes  due  to  the  lack  of  suitable  diagnostic  techniques.  A  vast  archive  of 
data  exists  from  insertion  of  thermocouples,  pressure  transducers,  and  gas  sampling 
probes  into  combustion  media.  However,  physical  probes  are  not  only  confronted 
with  survival  considerations,  they  are  of  questionable  utility  since  their  presence 
may  seriously  perturb  the  phenomena  under  study.  Optical  techniques  appear  ideally 
suited  to  diagnosing  such  phenomena,  but  many  are  inappropriate  due  to  poor  spatial 
resolution,  weak  signal  strength,  etc.  One  optical  technique,  coherent  anti-Stokes 
Raman  spectroscopy  or  CARS,  shows  considerable  promise  for  diagnosing  high  pressure 
combustion  phenomena  (Refs.  1-4).  CARS  is  „spatially-precise  (Ref.  5)  and  its 
coherent  or  beamlike  nature  is  quite  amenable  to  the  limited  optical  apertures  typical 
of  high  pressure  combustion  chambers.  Furthermore,  CARS  exhibits  a  nonlinear 
dependence  on  molecular  number  density  and  increases  rapidly  in  signal  strength  as 
the  gas  density  is  elevated.  CARS  has  already  been  shown  to  be  practically  applicable, 
having  been  demonstrated  in  gas  turbine  combustors  (Refs.  6-7),  internal  combustion 
engines  (Ref.  8),  propellant  flames  (Ref.  9),  and  shock  tubes  (Ref.  10). 

In  CARS,  incident  laser  beams  at  frequencies  and  (termed  the  pump  and 
Stokes,  respectively),  interact  through  the  third-order,  nonlinear  electric 
susceptibility,  x  *  to  produce  the  coherent  CARS  radiation  at  frequency  w-j  =  2w^  - 
ojo-  If  the  frequency  difference  <juj_  -  w2  coincides  with  a  Raman-active  vibrational 
or  rotational  mode  of  a  certain  species,  then  the  CARS  radiation  is  resonantly  en¬ 
hanced  and  uniquely  characteristic  of  that  species.  In  CARS,  measurements  of  tem¬ 
perature  are  performed  from  the  spectral  shape,  i.e.,  the  CARS  intensity  distribu¬ 
tion  with  frequency  (Ref.  11).  Species  density  can  be  determined  from  the  spectral 
shape  as  well  in  certain  concentration  ranges,  but,  in  general,  density  measurements 
are  based  upon  the  strength  of  the  CARS  signal  (Ref.  12).  The  intensity  of  CARS, 
however,  is  quite  linewidth  dependent.  As  the  gas  density  increases,  a  variety  of 
linewidth  phenomena  occur  such  as  pressure  broadening,  shifting,  motional  and  Dicke 
narrowing  which  complicate  the  calculation  of  the  CARS  spectrum.  Accurate  modeling 
of  high-pressure,  high  temperature  CARS  spectra  is  necessary  to  extract  temperature 
and  density  information  and,  thus,  is  fundamental  to  the  deployment  of  the  technique 
as  a  diagnostic. 

For  most  molecules  of  combustion  interest  at  atmospheric  pressure,  the  lines 
are  primarily  pressure  broadened.  The  isolated  linewidths  can  be  expected  to 
broaden  linearly  with  pressure.  At  pressures  where  significant  line  overlap  occurs, 
however,  the  character  of  the  broadening  changes,  and  the  phenomenon  of  collisional 
narrowing  sets  in  and  thereafter  governs  the  CARS  intensity  distribution.  Under 
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Army  Research  Office  sponsorship,  Contract  DAAG29-79-C-0008,  UTRC  has  been 
investigating  the  phenomenon  of  pressure-induced,  i.e.,  collisional,  narrowing  in 
N2  and  CO2.  In  N2,  where  most  of  the  effort  has  been  concentrated,  excellent 
progress  has  been  made  in  modeling  collisional  narrowing  in  CARS.  At  room  tem¬ 
perature,  CARS  spectra  of  N2  have  been  recorded  over  the  pressure  range  from  one 
to  one  hundred  atmospheres  in  a  specially-designed  high  pressure  cell  (Ref.  13). 

At  elevated  temperatures,  i.e.,  <v<  1600  K,  experiments  have  also  been  successfully 
performed  to  one  hundred  atmospheres.  Excellent  agreement  in  ^  has  been  obtained 
between  the  experimental  spectra  and  those  predicted  by  the  theory  developed  under 
the  contract.  The  CARS  spectrum  of  CO2  has  also  been  investigated  from  one  to 
sixty  atmospheres  at  room  temperature. 

It  is  also  of  interest  to  extend  CARS  studies  to  the  picosecond  time  domain 
and  exploit  potential  increased  signal  conversion  efficiency  and  possible  advan¬ 
tages  that  time  delay  techniques  afford  in  suppressing  nonresonant  background  con¬ 
tributions.  The  reason  this  approach  is  attractive  is  that  there  may  be  a  rela¬ 
tively  small  penalty  in  resonant  CARS  signal  reduction.  This  contrasts  with  sup¬ 
pression  by  polarization  orientation  in  which  there  is  a  resonant  signal  penalty 
of  a  factor  of  16  (Ref.  12). 

In  the  next  section  of  the  report,  the  results  of  high  pressure  CARS  studies 
are  presented  with  the  reader  referred  to  the  attached  appendices  for  details  and 
experimental  spectra.  A  discussion  is  also  presented  that  describes  the  colli¬ 
sional  narrowing  theory  of  CARS  that  was  developed  under  this  contract.  The  sub¬ 
sequent  section  reviews  the  objectives  and  results  of  the  picosecond  pulse  CARS 
studies.  The  work  at  UTRC,  sponsored  by  ARO,  has  been  received  with  considerable 
interest  by  the  technical  community  and  has  been  publicly  presented  and/or  published 
on  numerous  occasions  as  summarized  in  Appendix  B.  A  compilation  of  major  papers 
published  by  this  group  as  a  result  of  support  by  the  ARO  through  this  contract  are 
included  as  Appendices  C-I  at  the  end  of  the  report. 
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HIGH  PRESSURE  CARS  SPECTROSCOPY 


The  objective  of  work  performed  under  this  contract  was  to  extend  the 
applicability  of  CARS  as  a  diagnostic  tool  to  the  regime  of  high  pressure,  high 
temperature  combustion  media.  This  was  accomplished  by  the  acquisition  of  high 
quality,  moderate  resolution  CARS  spectra  which  were  used  to  modify  existing 
(isolated-line)  computer  models  for  the  prediction  of  temperature  from  the  CARS 
spectra  of  high  pressure,  high  temperature  gases.  Measurements  on  nitrogen  from 
one  to  100  atmospheres  from  300K  to  over  1600K  have  been  performed  as  well  as 
measurements  on  CO2  from  one  to  65  atmospheres. 

Experimental  Apparatus 

To  obtain  fundamental  CARS  spectroscopic  data  on  the  two  molecules,  a  static, 
internally-heated,  high  pressure  vessel  schematically  dep'cted  in  Fig.  1  was  utilized. 
It  is  rated  for  temperatures  to  1750K  and  pressures  to  5000  psig. 

Figure  2  is  a  layout  of  the  high  pressure  CARS  apparatus .  The  output  of  a  Quanta 
Ray  Nd:YAG  laser  was  frequency  doubled  to  generate  a  horizontally  polarized,  "pri¬ 
mary"  pump  beam  at  5320  A  (w^).  Residual  1.06  p  was  separated  from  the  ujj  beam  in 
a  prism  harmonic  separator  and  doubled  to  generate  a  secondary  beam  to  pump  the 
broadband  Stokes  dye  laser  oscillator.  Part  of  the  primary  beam  (30%)  was  split  off 
to  pump  the  amplifier  of  the  dye  laser,  Galilean  telescopes  were  provided  to 

control  beam  waists  and  the  focal  zone  locations  of  and  which  were  combined 
collinearly  on  the  dichroic  mirror  and  focused  inside  the  high  pressure  test  cell. 
CARS,  at  u)-}  was  generated  in  the  focal  volume,  and  all  three  frequencies  (u^,^,^) 
were  recollimated  after  exiting  the  cell.  A  second  dichroic  mirror  separated  W3 
from  003,  102  before  the  signal  was  incident  on  the  slit  of  a  double,  one-meter  mono¬ 
chrometer  with  resolution  of  ^  0.4  cm-^.  The  and  beams  were  trapped  after 
reflection  from  the  second  dichroic.  A  He-Ne  laser  is  shown  in  Fig.  2  whose  output 
was  coincident  with  the  003,  w 2  path  and  used  for  optical  alignment.  A  beamsplitter 
in  the  leg  provided  a  reference  signal  to  a  photomultiplier  which  was  used  to 
normalize  the  spectrally-dispersed  CARS  signal  from  the  monochrometer.  The  two  sig¬ 
nals  were  ratioed  in  a  boxcar  averager  to  account  for  variations  in  laser  intensity 
during  the  scan  of  the  CARS  spectrum.  A  photograph  of  the  L'TRC  High  Pressure 
Facility  is  shown  in  Fig.  3  with  many  of  the  components  identified  that  were  des¬ 
cribed  in  the  previous  discussion. 

Experimental  Results 

Using  the  high  pressure  facility,  studies  were  performed  from  1  to  100 
atmospheres  and  Appendix  D  contains  a  reprint  which  summarizes  the  300  K,  N2  CARS 
spectra  obtained  from  this  pressure  range.  ”he  important  result  of  this  work  was 
that  the  total  Q-branch  bandwidth  has  little  variation  with  pressure,  an  effect 
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FIG.  3 
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of  collisional  narrowing.  In  addition  to  the  discussion  of  pressure  induced 
narrowing  in  Appendix  D,  further  developments  of  this  theory  are  presented  in 
Appendices  C  and  C,  as  well  as  in  the  theoretical  section  of  this  report. 

Carbon  dioxide  was  also  studied  with  the  same  high  pressure  apparatus  over 
the  range  of  1  to  65  atmospheres.  Shown  in  Fig.  4  are  300  K  CARS  spectra  of  the 
v-p  symmetric  stretch  mode  of  CO2  at  one  and  sixty-five  atmospheres;  again,  little 
broadening  of  the  spectra  was  evident.  Due  to  the  very  small  rotational  spacing 
of  lines  in  a  vibrational  band  of  CO2,  collisional  narrowing  is  likely  to  be  quitp 
important  even  at  atmospheric  pressure.  This  was  supported  by  the  data  of  Fig.  5 
where  the  pressure  dependence  of  the  CO2  CARS  bandwidth  is  shown  to  be  very  weak. 

Also  shown  in  Fig.  5  is  the  isolated  line  prediction  of  the  bandwidth  dependence. 

This  model  assumes  isolated  rotational  lines  broaden  linearly  with  pressure  and  no 
interference  effects  occur  as  lines  overlap.  Since  bandwidth  is  used  as  a  parameter 
to  predict  temperature,  this  pressure  effect  clearly  has  a  large  impact  on  thermometry 
The  slight  discrepancy  between  the  rotational  diffusion  model  predictions  and  the 
experimental  bandwidth  is  commented  upon  in  the  section  on  Theoretical  Considerations, 
later  in  this  report. 

Studies  were  also  performed  at  high  temperature  over  the  pressure  range  pre¬ 
viously  stated.  Details  of  the  experiments  and  spectra  in  N2  at  8  and  30  atmos¬ 
pheres,  1630  K  are  described  in  Appendices  E,  F,  H.  The  theoretical  comparisons 
based  on  the  theory  of  collisional  narrowing  agree  very  closely  with  the  experimental 
data.  Also  shown  and  discussed  in  these  same  appendices  are  high  temperature  CO2 
spectra  taken  from  a  CO/air  flame  at  1620  K,  one  atmosphere.  The  high  temperature 
CO2  spectra  exhibited  several  "hot  bands",  strongly  dependent  on  temperature  and 
very  useful  for  thermometry.  The  agreement  between  theory  and  experiment  was  fairly 
good  and  has  improved  as  more  reliable  data  on  rotational  spectral  constants  have 
become  available. 

At  elevated  temperature,  and  pressures  above  30  atmospheres,  considerable  ex¬ 
perimental  difficulties  were  encountered  in  attempts  to  acquire  high  quality  CARS 
spectra.  These  difficulties  are  briefly  summarized  below,  together  with  the  modi¬ 
fications  necessary  to  overcome  them.  The  initial  design  of  the  high  pressure  cell 
did  not  include  the  sapphire  rods  shown  in  Fig.  1.  With  the  cell  mounted  horizontally 
density  gradients  induced  by  buoyancy  led  to  refractive  steering  of  the  optical  rays 
and  the  laser  beams  did  not  propagate  through  the  cell.  Instead,  they  were  steered 
into  the  cell  side  wall.  By  orienting  the  optical  axis  vertically,  i.e.,  aligned 
with  the  gravitational  field  which  is  responsible  for  the  buoyancy  effects,  the 
laser  beams  could  be  propagated  through  the  cell,  and  indeed,  one  could  look  through 
the  cell.  However,  buoyancy  forces  again  led  to  two  different  problems.  First, 
buoyancy-driven  turbulent  convection  cells  were  established  which  seriously  per¬ 
turbed  the  optical  phase  of  the  propagating  laser  beams  and  led  to  a  time-varying 
combination  of  steering,  focusing  and  defocusing.  At  high  pressures  and  high  tem¬ 
peratures,  the  ^2  in  the  cell  literally  took  on  the  appearance  of  a  stirred  "oil  and 
vinegar"  solution.  Much  of  the  laser  beam  coherence  and,  hence,  its  focusabilitv 
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was  destroyed.  Second,  a  strong  and  highly-peaked  temperature  distribution  was 
produced  whose  spatial  profile  was  pressure  dependent  as  determined  by  thermocouple 
measurements.  At  the  constant  pressure  prevailing  in  the  cell,  the  temperature  dis¬ 
tribution,  which  peaked  slightly  above  the  center  of  the  cell,  led  to  a  concomitant 
density  distribution.  The  density  was  highest  to  the  outer  edges  of  the  cell  near 
the  windows  and  lowest  near  the  middle  of  the  cell.  Ascertaining  the  temperature 
at  any  specific  location  in  the  cell  was  easily  solved  by  placing  a  stepped  ceramic 
liner  into  the  cell  (see  Fig.  1)  allowing  the  optical  pyrometer  to  focus  on  any 
particular  step  and,  thus,  a  particular  location.  However,  to  generate  CARS  at  a 
specific  temperature  location  required  the  CARS  measurement  to  be  spatially-precise. 
The  normally-employed  collinear  phase  matching  was  confronted  with  its  classic 
limitation  (Ref.  5).  Namely,  most  of  the  CARS  signal  was  generated  from  the  cooler, 
high  density  outer  regions  and  not  from  the  hot,  low  density  inner  zone  where 
measurements  at  high  temperature  were  desired.  Thus,  in  examining  high  temperature 
CARS  spectra  in  N2  with  increasing  pressure,  the  "hot  band"  in  the  N2  CARS  spectrum 
was  found  to  "disappear".  Although  this  could  by  physically  real,  it  was  suspected 
that  the  spatial  resolution  of  the  collinearly  phase-matched  CARS  was  degrading  with 
increasing  pressure.  This  smearing  out  of  the  resolution  resulted  in  (actually  from) 
greater  CARS  generation  from  the  cooler,  high  density  regions  and  a  spectrum 
not  representative  of  the  small,  high  temperature  region. 

Due  to  the  limited  optical  aperture  of  the  high  pressure  cell,  only  very  closely 
packed,  from  an  angular  sense,  crossed-beam  phase-matching  schemes  were  possible. 

A  two-beam,  folded  BOXCARS  approach  was  attempted  in  which,  in  the  near  field,  the 
Stokes  laser  beam  is  inserted  inside  the  annular  beam  (Ref.  14).  In  principle, 
this  would  avoid  w-^,  w 2  beam  overlap  except  in  the  focal  region,  obviating  the 
difficulties  of  collinear  CARS  and  leading  to  good  spatial  resolution.  Unfortunately, 
due  to  the  severe  refractive  index  fluctuations,  the  Stokes  laser  could  not  be  con¬ 
fined  within  the  inner  annular  region  of  the  pump.  Rather,  beam  overlap  occurred 
and  one  obtained  a  rather  haphazard  series  of  collinearly  phase-matched  interactions 
with  poor  spatial  resolution.  Indeed,  signal  strength  observations  did  not  demon¬ 
strate  a  pressure  dependence  anywhere  near  quadratic,  indicating  that  only  a  frac¬ 
tion  of  the  respective  beams  were  wave  mixing  from  a  geometric  sense.  Visual  ob¬ 
servations  of  the  laser  beams  emerging  from  the  high  pressure  cell  at  high  tempera¬ 
ture  also  illustrated  that  the  laser  beam  integrity  had  been  greatly  deteriorated 
with  the  laser  light  highly  dispersed  and  fluctuating  quite  severely  spatially  (of 
the  order  of  a  beam  diameter) . 

In  an  effort  to  retain  use  of  the  high  pressure  cell,  sapphire  rods  were  ob¬ 
tained  and  inserted  into  the  cell  at  locations  shown  in  Fig.  1  to  restrict  the 
gaseous  sample  to  approximately  7  cm.  By  restricting  the  optical  path  of  the  lasers 
through  the  gaseous  sample,  the  deterioration  of  the  laser  beams  was  greatly  reduced, 
and  hence,  the  spatial  precision  of  the  phase  matching  was  improved.  Additionally, 
the  gaseous  sample  should  have  been  considerably  more  uniform.  However,  experiments 
indicated  that  turbulence  still  inhibited  reliable  CARS  generation  with  the  cell  in 
the  vertical  position.  It  was  subsequently  found,  though,  that  with  the  cell  mounted 
horizontally  and  the  sapphire  rods  in  place,  the  "gas  lens"  effect  was  reduced 


12 


R82-954566-F 


significantly.  The  laser  beams  would  propagate  completely  through  the  cell  even 
for  the  most  severe  case  of  high  temperature  and  high  pressure,  and  CARS  generation 
was  observed.  However,  optical  damage  to  the  sapphire  rods  was  often  problematical. 

It  was  determined  that  the  optical  damage  was  the  result  of  time-dependent, 
thermal  lensing  in  the  absorptive  neutral  density  filters  used  to  attenuate  the  pump 
and  Stokes  lasers.  At  high  gas  pressures  it  is  necessary  to  attenuate  the  lasers 
to  avoid  gas  breakdown  and  stimulated  Raman  gain  which  would  interfere  with  and 
perturb  the  CARS  signal.  As  the  filters  absorb  energy,  they  heat  up  nonuni formly 
and  "lens"  with  the  result  that  the  focal  point  of  the  laser  beam  shifts.  Another 
approach  would  be  to  lower  the  flashlamp  energy  input  to  the  Nd:YAG  laser,  thus 
lowering  the  output  power;  however,  the  Nd:YAG  rods  also  exhibit  thermal  lensing 
and  the  laser  must  be  operated  at  a  constant  setting  to  maintain  unvarying  beam 
divergence  and  constancy  of  the  focal  location.  The  simple  solution  was  to  achieve 
the  desired  attenuation  through  a  Fresnel  reflection  (%  1.5%)  from  a  flat  surface. 

A  prism  was  used  for  this  purpose,  as  shown  in  Fig.  1,  to  avoid  the  back  surface  re¬ 
flection  which  would  result  from  a  plane  parallel  element.  Incorporating  this  modi¬ 
fication  and  the  cell  mounted  horizontally  with  the  sapphire  rods,  high  quality  CARS 
spectra  were  generated  at  high  pressure  and  high  temperature  as  shown  in  Fig.  6. 

This  was  the  first  time  spectra  have  been  acquired  at  these  elevated  pressures  and 
temperatures  from  a  static  cell.  The  103  atmosphere,  1630  K  data  showed  a  prominent 
hot  band,  confirming  the  suspicions  about  the  effect  of  turbulence  and  spatial 
resolution  discussed  previously.  It  also  allowed  verification  of  the  collisional 
narrowing  theory  in  the  regime  of  simultaneous  high  temperature  and  pressure.  As 
seen  in  Fig.  6,  the  theoretical  fit  is  excellent.  A  further  discussion  of  the 
collisional  narrowing  theory,  the  historical  precedence  in  Raman  spectroscopy,  its 
development  through  several  models,  including  'inverse  power  law',  Gordon  rotational 
diffusion,  etc.,  and  its  application  to  the  spectra  reported  here  are  detailed 
in  the  following  paragraphs. 


Theoretical  Considerations 

Prior  to  the  initiation  of  this  contract,  most  diagnostic  CARS  experiments 
and  associated  modeling  were  performed  at  low  pressures,  one  atmosphere  typically, 
where  it  is  a  reasonable  approximation  to  assume  that  the  CARS  signatures  consist 
of  a  collection  of  independent,  isolated  transitions.  Because  CARS  is  a  nonlinear 
process  whose  strength  is  proportional  to  the  squared  modulus  of  the  third-order 
nonlinear  electric  susceptibility,  it  was  known  that  interference  effects  arising 
from  line  overlap  could  make  very  important  contributions  to  the  CARS  signature.  If 
the  transition  line  positions  and  Raman  linewidths  are  known  to  a  reasonable  pre¬ 
cision,  however,  these  interference  effects  can  be  accurately  accounted  for.  At 
higher  pressures,  it  was  recognized  that  the  assumption  of  independently  broadened 
transitions  would  no  longer  suffice,  and  that  a  reformulation  of  the  expression  for 
the  resonant  contribution  to  the  nonlinear  susceptibility  would  be  necessary. 
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FIG.  6 
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It  has  long  been  known  in  such  diverse  spectroscopies  as  infrared,  Raman, 
microwave,  atomic  hyperfine,  and  ultrasonics,  that  the  mathematical  treatment  of 
overlapping  spectral  lines  is  not  adequately  given  in  terms  of  isolated  lines 
undergoing  ordinary  pressure  broadening  (Ref.  15).  The  mathematical  formulation 
of  the  overall  bandshape  resulting  from  line  overlap  was  first  given  by  Baranger 
(Ref.  16)  and  Kolb  and  Griem  (Ref.  17).  The  central  result  of  these  analyses  is  the 
need  to  specify  a  so-called  relaxation  or  transition  rate  matrix  whose  diagonal 
elements  are  related  to  the  width  and  shift  of  isolated  lines,  and  whose  off- 
diagonal  elements  control  the  rate  of  intensity  transfer  between  adjacent  lines. 

When  the  pressure  has  increased  to  the  point  where  adjacent  lines  are  substantially 
overlapped,  these  off-diagonal  elements  give  rise  to  a  collapse  or  coalescence  of  the 
collection  of  lines  into  an  unresolved  bandshape  whose  width  can  actually  decrease 
with  increasing  pressure.  This  is  the  well-known  phenomenon  of  collisional,  or 
motional  narrowing. 

Under  this  contract,  the  existing  analyses  of  the  collisional  narrowing  effect 
were  extended  to  describe  the  effect  in  spectroscopies  derived  from  the  third-order 
electric  susceptibility.  A  generalized  expression  for  the  susceptibility  was  derived 
which  remains  valid  for  substantial  line  overlap.  For  the  vibration-rotation  transi¬ 
tions  of  main  interest  in  combustion  research,  it  is  most  likely  a  good  assumption 
to  ignore  vibrational  processes  such  as  vibrational  dephasing  and  vibrational  relaxa¬ 
tion,  and  furthermore  to  neglect  any  variation  of  the  rates  of  inelastic  rotational 
energy  transfer  with  vibrational  quantum  number.  Given  these  approximations,  the 
important  off-diagonal  linewidth  parameters  are  given  simply  by  the  first-order 
rate  coefficients  for  rotationally  inelastic  energy  transfer. 

Certain  relationships  involving  the  off-diagonal  relaxation  rate  parameters 
are  known;  for  example,  scattering  matrix  unitarity  requires  that  the  sum  of  off- 
diagonal  elements  in  a  given  column  equal  the  isolated  linewidth  in  the  column,  and 
detailed  balance  relationships  connect  the  values  for  forward  and  reverse  processes. 
However,  this  is  not  generally  enough  information  to  specify  the  entire  matrix. 
Unfortunately,  it  is  not  presently  possible  to  calculate  these  off-diagonal  elements 
on  an  ab  initio  basis  with  the  same  semi-classical  theories  that  have  been  very 
successful  in  predicting  the  widths  of  isolated  lines  (Ref.  18).  While  this  for¬ 
midable  problem  is  being  resolved,  it  lias  been  necessary  to  resort  to  parameteriza¬ 
tion  of  the  off-diagonal  elements.  In  Ref.  13,  it  was  shown  that  good  results  could 
be  obtained  by  assuming  an  "exponential  gap"  relations! ip  between  relaxation  matrix 
element  and  the  rotational  energy  defect  for  the  associated  inelastic  energy  trans¬ 
fer  process. 

Use  of  the  generalized  expression  for  the  third-order  electric  susceptibility 
does  involve  some  calculational  difficulties,  because  at  each  frequency  mesh  point 
of  interest  a  so-called  "G-matrix"  related  to  the  aforementioned  relaxation  matrix 
must  be  inverted.  This  process  can  be  quite  cumbersome  and  costly  in  computer  run 
time  for  those  molecules  whose  CARS  spectra  consist  of  very  large  numbers  of  Q-branch 
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transitions.  Example  of  the  latter  are  CC^,  H20,  and  even  N2  when  the  temperature 
is  quite  high.  A  search  for  a  somewhat  simpler  approach  to  be  applied  to  high  pres¬ 
sure  H2O  CARS  spectra  led  to  investigation  of  the  Gordon  "Rotational  Diffusion" 
model  (Ref.  19).  This  research  was  not  directly  supported  under  this  contract, 
but  since  it  represents  just  a  special  case  of  the  more  generalized  susceptibility 
expression,  and  has  been  used  in  the  analysis  of  spectra  obtained  under  this  con¬ 
tract,  it  forms  an  important  part  of  any  discussion  of  results.  This  model  is 
based  on  a  picture  of  molecular  collisional  processes  in  which  strong  collisions 
simply  cause  a  thermal  randomization  of  the  rotational  angular  momentum.  Based  on 
this  assumption,  the  relaxation  rate  off-diagonal  parameters  are  specified,  with 
no  adjustable  parameters  required.  Further,  it  is  not  necessary  to  invert  the  G- 
ir.ntrix  for  this  particular  case.  The  effect  of  narrowing  is  accounted  for  in  the 
susceptibility  algorithm  by  a  simple  summation  and  subsequent  division.  This  makes 
possible  a  drastic  simplification  of  the  numerical  task  involved  in  narrowing  cal¬ 
culations,  and  makes  possible  the  modeling  of  C02  and  H20  CARS  spectra  that  would 
otherwise  be  too  expensive. 

The  rotational  diffusion  model  gives  very  good  agreement  with  N2  CARS  spectra 
at  300  K  and  1630  K.  It  also  agrees  quite  satisfactorily  with  the  measured  pressure 
dependence  of  the  room  temperature  CARS  bandwidth  of  C02 ,  as  can  be  seen  in  Fig.  5. 

The  slight  disagreement  at  the  higher  pressures  is  most  likely  due  to  uncertainty 
in  the  Q-branch  linewidths  for  CO2;  no  data  on  these  widths  are  presently  available, 
and  the  assumption  has  been  made  that  the  Q-branch  linewidths  are  equal  to  those  of 
the  pure  rotational  S-branch.  Use  of  slightly  smaller  values  for  the  Q-branch 
linewidths,  which  would  be  physically  realistic,  would  resolve  the  disagreement  at 
high  pressure. 

The  one  criticism  that  can  be  made  about  use  of  the  rotational  diffusion  theory 
is  that  for  N2  it  does  not  give  quite  enough  narrowing  at  low  pressure  where  the 
collisional  narrowing  effect  is  just  beginning  to  become  important.  This  is  a 
small  defect,  but  it  can  be  remedied  by  parameterization  of  the  relaxation  rate 
parameters,  and  resort  to  the  full  G-matrix  inversion.  In  the  latter  stages  of 
this  contract,  further  research  was  performed  on  the  subject  of  rotational  relaxa¬ 
tion  and  particularly  the  semi-empirical  models  for  rotational  transfer  rates.  One 
model  for  rotationally  inelastic  rate  coefficients  that  has  had  a  great  deal  of 
success  in  correlating  the  vast  amount  of  experimental  data  accumulating  on  this 
subject  is  based  on  an  inverse  power  law  relationship  between  energy  defect  and 
inelastic  rate  coefficient  (Ref.  21).  That  is,  the  inelastic  energy  transfer  rate 
is  assumed  to  be  proportional  to  the  absolute  value  of  the  inelastic  energy  defect 
raised  to  a  negative  power  which  is  an  unknown  parameter.  In  this  way,  the  theory 
can  be  fitted  to  the  experiment,  and  it  has  been  found  that  a  value  of  the  power  fit 
parameter  equal  to  -.5  results  in  very  good  agreement  with  the  300  K  data  of  Ref.  13. 
Agreement  at  the  highest  pressures  (p  -  100  atm)  is  optimized  by  assuming  a  physically 
realistic  contribution  of  vibrational  processes  to  the  line  broadening  of  4%.  The 
same  values  of  the  fit  parameters  give  good  agreement  with  the  hot,  high  pressure  N2 
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spectra,  as  shown  in  Fig.  6.  Also  shown  in  Fig.  6  is  a  calculation  based  upon 
isolated  line  theory;  it  is  plain  that  any  attempt  to  analyze  spectra  at  these 
extremes  of  temperature  and  pressure  with  isolated  line  calculations  would  lead  to 
highly  inaccurate  inferred  temperatures.  The  fact  that  the  first  vibrational  hot 
band  appears  in  the  N2  spectrum  even  at  103  atmospheres  augurs  well  for  CARS 
thermometry  of  high  pressure  gases,  and  is  a  confirmation  of  the  assumption  tlut 
vibrational  rate  processes  remain  slow  relative  to  rotational  relaxation. 

The  inverse  power  law  model  has  also  given  very  good  agreement  with  inverse 
Raman  spectra  of  N2  obtained  by  G.  Rosasco  and  W.  Lempert  of  NBS  (Ref.  22).  These 
spectra,  obtained  with  extremely  high  resolution  using  cw  laser  sources  at  pres¬ 
sures  between  1/2  and  2  atmospheres,  are  very  revealing  of  departures  from  isolated 
line  behavior  caused  by  collisional  narrowing. 

Based  on  the  investigations  of  this  contract,  it  is  fair  to  state  that  the 
modeling  of  high  pressure  CARS  spectra  has  proved  to  be  tractable;  theoretical 
spectra  are  in  good  agreement  with  experiment  for  N2  over  a  wide  range  of  tempera¬ 
ture  and  pressure.  Satisfactory  agreement  can  be  obtained  using  a  very  simple 
model  of  the  rotational  relaxation  processes,  the  Gordon  "Rotational  Diffusion"  model 
No  adjustable  parameters  are  involved,  and  the  expression  for  the  third-order  sus¬ 
ceptibility  simplifies  to  one  involving  no  more  computational  expense  than  an 
isolated  line  calculation.  For  N2,  theory-experiment  agreement  in  the  low  pressure 
regime  can  be  improved  by  resort  to  a  semi-empirical  expression  for  the  off-diagonal 
relaxation  matrix  parameters;  the  inverse  power  law  relationship  is  probably  the  best 
available  at  present,  and  gives  good  agreement  with  the  N2  results.  In  C02,  the 
very  large  number  of  Q-branch  transitions  comprising  the  CARS  signature  requires 
that  the  rotational  diffusion  theory  be  employed  for  reasons  of  efficiency;  for¬ 
tunately,  the  model  also  explains  the  CC>2  symmetric  stretch  mode  pressure  dependence 
quite  accurately. 
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PICOSECOND  CARS  INVESTIGATIONS 


A  perceived  limitation  of  CARS  has  been  that  nonlinear  mixing  produces  a 
background  signal  which  interferes  with  measurements  when  the  concentration  of 
the  target  species  is  low.  Polarization  orientation  techniques  (Ref.  12)  have 
been  employed  to  suppress  the  nonresonant  background  but  the  resonant  signal  is 
decreased  as  well,  a  factor  which  can  be  critical  in  the  detection  of  minority 
species.  Kamga  and  Sceats  (Ref.  23)  introduced  a  technique  which  exploits  the 
difference  between  the  temporal  responses  of  the  target  species  and  the  background. 
They  delayed  a  second  pump  beam  to  probe  the  coherent  excitation  produced  by  pre¬ 
vious  temporally  synchronous  pump  and  Stokes  beams.  A  reduction  in  the  background 
signal  in  a  mixture  containing  10  percent  CS2  dissolved  in  toluene  in  the  liquid 
phase  was  shown.  Dzhidzhoev,  et  al.  (Ref.  24)  demonstrated  background  suppression 
in  gaseous  NH^  using  a  mode-locked  Nd:YAG  laser  at  1.06  p  along  with  an  optical 
parametric  oscillator  to  generate  a  coherent  excitation  in  the  vibrational  mode 
(3334  cm~l) .  A  portion  of  the  Nd:YAG  output  was  frequency-doubled  and  delayed  to 
suppress  the  background. 


The  investigations  reported  here  were  undertaken  to  study  background  suppres¬ 
sion  in  gaseous  N2  using  delayed  pump  techniques.  A  broadband  Stokes  laser  was 
selected  to  generate  the  entire  CARS  spectrum  simultaneously,  the  approach  typically 
employed  for  time-resolved  diagnostics.  The  following  section  briefly  reviews  pre¬ 
vious  experiments  using  a  delayed  pump  to  suppress  the  nonresonant  background  in 
CARS  measurements.  Subsequently,  the  experimental  apparatus  and  the  measurements 
are  described.  The  theory  of  transient  coherent  Raman  spectroscopy  is  outlined 
and  a  simple  model  describing  the  measurements  is  presented  and  used  to  discuss 
the  results.  Finally,  implications  for  diagnostics  are  assessed. 

CARS  Nonresonant  Susceptibility  Effects 


The  nonlinearity  of  CARS  produces  mixing  which  results  in  constructive  and 
destructive  interferences  between  nearby  resonant  transitions  and  between  the 
resonant  part  of  the  susceptibility  and  the  background.  This  can  be  seen  by 
writing  the  resonant  contribution  as  a  sum  of  real  and  imaginary  parts  and  expanding 


lx(3)  1 2  =  |  ?  (x!  +  ix'.’)  +  xnr|2 
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The  real  part  of  the  susceptibility  is  an  anti-symmetric  function  of  the  detuning 
frequency  (Am  =  wj  -  w-,  +  uij)  where  icj  is  the  Raman  transition  frequency.  The 
i*'-. inary  part  is  symmetrical,  and  the  nonresonant  susceptibility  is  largely  real 
25).  For  a  minority  resonant  species,  xnr  ‘C  EXj  »  and  the  spectrum  consists 
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of  a  constant  background  with  a  dispersion  shaped  modulation  at  the  resonance 
frequency.  When  the  resonant  species  mole  fraction  is  sufficiently  low  the  modula¬ 
tion  can  be  'lost'  in  the  background,  thereby  setting  a  maximum  detection  limit. 

Several  techniques  have  been  devised  to  suppress  the  nonresonant  background. 

The  most  widely  employed,  polarization  orientation,  makes  use  of  the  difference 
between  the  tensor  properties  of  the  electronic  contribution  to  the  third  order 
susceptibility  and  that  due  to  the  resonant  nuclear  response.  Bunkin,  et  al . 

(Ref.  26)  showed  that  the  nonresonant  electronic  contribution  is  suppressed  for 
certain  orientations  of  the  polarization  vectors  of  the  laser  fields  and  an  analyzer 
in  the  CARS  beam.  It  has  also  been  shown  however,  that  the  suppression  of  the  non¬ 
resonant  signal  is  accompanied  by  a  factor  of  sixteen  loss  in  the  resonant  mode 
signal  for  isotropic  Raman  modes  (Ref.  25).  Eckbreth  and  Hall  (Ref.  12)  found  in 
background  suppression  experiments  with  carbon  monoxide  that,  under  some  circum¬ 
stances,  the  CARS  signal  can  be  reduced  by  over  two  orders  of  magnitude. 

There  are  other  methods  for  background  suppression  (or  resonant  contribution 
enhancement).  Lotem,  et  al.  (Ref.  27)  demonstrated  a  technique  which  enhanced 
the  resonance  spectra  of  cyclohexane  with  respect  to  a  benzene  background.  Their 
technique  uses  a  third  input  frequency  u0  tuned  so  that  the  difference  frequency 
oj0-w2  is  resonant  with  another  molecular  resonance.  This  technique  has  not  been 
demonstrated  in  the  gas  phase  and  does  not  offer  complete  background  suppression. 
Electronic  resonance  CARS  has  also  been  used  to  enhance  Raman  resonances  (Ref.  28). 
This  technique  uses  lasers  which,  in  addition  to  having  a  Raman  resonance  at 
also  are  resonant  with  an  electronic  transition  in  the  probed  molecular  species 
through  a  one  photon  absorption.  This  technique  may  prove  useful  for  a  select  few 
molecular  species,  however  most  major  molecules  have  electronic  transitions  in  the 
relatively  inaccessible  spectral  region  below  2000  A  (Ref.  22). 

Pulse-Sequenced  CARS 

Kamga  and  Sceats  (Ref.  23)  first  suggested  the  pulse-sequenced  technique  for 
background  suppression.  This  technique  exploits  the  difference  in  the  response 
times  for  electronic  and  resonant  vibrational  contributions  to  the  third  order 
susceptibility.  Three  beams  are  used.  Two  at  frequencies  and  u>2  are  applied 
simultaneously  to  the  medium.  A  transient  coherent  excitation  in  the  nuclear  and 
electron  distributions  is  excited  by  these.  Collisions  dephase  the  coherent  exci¬ 
tation,  but  the  dephasing  is  much  more  rapid  for  electronic  excitation  than  for  the 
vibrational  excitation.  The  third  beam  therefore  is  delayed  until  the  electronic 
excitation  has  dephased,  the  vibrational  excitation  remains,  so  the  background  is 
effectively  suppressed.  The  dephasing  time  for  electronic  excitations  is  'v  10-1- 
10-2  picoseconds  (Ref.  29),  whereas  the  vibrational  dephasing  time  for  gases  at 
atmospheric  pressure  is  on  the  order  of  10^  picoseconds. 
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Kamga  and  Sceats  demonstrated  background  suppression  in  the  liquid  phase  using 
the  656  cm--*-  vibrational  Raman  mode  of  CS2.  Two  dye  lasers  were  synchronously- 
pumped  by  an  actively  mode-locked  argon  ion  laser.  The  Stokes  frequency  was  pro¬ 
vided  by  a  rhodamine  6G  dye  laser  tuned  to  17452  cm-^.  The  pulse  duration  was  set 
by  etalons  at  about  30  psec .  The  second  dye  laser  using  sodium  fluorescein  provided 
a  3-5  psec  pulse  tuned  to  18109  cm"^  that  was  used  to  provide  both  pump  and  probe 
pulses  using  a  beamsplitter  and  optical  delay.  The  three  beams  were  crossed  to 
satisfy  phase  matching  (Ref.  5),  and  spatial  filtering  was  used  to  select  the  CARS 
signal  due  to  the  delayed  probe  alone.  Kamga  and  Sceats  measured  the  CARS  signal 
as  a  function  of  time  delay  for  several  detuning  frequencies.  These  decay  curves 
were  used  to  synthesize  CARS  spectra  at  a  fixed  time  delay.  A  substantial  reduc¬ 
tion  of  the  background  was  observed  at  a  delay  of  20  psec  with  respect  to  the  zero 
time  delay  spectrum. 

Recently  Dzhidzhoev,  et  al.  (Ref.  24)  used  delay  techniques  to  demonstrate 
nonresonant  background  cancellation  in  gaseous  ammonia.  Coherent  excitation  in 
the  totally  symmetric  vibration  mode  at  3334  cm-1  was  generated  with  the  funda¬ 
mental  (1.06  p)  of  a  passively  mode-locked  Nd:YAG  laser  mixed  with  the  infrared 
output  of  an  optical  parametric  oscillator.  A  portion  of  the  Nd :YAG  beam  was 
doubled  and  delayed  to  probe  the  coherent  excitation. 

For  a  slightly  different  reason,  Magnitskii  and  Tunkin  recently  used  the  same 
approach  as  Dzhidzhoev,  et  al.  to  measure  transient  CARS  in  H2  from  1  to  10  atmos¬ 
pheres  (Ref.  30).  An  optical  parametric  oscillator  pumped  by  a  mode-locked  Nd:YAG 
laser  was  used  to  excite  H2  molecular  vibrations  (4155  cm-^)  by  tuning  the  0P0  to 
1.9  pm  wavelength.  The  0P0  had  a  bandwidth  of  100  cm~^.  The  Q(l)  transition  was 
isolated  with  a  monochromator  and  delay  curves  were  measured  for  various  pressures. 
Magnitskii  and  Tunkin  determined  linewidths  from  the  exponential  decay  rates,  and 
demonstrated  that  the  measured  linewidths  agreed  with  previous  measurements  showing 
Dicke  narrowing  in  H2  (Ref.  31). 

The  goal  of  the  investigation  reported  here  was  to  demonstrate  transient 
coherent  Raman  techniques  in  gaseous  N2  at  one  atmosphere  and  to  evaluate  the 
attendant  resonant  mode  signal  penalty.  These  experiments  are  reported  in  the 
following  section. 


N2  Delayed  CARS  Experiments 

The  apparatus  used  to  generate  broadband  CARS  with  a  mode-locked  train  of 
pulses  is  shown  in  Fig.  7.  The  major  components  of  this  system  are  a  mode-locked 
Nd:YAG  laser  and  a  synchronously-pumped  dye  laser.  The  Nd:YAG  output  is  frequency- 
doubled  and  used  to  provide  the  pump  frequency  and  to  synchronously-pump  a  broad¬ 
band  dye  laser,  which  provides  the  Stokes-shifted  frequency  e^. 
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Mode-Locked  Nd:YAG  Laser 

The  Nd:YAG  used  in  these  experiments  is  a  standard  Quanta-Rav  laser  (DCR1-A), 
that  was  modified  at  UTRC  for  these  experiments.  The  Q-switch  and  polarizer  were 
removed  and  a  contacted  dye  cell  was  incorporated  with  the  back  reflector.  Eastman- 
Kodak  dye  #9860  dissolved  in  1,  2  diehloroethane  is  used  in  this  cell  to  passively 
mode-lock  the  oscillator.  The  dye  solution  is  circulated  in  a  laminar  flow  through 
the  1  mm  thick  dye  cell.  A  Brewster  polarizer  is  positioned  between  the  back 
reflector  and  the  oscillator  rod. 

The  unstable  resonator  of  the  Quanta-Rav  laser  design  (Ref.  32)  was  retained. 
Other  experiments  have  reported  the  operation  of  mode-locked  unstable  resonator 
Nd:YAG  lasers  (Refs.  33-35).  The  resonator  investigated  by  Liu  and  Yen  (Ref.  34) 
was  quite  close  to  the  Quanta-Ray  design  having  comparable  values  for  the  equiva¬ 
lent  Fresnel  number,  magnification  and  cavity  length.  Reali  (Ref.  35)  noted  that 
the  dye  cell  in  contact  with  the  back  reflector  forms  a  lens  and  its  effect  must 
be  taken  into  account  for  proper  operation  of  the  laser.  Accordingly,  the  original 
4-m  radius-of-curvature  back  reflector  was  replaced  with  a  10-m  mirror. 

The  entire  mode-locked  train  is  amplified  in  the  (unmodified)  Quanta-Ray  amp¬ 
lifier  and  doubled  in  the  standard  HGR-2  doubler.  Tie  design  parameters  and  per¬ 
formance  are  summarized  in  Table  1.  Figure  8  shows  a  photograph  of  a  typical  mode- 
locked  train  of  pulses  measured  with  the  output  of  an  HP  5082-4220  PIN  Diode  (Ref.  3b  ) 
displayed  on  a  Tektronix  7104  oscilloscope.  The  pulse  duration  of  the  Nd:YAG  funda¬ 
mental  was  measured  using  two  photon  fluorescence  (Ref.  37)  in  rhodamine  6G.  A 
pulse  duration  of  35  ps  was  indicated,  similar  to  the  duration  measured  by  other 
investigators  (Ref.  33-35). 

The  performance  of  the  mode-locked  Nd:YAG  laser  can  be  compared  to  the  com¬ 
mercial  mode-locked  YAG  lasers  that  have  recently  come  onto  the  market.  Two 
suppliers  currently  offer  such  lasers,  Quantel  and  Quantronix.  The  Quantel  laser 
has  a  polarization  coupled  stable  resonator  cavity  and  features  a  design  that  per- 
ments  the  laser  to  be  readily  changed  from  Q-switched  (’v  15  ns)  operation  to 
passively  mode-locked  operation  by  inserting  a  TIR  steering  prism  into  the  cavity. 

The  Model  471  is  comparable  in  energy  to  the  Quanta-Ray  (DCR-1A) .  The  performance 
of  this  model:  pulse  duration  =  35  ps,  E(1.06  pm)  =  150  mJ,  E(5230  A)  =  %  16  mJ, 
and  7-9  pulses  per  train  at  10  Hz  repetition  rate  is  similar  to  the  performance  of 
the  Quanta-Ray  laser  mode-locked  in  these  experiments.  The  Quantronix  laser  is  not 
passively  mode- locked  but  rather  uses  an  acousto-optic  modulator  to  actively  mode- 
lock  the  laser.  The  flashlamps  of  this  laser  are  on  continuously.  The  output  con¬ 
sists  of  a  train  of  pulses  at  100  MHz  repetition  frequency  with  an  individual  pulse 
duration  of  v  100  ps.  This  type  of  operation  would  not  be  appropriate  for  these 
experiments  because  of  the  long  duration  of  the  pulses. 
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TABLE  1 

MODE-LOCKED  Nd: YAG  OSCILLATOR/AMPLIFIER 

Equivalent  Fresnel  No.,  Neq  =  1.5 

Magnification,  M  =  1.5 

Cavity  Length,  L  =  75  cm 

Pulse  duration  (1.06  pm)  %  35  ps 

Bandwidth  %  0.5  cm-^ 

tp  Av  =  18  ps-cm-^ 

Number  of  pulses  in  train  =  9-12 

Beam  divergence  %  1  milliradian  (60/70J) 

Energy  per  mode-locked  train  (mJ)  (10  Hz) 


osc/amp) 

e1.06u 

E0 . 532y 

60/0 

50 

3.5 

60/10 

60 

4 

60/40 

130 

9 

60/70 

200 

20 

R82-954566-F 


FIG  8 


Nd:  YAG  MODE-LOCKED  TRAIN 
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Synchronously-Pumped,  Broadband  Dye  Laser 

The  synchronously-pumped  dye  laser  is  shown  in  Fig.  9.  The  resonator  is 
stabilized  with  three  quartz  rods.  The  length  is  tuned  by  a  micrometer  adjustment 
at  the  back  reflector  end  of  the  cavity.  The  dye  is  longitudinally  pumped,  slightly 
off  axis  with  approximately  25  percent  of  the  Nd : YAG  output.  Rhodamine  640  (Exciton) 
dye  dissolved  in  methanol  at  a  concentration  of  7.1  x  10-^  moles/liter  is  circulated 
through  the  dye  cell.  The  output  of  the  dye  oscillator  is  amplified  in  a  longitudinally 
pumped  dye  cell  which  is  series  connected  with  the  flow  circulation  through  the 
oscillator.  Approximately  38  percent  of  the  total  output  of  the  Nd:YAG  laser  is 
used  to  pump  the  amplifier.  A  delay  line  is  used  to  synchronize  the  amplifier  pump 
pulse  train  to  the  synch-pumped  oscillator  output  pulse  train.  The  delay  is  ad¬ 
justed  to  produce  the  maximum  amplification.  Figure  10  shows  an  oscilloscope  photo¬ 
graph  of  the  2x  Nd:YAG  and  dye  laser  mode-locked  trains.  The  optical  and  cable 
delays  in  the  two  legs  were  adjusted  and  checked  so  that  Fig.  10  is  a  true  represen¬ 
tation  of  the  temporal  development  of  the  dye  pulse  train  relative  to  the  Nd:YAG 
laser  train.  The  dye  laser  output  builds  up  rapidly  and  matches  fairly  well  the 
YAG  output. 

The  synchronously-pumped  dye  laser  energy  was  about  0.15  mJ  for  the  entire  mode- 
locked  train  for  the  oscillator/amplifier  combination.  The  dye  laser  spectral  width 
was  about  130  cm-1  FWHM.  The  dye  laser  pulsewidth  was  not  measured,  but  col linear 
CARS  measurements  in  which  the  frequency  doubled  Nd:YAG  laser  output  was  delayed 
with  respect  to  the  dye  pulse  train  indicate  the  dye  pulse  duration  is  comparable  to 
the  pump  pulse. 

Flame  N2  Experiments 


Spectrally-resolved  picosecond  CARS  measurements  were  undertaken  with  a  collinear 
beam  configuration,  and  no  pump  delay.  These  tests  were  conducted  at  atmospheric 
pressure  in  room  air  and  in  a  small  laboratory  burner.  The  burner  used  in  the  tests 
was  a  7.5  cm  diameter  premixed  tube  bundle  design  burning  a  premixed  stoichiometric 
flame.  The  measured  flame  temperature  in  this  burner  is  around  1700  K.  For  these 
tests  the  collinear ly  generated  CARS  beam  was  separated  with  colored  glass,  and 
narrowband  interference  filters,  dispersed  in  a  0.5  m  monochromator  (Spex  1870)  and 
detected  with  a  photomultiplier  (RCA  31034A) .  The  spectral  resolution  was  3  cm 


There  may  be  differences  between  CARS  spectra  generated  with  conventional  Q- 
switched  multimode  and  mode-locked  lasers.  The  reason  for  this  is  the  difference 
in  the  phase  distribution  between  these  sources.  The  CARS  spectrum  results  from  a 
convolution  of  the  pump  laser  and  the  Stokes  laser.  For  a  conventional  multimode 
laser  the  modal  phases  are  randomly  distributed  and  in  computing  the  spectrum  the 
various  modal  contributions  must  be  added  up  incoherently.  That  is 
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FIG.  10 


COMPARISON  OF  2x  Nd:  YAG  AND  SYNCHRONOUSLY  PUMPED  DYE  LASER 

MODE-LOCKED  TRAINS 
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On  the  other  hand,  if  the  pump  and  Stokes  lasers  are  mode-locked,  then  the  phases 
are  determinate  and  can  be  set  equal  to  zero.  In  this  case,  a  coherent  addition 
of  the  various  polarization  amplitudes  must  be  performed,  so  that 
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These  different  convolutions  lead  to  differences  in  the  calculated  spectrum.  This 
could  not  be  confirmed  in  the  spectra  measured  here,  perhaps  because  of  the  limited 
spectral  resolution. 


Time  Delay  Measurements 

Measurements  were  made  in  which  the  difference  frequency  uh  -  oun  was  tuned  to 
the  N2  vibrational  Raman  resonance  at  2331  cm  ^ .  The  apparatus  used  is  shown  in 
Fig.  11.  The  synchronously-pumped  dye  laser/amplifier  output  is  beam  expanded  and 
sent  through  a  dichroic  mirror  where  it  is  combined  with  the  undelayed  pump  beam 

A  rotatable  flat  (RF)  permits  the  phase  matching  to  be  optimized  (Ref.  5). 

The  pump,  after  portions  of  it  have  been  split  off  to  pump  the  dye  oscillator  and 
amplifier,  is  sent  through  a  telescope  to  contract  the  beam  somewhat  and  then  is 
split  into  two  equal  beams.  One  beam  is  sent  through  an  optical  delay  which  remains 

fixed  once  the  proper  timing  with  respect  to  the  dye  laser  pulse  is  determined.  The 

other  portion  of  the  Nd:YAG  beam  (03)  is  sent  through  a  variable  delay  which  permits 

delays  up  to  330  psec.  The  delayed  pump  beam  also  reflects  from  the  dichroic  mir¬ 

ror  and  is  focused  into  the  test  region.  The  Stokes  and  pump  beams  are  blocked 
with  a  glass  filter  and  the  remaining  CARS  beams  are  sent  to  photomultiplier  tubes 
for  detection.  Interference  filters  and  cutoff  filters  in  front  of  the  photomulti¬ 
plier  tubes  provide  additional  isolation  of  the  CARS  signals.  The  PMT  signals  are 
processed  by  a  dual-channel,  digital  boxcar  averager  which  is  triggered  from  a  PIN 
diode  viewing  a  portion  of  the  dye  laser  output.  The  beams  arising  from  collinear 
phase  matching  and  BOXCARS  (Ref.  5)  phase  matching  are  detected  using  spatial 
filtering  to  separate  them.  The  collinear  CARS  is  generated  by  the  undelayed  pump 
and  the  Stokes,  while  the  BOXCARS  signal  was  generated  by  the  delayed  pump.  The 
BOXCARS  signal  vanishes  when  either  of  the  three  beams  is  blocked. 

The  measured  CARS  signals  as  a  function  of  time  delay  are  shown  in  Figs.  12  and 
13  for  pressures  of  13  and  67  atm  in  N£.  The  delayed  BOXCARS  signal  is  normalized 
to  the  undelayed  collinear  CARS  signal.  The  zero  in  the  time  delay  is  arbitrarily 
selected  corresponding  to  a  fixed  micrometer  setting,  and  has  been  maintained  con¬ 
stant  throughout  the  measurements.  For  these  gas  densities,  the  dephasing  time  is 
shorter  than  the  pulse  duration,  and  therefore  there  is  little  difference  to  be  ex¬ 
pected  between  the  two  responses.  The  two  measured  dependences  indeed  are  identical 
within  the  limits  of  scatter  in  the  data.  Except  for  the  absolute  signal  amplitude, 
a  similar  response  is  expected  for  the  nonresonant  electronic  contributions  to  the 
third  order  susceptibility.  Experiments  with  carbon  dioxide  confirmed  this. 
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CARS  delay  measurements  in  N2  at  1  atm  are  shown  in  Fig.  14.  The  data  labeled 
w/o  etalon  represent  the  Raman  signal  measured  as  a  function  of  time  delay  with  the 
synchronously-pumped  dye  laser  operated  broadband.  The  decay  rate  is  less  rapid 
than  the  higher  pressure  data  and  is  exponential  for  time  delays  between  20  and  80 
picoseconds.  This  contrasts  with  the  high  pressure  data  which  is  Gaussian  in  shape 
as  a  function  of  time  delay.  The  decay  rate  of  the  w/o  etalon  curve  in  Fig.  14 
corresponds  to  a  time  constant  of  about  30  psec,  much  less  than  the  rate  expected 
from  a  dephasing  time  of  100  psec,  i.e.,  T2/2  =  50  psec,  corresponding  to  the  Raman 
linewidth  (Ref.  38).  In  order  to  explain  this  apparent  discrepancy,  it  is  necessary 
to  examine  the  theory  of  transient  coherent  Raman  spectroscopy. 


Theory  of  Transient  Coherent  Raman  Spectroscopy 

Transient  Raman  scattering  has  been  the  subject  of  considerable  interest, 
particularly  for  the  probing  of  the  very  rapid  vibrational  interactions  occurring 
in  condensed  phases.  Coherent  pump  and  probe  techniques  for  the  transient  probing 
of  vibrational  excitations  in  liquids  were  described  in  an  excellent  review  by 
Laubereau  and  Kaiser  in  1978  (Ref.  39).  The  methods  developed  by  these  authors 
employ  an  intense  pump  beam  of  short  duration  to  excite  a  coherent  vibrational 
excitation  by  stimulated  Raman  scattering.  After  the  passage  of  the  exciting 
optical  pulses  the  remaining  coherent  vibrational  excitation  decays  by  collisions. 
The  rate  of  decay  is  measured  with  a  weak  probe  pulse  which  is  delayed  with  respect 
to  the  stronger  excitation.  The  measured  decay  rate  then  is  related  to  molecular 
processes  occurring  in  the  sample.  Laubereau  and  Kaiser  described  several  methods 
of  excitation,  with  the  most  generally  applicable  being  stimulated  Raman.  This 
technique  is  similar  to  that  described  in  this  report,  but  differs  in  the  important 
aspect  of  the  method  by  which  the  Stokes  frequency  is  generated.  In  this  report, 
the  Stokes  wave  is  generated  by  the  synchronously-pumped  dye  laser.  In  contrast, 
Laubereau  and  Kaiser  use  the  Raman-shifted  wave  generated  by  stimulated  Raman 
scattering  to  generate  the  Stokes  frequency.  The  Stokes  wave  builds  up  from  the 
spontaneous  Raman  scattering  to  relatively  high  levels  due  to  the  large  gain 
coefficient  for  intense  fields  in  concentrated  molecular  systems.  Gains  of  10^  - 
1011  are  not  unusual  for  intensities  of  10^  -  lO^  W/cm^.  For  gases,  however,  the 
gain  is  not  high  enough  for  this  approach  to  be  practical,  hence  the  synchronously- 
pumped  dye  laser. 

Laubereau  and  Kaiser  in  their  review  derive  the  coupled  differential  equations 
governing  the  evolution  of  the  Stokes  amplitude  Eg,  the  coherent  vibrational  ampli¬ 
tude  Q,  and  the  vibrational  excitation  (population)  n,  as  follows: 
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here  x'  and  t '  are  the  coordinate  and  time  in  a  frame  moving  at  velocity  v,  the 
group  velocity  of  the  Stokes  pulse.  The  coupling  coefficients  v and  s  *  arc 
given  by 


and 


E^  is  the  amplitude  of  the"  pump  laser,  and  T^  and  T2  represent  the  population  life¬ 
time  (energy  relaxation  time)  and  the  dephasing  time,  respectively.  N  is  the  total 
number  density  of  the  probed  species,  a  is  the  isotropic  part  of  the  Raman  polariz¬ 
ability,  n  is  the  fractional  population  of  the  upper  vibrational  level  at  thermal 
equilibrium  and  00  and  k  denote  respectively  the  circular  frequency  and  the  wave- 
number  of  the  waves,  denoted  by  the  subscripts  L  and  S. 

Carman,  et  al.  (Ref.  40)  gave  a  general  closed  form  solution  for  these  equa¬ 
tions.  This  solution  is  also  described  in  the  Laubereau  and  Kaiser  review  (Ref.  39) 
in  which  the  form  of  the  initial  Stokes  field  is  also  considered.  Carman,  et  al. 
and  Laubereau  and  Kaiser  found  the  solutions  for  high  gain  to  depend  only  weakly 
on  the  absolute  magnitude  of  the  initial  Stokes  fields.  In  the  case  considered  here 
the  Stokes  field  is  large,  and  will  be  taken  to  be  invariant  with  x'.  The  solution 
of  Eqs.  (4-6)  simplifies  considerably  in  this  case.  This  simplicatior.  is  justified 
since  the  gain  at  the  Stokes  frequency  is  small.  These  calculations  will  be  out¬ 
lined,  but  it  is  instructive  to  consider  first  the  physical  picture  of  delayed  CARS. 

The  model  of  the  delayed  pump  probing  of  Raman  anti-Stokes  resonances  can  be 
described  as  follows.  The  pump  and  Stokes  waves  incident  on  the  medium  in  the  focal 
volume  excite  the  resonant  molecular  vibrational  level.  Since  the  lasers  are  mode- 
locked,  they  possess  a  very  narrow  distribution  of  phases,  and  a  coherent  wave  is 
generated  in  the  material  coordinate  with  a  well-defined  phase.  The  coherence  of 
the  material  excitation  decays  by  collisions  as  determined  by  the  dephasing  time  T2. 
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The  excited  population  decay  on  the  other  hand  is  determined  by  the  vibrational 
relaxation  rate,  T^.  The  coherent  wave  then  builds  up,  from  an  initially  small 
value,  at  a  rate  determined  by  the  pulse  duration  and  intensity  of  the  pulses,  and 
proceeds  to  decay  with  a  time  constant  proportional  to  T2.  If  T2  <  Tp,  the 
coherent  excitation  decays  essentially  with  the  laser  pulses. 

The  short  time  constant  of  molecular  processes  dictates  the  use  of  mode-locked 
lasers.  The  pulse  duration  of  a  mode- locked  laser  is  determined  by  the  bandwidth 
of  the  laser  medium  and  the  number  of  cavity  modes  under  the  gain  curve.  There  is  a 
fundamental  relationship  between  the  pulse  duration  and  the  gain  bandwidth,  since 
these  variables  are  Fourier  transform  pairs.  The  product  of  the  bandwidth,  Av, 
and  the  pulse  duration,  Tp,  is  given  by  Av  Tp  -  K  (Ref.  41).  The  constant  K  depends 
on  the  shape  of  the  pulses  (Ref.  41).  For  Gaussian  pulses  K  is  0.441.  Therefore 
for  the  gain  bandwidth,  Ak,  in  wavenumbers,  Ak  Tp  -  15  ps-cm  ^ ,  for  Gaussian  pulses. 

The  finite  linewidth  laser  means  that  several  frequencies  in  the  molecule  are 
excited  simultaneously.  These  frequencies  can  interfere,  producing  beating  on  a 
time  scale  of  (v^  -  Vj)-^,  where  v^  and  Vj  are  the  frequencies  corresponding  to  the 
i-th  and  j-th  Raman  resonances.  Laubereau,  Wochner  and  Kaiser  (Ref.  42)  observed 
this  beating  in  the  coherent  excitation  and  probing  of  the  tetrahedral  Alg  vibra¬ 
tion  of  CCl^  around  459  cm-1.  Isotope  shifts  in  the  two  chlorine  isotopes  split 
the  Raman  line  into  four  main  frequency  components  separated  by  approximately 
3  cm"^.  Beating  was  observed  in  the  coherent  probe  signal  as  a  function  of  delay 
time  on  a  time  scale  of  'v  11  psec  in  good  agreement  with  the  frequency  splitting 
observed  in  the  spontaneous  Raman  spectrum. 

In  the  practice  of  CARS  spectroscopy,  there  are  two  approaches  to  generating 
the  spectrum  of  a  molecule.  One  utilizes  a  relatively  narrow  bandwidth  Stokes 
laser  that  is  scanned  in  frequency  to  generate  successively  the  molecular  spectrum. 
The  advantage  of  this  technique  is  that  very  high  spectral  resolutions  can  be  ob¬ 
tained.  The  other  approach  is  to  use  a  broadband  Stokes  laser  to  generate  the 
entire  molecular  spectrum  during  each  laser  pulse.  The  latter  approach  is  desirable 
for  practical  diagnostic  applications,  where  rapid  measurements  are  sometimes 
required.  The  broadband  approach  was  selected  here  in  order  to  compare  directly 
picosecond  CARS  diagnostics  to  conventional  CARS  measurements  with  0  (10  sec) 
pulses.  The  effect  of  beating  on  the  picosecond  CARS  spectrum  of  nitrogen  must 
also  be  considered  because  of  the  broad  Stokes  bandwidth.  The  fundamental  vibra¬ 
tional  Q-branch  band  (AJ  =  0,  v  ■  0  -*•  1)  is  split  by  centrifugal  distortion  and  the 
vibrational  frequencies  are  given  by  (Ref.  43) 

«nl(J)  *  ~  a  J(J+1>  (9) 

01  0  e 

The  effect  of  the  distribution  of  frequencies  can  be  included  by  calculating  the 
sum  of  the  various  vibrational  components.  This  sum  is  calculated  according  to 
the  following  equations. 
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where  E^  and  Eg  are  the  time  dependent  amplitudes  of  the  pump  and  the  Stokes  lasers 
respectively,  and  Qj  is  tne  vibrational  excitation  of  the  j-th  rotational  level. 

qj-  is  the  instantaneous  coordinate  of  the  j-th  wave  produced  by  F.;  and  Eg,  with 
anti-Stokes  frequency  Wj  and  phase  ij .  The  phase  constant  can  be  arbitrarily  set 
to  zero  without  loss  of  generality.  The  total  amplitude  of  the  coherent  wave  in¬ 
cluding  all  rotational  components  j  is  given  by  the  sum  in  Eq .  (12)  where  fj  is  the 
fractional  population  of  the  j-th  rotational  level. 

2 

The  total  amplitude  |Qtot|  is  not  directly  measured  in  pulse  delay  techniques 
The  ultrashort  laser  pulses  are  too  short  for  direct  observation  with  standard 
electronic  detectors.  Instead  the  time-integrated  scattering  signal  is  observed. 
This  is  given  in  Ref.  (39)  as  a  function  of  the  time  delay,  tD,  between  the  pump 
laser  pulse  and  the  delayed  probing  pump  pulse  (1^) 
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where  nAS  ls  the  index  of  refraction  at  the  anti-Stokes  frequency,  f (Ak^g)  is  a 
function  of  the  wavevector  mismatch  (Ref.  39),  and  Ej^  is  the  amplitude  of  the 
delayed  pump  signal.  The  convolution  in  this  integral  has  relatively  little  effect 
in  the  calculations  reported,  therefore,  the  coherent  excitation  description  is 
retained  in  the  results  shown  later. 


The  time  dependence  of  Qj  has  been  calculated  by  Runge-Kutta  integration  of 
Eq.  (10)  assuming  Gaussian  pulses  for  the  pump  and  Stokes  fields.  The  initial 
value  of  Qj  was  assumed  to  be  zero.  The  summation  in  Eq .  (12)  was  formed  using 
Eq.  (11)  to  calculate  the  time  dependent  waves.  The  frequencies  uij  contained  in 
Eq .  (11)  are  much  higher  than  the  time  step  At  used  in  the  calculations,  as  a  time 
step  At  *  ljjj~^  would  not  be  practical.  Note  however,  that  the  large  number  of 
frequencies  in  the  summation  smooth  the  envelope  of  the  resultant  function.  The 
result  then  is  a  function  which  has  high  frequencies  which  the  calculation  samples 
at  the  rate  of  At-*.  A  number  of  calculations  were  performed  where  the  time  step 
was  varied,  effectively  changing  the  sampling  frequency,  and  it  was  found  that  the 
substructure  varied  considerably  but  the  envelope  shape  did  not.  Accordingly,  the 
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envelope  of  the  resulting  waves  was  determined,  and  it  is  that  that  is  shown  in 
the  calculations.  Figure  15  shows  the  normalized  coherent  excitation  as  a  function 
of  T/T2>  the  time  after  the  pulse  divided  by  the  dephasing  time,  with  the  spectral 

width  of  the  dye  laser  as  a  parameter.  The  curve  for  A<g  =  0.01  cm  *  is  shown 

because  it  represents  the  temporal  behavior  of  the  individual  rotational  components. 
It  is  valid  to  give  this  curve  this  interpretation  only  becuase  the  .01  cm-^  dye 
laser  bandwidth  is  less  than  that  for  a  transform  limited  pulse  of  25  psec  as 
assumed  in  this  calculation.  The  ratio  T2^p  e9uals  4.0  indicative  of  a  dephasing 
time  of  100  psec.  This  dephasing  time  corresponds  to  the  spontaneous  Raman  line- 
width  of  0.1  cm-'*'  at  300  K  and  1  atm  (Ref.  38). 

The  dye  laser  width  is  taken  to  be  Gaussian  characterized  by  A<g  (FWHM) .  The 

effect  of  larger  dye  laser  spectral  widths  is  to  excite  more  ro-vibrational  fre¬ 
quencies  which  interfere  causing  the  total  coherent  excitation  to  decay  very 
rapidly  at  first.  At  longer  times,  since  the  individual  amplitudes  remain  rela¬ 
tively  large,  constructive  and  destructive  interferences  cause  beating.  The  beating 
period  is  irregular  because  of  the  nonlinear  dependence  of  the  frequencies  on  j . 

The  initial  decay  for  A<s  -  1  cm-'*'  very  nearly  follows  the  input  pulsewidth.  The 
initial  decay  and  interference  are  even  more  dramatic  for  the  case  where  A<c  = 

100  cm-1  as  shown  in  Fig.  16.  This  width  is  more  nearly  typical  of  the  broadband 
dye  laser  used  in  these  experiments.  The  initial  decay  is  quite  rapid  and  inter¬ 
ference  produces  irregular  beating  for  T/T2  >  0.2. 

The  temporal  dependence  of  the  coherent  excitation  shown  in  Fig.  16  is  repre¬ 
sentative  of  the  broadband  Stokes  laser  case  shown,  in  Fig.  14,  by  the  curve  labeled 
w/o  etalon.  Both  curves  show  a  rapid  decrease  in  the  signal  as  a  function  of  time 
following  the  application  of  the  excitation  pulses.  To  check  the  calculated  effect 
of  restricting  the  number  of  excited  Raman  transitions,  the  Stokes  bandwidth  was 
narrowed  by  inserting  an  etalon  into  the  synchronously  pumped  laser  cavity.  The 
etalon  had  a  free  spectral  range  of  11.4  cm-*-  so  the  dye  laser  output  spectrum  con¬ 
sisted  of  a  series  of  spikes  with  this  frequency  spacing.  An  analyzing  etalon  with 
a  FSR  of  0.67  cm-*-  indicated  that  the  frequency  spread  of  an  individual  spectral 
peak  was  less  than  this  value.  Since  the  N2  Raman  spectrum  is  less  than  10  cm  ^  at 
T  =  295  K,  only  a  narrowband  of  Raman  frequencies  is  resonant  with  the  5320  A  pump. 
The  delay  curve  for  this  case  is  shown  in  Fig.  14  as  the  data  labeled  w/etalon.  The 
decay  is  clearly  slower  in  this  case  and  is  not  exponential,  as  predicted  by  the 
calculations.  For  a  time  delay  greater  than  50  psec,  the  1  atm  signal  decays 
slowly.  This  curve  should  be  compared  with  the  high  pressure  data,  which  is  expected 
to  be  similar  to  the  decay  of  a  nonresonant  signal,  since  T2  <  Tp  in  both  cases. 
Suppose  that  N2  is  present  as  a  minority  species  in  a  mixture  such  that  at  zero  time 
delay,  the  nonresonant  signal  contribution  is  equal  to  the  resonant  contribution. 

Then  at  a  delay  of  100  psec  these  data  would  indicate  that  the  nonresonant  signal 
should  be  two  or  three  orders  of  magnitude  smaller  than  the  resonant  signal  which 
is  also  decaying  but  at  a  slower  rate.  For  this  time  delay  the  resonant  signal  would 
decrease  by  a  factor  of  about  5.  The  low  power  of  the  Stokes  laser  precludes  this 
from  being  demonstrated,  because  the  CARS  signal  is  near  the  limit  of  detection  in 
our  experiments. 
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Discussion  of  Results 


The  simple  theory  developed  for  N2  illustrates  the  importance  of  interference 
effects  in  gaseous  nitrogen  transient  coherent  Raman  spectroscopy.  This  simple 
theory  has  several  shortcomings,  however.  The  theory  does  not  consider  spatial 
variations  in  the  anti-Stokes  wave,  nor  does  it  include  what  properly  should  be  a 
convolution  between  the  Stokes  laser  and  pump  bandwidths.  Furthermore,  the  Raman 
transitions  are  treated  as  isolated,  a  fact  that  certainly  is  not  valid  for  the 
lowest  rotational  transitions  even  at  one  atmosphere.  At  higher  pressure  colli- 
sional  narrowing  effects  become  more  important  (Refs.  13  and  19)  throughout  the 
spectrum,  but  then  the  dephasing  time  will  be  less  as  well,  and  the  technique  would 
not  be  applicable. 
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The  Raman  spectrum  in  the  liquid  phase  consists  of  a  single  sharp  resonance.  Inter¬ 
ference  effects  are  not  expected  for  liquid  samples  unless  discrete  frequencies 
remain  as  in  the  case  of  the  CCl^  isotopes  as  studied  by  Laubereau,  Wochner  and 
Kaiser  (Ref.  42). 


Transient  coherent  Raman  spectroscopy  of  gaseous  ammonia  was  examined  by 
Dzhidzhoev,  et  al.  (Ref.  24).  The  rotational  structure  of  the  NH3  vibrational  Raman 
band  is  expected  to  consist  of  small  number  of  closely  spaced  transitions.  There¬ 
fore,  interference  effects  would  not  be  expected  to  be  strong  in  ammonia  even  though 
a  relatively  broadband  Stokes  source  (30  cm-''')  was  used. 

Magnitskii  and  Tunkin  (Ref.  30)  observed  delayed  probe  signals  in  H2  using 
as  a  Stokes  source,  an  optical  parametric  oscillator,  having  a  bandwidth  of  100  cm-'. 
The  Raman  Q-branch  spectrum  of  H2  consists  of  several  lines  of  alternating  intensity 
(Ref.  43)  which  would  be  excited  by  the  spectrally  broad  Stokes  beam,  so  that  beating 
is  expected  to  be  important.  This  point  though  brings  out  something  that  is  important 
to  understand.  The  coherent  excitation  of  each  individual  transition  is  the  result 
of  a  rapid  buildup  along  with  the  excitation  pulses  followed  by  an  exponential  decay 
at  a  rate  determined  by  the  dephasing  time.  Beating  arises  when  a  detector  is  used 
to  observe  simultaneously  all  the  Raman  frequencies.  The  completely  isolated  trans¬ 
ition  will  not  show  beating.  Magnitskii  and  Tunkin  observed  no  beating  because  they 
isolated  a  single  transition  (Q(l)).  This,  as  a  practical  matter,  is  more  difficult 
in  N2. 
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Conclusions 

If  a  stronger  narrowband  Stokes  laser  had  been  available,  the  measurements 
show  it  should  be  possible  to  suppress  the  nonresonant  background  by  approximately 
two  orders  of  magnitude  relative  to  the  resonant  mode  signal  near  one  atmosphere 
pressure.  At  the  same  time,  the  loss  in  the  resonant  signal  would  be  less  than  an 
order  of  magnitude.  The  short  dephasing  time  at  higher  pressures  however,  would 
limit  the  utility  of  the  transient  Raman  pump  and  probe  technique  in  high  pressure 
media  for  nitrogen-like  molecules,  due  to  the  fundamental  limitation  of  the  pulse 
duration- frequency  bandwidth  product. 

For  optimum  background  cancellation  the  pump  and  Stokes  lasers  should  both  be 
operated  with  a  narrow  bandwidth  and  a  short  pulse  duration  to  minimize  the  number 
of  Raman  frequencies  excited.  This  means  the  CARS  sources  should  be  transform- 
limited.  In  this  regard,  the  approach  using  an  optical  parametric  oscillator  may 
give  a  higher  energy  conversion  to  the  Stokes  frequency,  but  the  pulses  are  not 
generally  transform-limited,  and  therefore  may  not  be  appropriate  for  N2-like 
molecules. 

The  transient  Raman  pump  and  probe  technique  promises  to  produce  excellent 
results  for  liquid  materials  as  has  been  demonstrated  by  Laubereau  and  Kaiser  and 
their  co-workers.  The  technique  may  also  be  useful  for  selected  gases  at  atmos¬ 
pheric  pressure.  At  high  pressure  the  technique  is  not  generally  applicable. 

On  the  other  hand,  the  signal  loss  associated  with  polarization  orientation  sup¬ 
pression  would  be  more  than  compensated  by  the  signal  increase  at  elevated  pressures. 
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ABSTRACT  -  The  CARS  spectrum  of  N^  has  been  investigated  experimentally  over  the  pressure  range  of 
one  to  approximately  100  atmospheres.  The  spectra  show  evidence  of  collisional  narrowing,  and  are 
in  good  agreement  with  theory. 

__  '  ff  t  t 

RESUME  -  Le  spectre  CARS  de  N2  a  ete  recherche  par  experience  pour  des  pressions  entre  une  et  cent 
atmospheres.  Les  spectres  demont  retrecissement  a  cause  des  collisions  et  sont  bien  d'accord  avec 
la  theorie. 


By  solving  the  equation  of  change  for  the 
density  operator  with  the  interaction  energy  -  4oE2 , 
and  with  a  collision  integral  representing  line- 
width  effects*,  the  third-order,  nonlinear  electric 
susceptibility  governing  resonant  CARS  generation 
from  vibrational  Q-branch  transitions  is  shown  to 
be  given  by 


X.  -  H  X/i 
3  ft  r  r  s  s  s 


Co) 


(1) 


where  a  is  isotropic  polarizability:  E  is  electric 
field  strength;  N  is  active  molecule  number  den¬ 
sity,  and  Ao*°*  is  the  normalized  population 
difference  for  vibration-rotation  Q-branch  Q(s). 
The  elements  of  the  G-matrix  are  given  by 


i(u  — ui  _uj  )  6 
1  *  r  rs 


♦  Op  iAr)«rS+Yrsa-{rs>  <2> 


where  and  are  the  frequencies  of  the  pump 
and  Stokes  sources,  respectively:  «r,  T  ,  and  &r 
are  the  vacuum  frequency,  isolated  linewidth,  and 
shift,  respectively,  of • transition  Q(r),  and  the 
off-diagonal  linewidth  elements  y  govern  the 
rates  of  collisional  energy  transfer  between  lines. 
The  presence  of  the  yrs  causes  a  change  in  the 
broadening  of  the  CARS  signatures  at  pressures 
where  individual  N2  Q-branches  have  begun  to  over¬ 
lap.  Further  increases  in  pressure  are  expected 
to  result  in  a  continuous  narrowing  01  ~i e  overall 
CARS  signature,  a  manifestation  of  the  phenomenon 
of  collisional,  or  pressure-induced  narrowing.*’’2 
Experimentally,  a  high  pressure  vessel  with  co¬ 
axial  optical  access  has  been  employed  to  study 
the  300°K  CARS  spectrum  of  No  over  the  pressure 
range  one  to  95  atmospheres.  The  pump  source 
was  generated  by  frequency-doubling  the  output  of 
a  pulsed  Nd:YAG  laser;  the  broadband  Stokes  source 
was  provided  by  a  dye  oscillator-amplifier  driven 
by  a  portion  of  the  pump  output.  The  pump  laser 
and  dye  laser  beams  were  combined  by  means  of  a 


dichroic  mirror  and  focussed  in  the  high  pressure 
test  cell,  generating  coaxial  CARS.  The  pump 
beam  was  attenuated  with  neutral  density  filters 
in  order  to  eliminate  stimulated  Raman  scattering, 
which  can  be  a  problem  at  high  pressures  with 
unattenuated,  high  intensity  beams.  At  the 
higher  pressures,  the  CARS  signature  was  readily 
visible  as  a  bright  blue  spot.  Detection  of  the 
CARS  signal  was  carried  out  using  a  one  meter 
Jobin-Yvon  monochromator  fitted  with  an  optical 
multichannel  analyzer,  providing  an  instrumental 
resolution  of  about  1.5  cm"*. 

As  a  function  of  increasing  pressure,  the 
width  of  the  N2  CARS  signature  was  found  to 
decrease  continuously  from  one  to  about  50 
atmospheres,  with  little  variation  in  the  range 
50 -95  atmospheres.  Calculations  based  on 
Equation  (1)  were  found  to  be  in  good  agreement 
with  these  results  over  the  entire  pressure 
range.  Figure  1  displays  theory-experiment 
comparisons  at  one  and  95  atmospheres.  The 
significant  role  that  collisional  narrowing  plays 
at  95  atmospheres  can  be  appreciated  by  consider¬ 
ing  the  fact  that  the  Isolated  linewldths  for  N2 
Q-branches  are  0(10  cm-*)  at  this  pressure^,  and 
the  overall  width  of  the  CARS  spectrum  would  be 
at  '.east  10  cm~*  if  the  narrowing  mechanism  were 
not  present.  As  seen,  the  width  at  95  atmospheres 
is  about  75 Z  of  the  one  atmosphere  value.  In  the 
theoretical  calculations  of  Figure  1.  the 
Isolated  Q-branch  llnewidths  of  Ref.  3  were 
scaled  linearly  with  pressure.  The  selection  rule 
AJ  »  *  2  was  assumed  in  calculating  the  Yrs: 
together  with  detailed  balance  and  S-matrix 
unitarlty  relationships,  this  results  in  a 
determinate  set  of  yrg  containing  no  adjustable 
parameters. 
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Figure  1.  Comparison  of  experimental  and  theoret¬ 
ical  N,  CARS  signatures  at  one  and  95  atmospheres. 
The  white  dots  are  experimental  points,  and  the 
theoretical  predictions  are  given  by  the  solid 
white  lines.  The  pump  laser  bandwidth  is  .4  cm 
and  the  instrumental  resolution  is  1.5  cm 
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The  coherent  anti-Stokes  Raman  spectroscopy  (CARS)  spectrum  of  N2  has  been  recorded  experimental!*  over  the  pres¬ 
sure  range  from  one  to  100  atmospheres.  The  pressure  dependence  of  the  300  K  vibrational  Q-branch  signature  is  strongly 
affected  by  collisional  narrowing.  A  theoretical  model  of  CARS  generation  for  overlapping  spectral  lines  gives  good  agree¬ 
ment  with  the  experimental  results  over  the  entire  pressure  range. 


1 .  Introduction 

Coherent  anti-Stokes  Rantan  spectroscopy  (CARS) 
is  a  nonlinear,  three-wave,  optical  mixing  technique 
that  is  being  successfully  applied  as  a  diagnostic  tool 
in  a  number  of  different  fields  [1—4] .  Of  particular 
promise  is  the  application  of  CARS  for  combustion 
diagnostics,  where  its  coherence  and  high  conversion 
efficiency  make  it  well  suited  for  spatially  and  tempo¬ 
rally  precise  probing  of  the  high  interference  environ¬ 
ments  characteristic  of  practical  devices.  The  useful¬ 
ness  of  CARS  has  been  demonstrated  by  temperature 
measurements  in  a  highly  sooting  diffusion  flame  [5], 
an  internal  combustion  engine  [6] ,  and  in  liquid-fu¬ 
elled  combustors  [7-9],  Single  pulse  spectra  have  al¬ 
so  been  obtained  in  turbulent  combustion  media  [8,9], 
offering  the  potential  for  determination  of  the  spec¬ 
trum  of  turbulent  temperature  fluctuations.  Because 
N2  is  the  dominant  constituent  of  air-fed  combustion, 
it  has  therefore  received  the  most  attention  in  the 
foregoing  experiments. 

While  most  of  these  experiments  have  been  carried 
out  at  one  atmosphere  pressure,  there  are  strong  prac¬ 
tical  motivations  for  extending  CARS  diagnostic  capa¬ 
bilities  to  higher  pressures  because  high  pressure  com¬ 
bustion  is  very  important  in  such  contexts  as  air- 
breathing  and  rocket  propulsion,  and  ballistic  com- 

*  Research  supported  in  part  by  the  Army  Research  Office, 
Durham,  N.C.  under  Contract  #DAAG29-79-C-0008, 


pression.  As  an  example,  better  efficiency  is  obtained 
from  gas  turbines  if  they  are  operated  at  higher  pres¬ 
sure.  In  addition,  studies  of  the  pressure  dependence 
of  overlapping  spectral  lines  can  reveal  information 
about  intermolecular  collisional  energy  transfer  pro¬ 
cesses,  as  will  be  seen. 


2.  Collisional  narrowing 

In  N2  CARS  thermometry,  the  signature  is  domi¬ 
nated  by  vibrational  (7-branch  transitions.  At  low  pres¬ 
sure,  where  the  homogeneous  Raman  linewidths  in 
N2  are  small  (~0.1  cm-1  atm-1)  [10],  the  CARS  in 
tensity  distribution  will  be  determined  by  frequency 
splitting  between  lines  due  to  vibration-rotation  inter¬ 
action.  The  isolated  linewidths  can  be  expected  to 
broaden  linearly  with  pressure,  but  at  pressures  where 
significant  line  overlap  occurs  the  character  of  the 
broadening  changes,  and  the  phenomenon  of  collision¬ 
al  narrowing  [11-19]  will  thereafter  govern  the  CARS 
intensity  distribution.  Because  the  interaction  of  vi¬ 
bration  and  rotation  is  so  small  in  N2,  the  vibrational 
2-branches  are  very  closely  spaced,  and  substantial 
line  overlap  will  occur  even  at  moderate  pressures. 

Collisional  narrowing  was  first  observed  experimen¬ 
tally  by  Bloembergen  et  al.  [1 1  ]  in  studies  of  nuclear 
magnetic  resonance.  Alekseyev  and  Sobelman  [12] 
showed  theoretically  that  collisional  energy  transfer 
between  overlapping  Raman  lines  would  give  rise  to  a 
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narrowing  of  the  Raman  signature  with  increasing  pres¬ 
sure.  The  pressure  dependence  of  the  spontaneous 
Raman  spectrum  of  N2  was  investigated  by  Mikhailov 
[13]  and  by  May  et  al.  [14],  The  latter  high  resolution 
experiments  showed  that  the  fundamental  vibrational 
band  of  N2  undergoes  a  continuous  contraction  in 
width  as  the  pressure  increases  from  about  30  to 
roughly  350  atmospheres.  This  result  was  confirmed 
by  Wang  and  Wright  [15],  Collisional  narrowing  is  also 
an  explanation  for  the  fact  that  the  Raman  spectrum 
of  liquid  N2  is  even  sharper  than  it  is  in  the  gaseous 
state  [16], 

In  the  gas  phase,  the  broadening  of  isolated  N2  Q- 
branches  arises  mainly  from  lifetime  limiting,  rota- 
tionally  inelastic  collisions,  and  as  the  pressure  is  in¬ 
creased,  this  collisional  energy  transfer  becomes  in¬ 
creasingly  frequent.  Eventually  molecules  will  under¬ 
go  so  many  changes  of  rotational  quantum  number 
during  a  scattering  process  that  they  will  tend  to  radi¬ 
ate  at  an  average  frequency  rather  than  that  of  individ¬ 
ual  lines.  The  intensity  distribution  will  then  collapse 
toward  a  center  of  gravity  determined  by  the  frequen¬ 
cy  of  the  transition  associated  with  the  largest  Boltz¬ 
mann  population.  The  narrowing  effect  can  also  be 
thought  to  arise  from  a  reduced  rotationally  inelastic 
collision  frequency  that  is  caused  by  the  coalescence 
of  lines  [20]. 

3.  Theory  of  collisional  narrowing  in  CARS 

In  CARS,  incident  laser  beams  at  frequencies  Wj 
and  u>2  (termed  the  pump  and  Stokes  sources,  re¬ 
spectively)  interact  through  the  third-order,  nonlinear 
electric  susceptibility,  X3. t0  produce  polarization  at 
Wj  =  2wj  -  co2.  If  the  frequency  difference  wj  -«2 
coincides  with  the  frequency  of  a  Raman-active  vibra¬ 
tional  mode,  a  coherent  oscillation  is  excited  in  the 
molecule,  and  the  CARS  polarization  appears  reso¬ 
nantly  enhanced  as  the  upper  side  band  of  a  polari¬ 
zability  modulated  at  this  difference  frequency. 

The  interaction  energy  governing  the  nuclear  re¬ 
sponse  may  be  expressed  as  -  -j  aE2  where  a  is  mo¬ 
lecular  polarizability  and  E  is  electric  field  strength. 
For  two  monochromatic  sources 

^i(0  =  t(^i  e-iw,r  +c.c.), 

£2(0=l(>l2e"iW,'+cc.),  (1) 


E2  will  contain  the  terms 

2*lA2e  1  2  +  c.c.,  (2) 

which  will  be  responsible  for  the  resonant  nuclear  con¬ 
tribution  to  the  polarization.  In  terms  of  the  molecu¬ 
lar  density  operator,  p(t),  the  third-order  polarization 
•Pj(f)  may  be  expressed  as: 

P3(t)  =  NTHp<*)E]U),  (3) 

where  JV  is  the  number  density  of  Raman-active  mole¬ 
cules.  Fourier-transforming  (3)  gives: 

P3(oj)  =  N  J Tr  [p(cj-co’)  a]  dto'.  (4) 

where  the  ~  denotes  Fourier  transform.  The  equation 
of  motion  for  the  density  matrix  is  now  solved  with 
the  interaction  energy  =  and  with  a  col¬ 

lision  integral  representing  linewidth  effects  [12].  This 
equation  is: 

P  =  -  jr  [H0,p\ -  -jr  [h(t), p]  + J(S+pS-p )  Fig)  dg, 

(S) 

where  H0  is  the  unperturbed  molecular  hamiltonian, 

S  is  the  collisional  scattering  matrix,  and  F(g)  is  the 
rate  of  g-type  collisions.  The  collision  integral  assumes 
the  validity  of  the  impact  approximation  and  thus  this 
analysis  is  limited  to  the  binary  collision  regime.  Eq. 
(5)  is  now  solved  in  the  usual  manner  [21  ]  by  a  per¬ 
turbation  expansion  in  powers  of  the  interaction  ener¬ 
gy  h  it).  If  the  resulting  equation  of  change  for  the 
first  order  correction  to  (5)  is  Fourier-transformed, 
and  the  operations  of  eq.  (4)  carried  out,  then  X3  is 
the  coefficient  ofjA^Aj  in  the  expression  for  the 
Fourier  components  of  the  polarization.  The  resulting 
expression  for  X3  is: 

Xj^YStt.ApflG]^,  (6) 

t  n  s 

where  the  t  and  s  indices  denote  the  collection  of 
quantum  numbers  (vibrational-rotational)  belonging 
to  particular  2-branch  transitions,  A pf^  is  the  unper¬ 
turbed  fractional  population  difference  for  transition 
s,  and  G,s  is  the  "G  matrix”  that  is  given  by: 

Gts  =  Kwi-w2~  wr)5» 

+  (Tf/2  -  iAf)6„  +  7W(1— 6W).  (7) 
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Here  cc,.  f'  and  A,  are  t he  vacuum  frequency,  isolat¬ 
ed  linewidth,  and  shift,  respectively,  of  transition  r; 

Srs  is  the  Dirac  delta  function;  and  yn  is  an  off-diago¬ 
nal  linewidth  parameter.  1',,  Af  and  yts  are  in  turn  re¬ 
lated  to  the  scattering  matrix  S  by  the  relationships: 

T  l'r  -  lAr  =/(  1  -  S*tiStil2)F(g)  dg  , 

y»  =  -ft;,s,Ss2, (») 

where  the  subscripts  1  and  2  refer  to  the  initial  and 
final  quantum  states,  respectively,  of  the  transitions  t 
and  s.  Vibrationally  inelastic  collisions  (Auv^O)  can 
be  ignored  in  N2  at  room  temperature,  so  that  to  a 
good  approximation 

~^vj,vj'^Vu'  ■ 

If  the  influence  of  vibration-rotation  interaction  and 
vibrational  dephasing  collisions  are  also  neglected,  so 
that 

Vr  ~ 

one  gets  the  result  that  for  Q-branches 

ir,=/(l-|S„|2)F(s)d*, 

yt,m-f^J2F(g)ig,  Ar=0-  00 

where  the  r  and  s  now  refer  to  the  rotational  quan¬ 
tum  number  of  the  transition,  i.e.,  Q(t).  Thus,  in  the 
absence  of  vibrational  effects,  the  shifts  vanish  and 
the  y[s  are  real.  In  addition  to  the  known  inefficiency 
of  N2  vibrational  relaxation,  the  fact  that  the  shifts 
in  N2  are  relatively  small  [14,15]  and  that  the  isolated 
linewidths  display  essentially  no  vibrational  depen¬ 
dence  nor  evidence  of  vibrational  dephasing  [10] lend 
support  to  these  approximations.  Thus,  our  calcula¬ 
tions  have  been  made  neglecting  vibrational  depen¬ 
dences  of  linewidth  parameters.  One  consequence  of 
eq.  (10)  is  that  the  yts  will  only  couple  Q-branches 
belonging  to  the  same  vibrational  band. 

The  unitarity  principle  of  the  S-matrix 

ElS„l2  =  l,  (12) 

S 

leads  to  the  following  relationship  between  T,  and 
the7„ 


E7„=  r,/-- 


Further,  the  yts  will  obey  detailed  balance  lelation- 
ships: 


(14) 


where  p(0'  denotes  the  normalized  rotational  Boltz¬ 
mann  population. 

If  the7,s  are  ignored  in  the  (/'-matrix  equation  (7), 
then 


|z|  - 1  =  " 

1  '»  i(u,-u2  -w,)  +  r,/2 

and  eq.  (6)  becomes 


2  N  sp 


n,2Ap<0' 


2(Aco12  -  co,)  -  if, 


(15) 


1 1  hi 


This  conventional  sum  representation  is  therefore 
strictly  valid  only  for  isolated  lines. 

For  moderate  pressures  where  the  line  overlap  is 
small,  one  can  expand  G,~'  as  follows' 

[G]  -  ]G}?>]  ->  6„  ~<G«?'r 1  7„(C®>  r  ■ 1  ( 1  6,j), 

(17) 

where 

G®)  =  i(«1~w2  co,)  +  I',/2. 

Substituting  into  eq.  (6 )  gives 

_2N  V  atAp(,0) 

*3  =  n  ,  2Aw,-iT, 

i(V^  a,  yx 

*  i  iAw, +T,/2  j  iAcOj  +  Fj/2'  (18) 


This  expression  is  valid  in  the  limit  of  small  overlap¬ 
ping.  In  general,  however,  the  complex  G-matrix  has 
to  be  numerically  inverted. 

Eq.  (11)  shows  that  the  off-diagonal  linewidth  pa¬ 
rameters  yts  govern  the  rates  of  colhsional  energy 
transfer  between  lines  Q(r)  and  Q(s)  While  the  F, 
and  A,  determine  the  width  and  shift  of  isolated  lines, 
the  off-diagonal  linewidth  parameters  control  the  col¬ 
lapse  or  coalescence  of  overlapping  lines  which  occurs 
in  collisional  narrowing.  If  vibrationally  inelastic  and 
dephasing  collisions  are  slow,  then  the  yrs  are  equiv¬ 
alent  to  minus  the  first-order  rate  constants  for  rota- 
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tional  cross-relaxation.  For  N2 ,  which  consists  of 
ortho-  and  para-modifications,  the  selection  rule  on 
rotational  quantum  number  for  these  processes  will 
be  ±  2,  ±  4 . Unless  the  selection  rule  ±  2  is  assum¬ 

ed,  however,  the  number  of  unknown  linewidth  pa¬ 
rameters  exceeds  the  number  of  relationships  between 
these  (eqs.  (13),  ( 14),  and  adjustable  parameters  will 
have  to  be  introduced  if  the  yt  cannot  be  calculated 
on  an  a  priori  basis. 

In  the  calculations  to  follow,  a  provisional  model 
for  the  yn  has  been  employed  in  which  the  rate  con¬ 
stant  for  the  transition./-*  /'  is  assumed  to  be  propor¬ 
tional  to  the  exponential  of  the  rotational  energy  de¬ 
fect  [22.23] .  That  is,  for  downward  transitions 
( J'<J)yj  j  is  given  by: 

Uj  =dj.Aj  exp(-Cj\Atjj  \/kT),  (19) 

where  dj.  is  the  rotational  degeneracy  of./',  Aj  and  Cj 
are  parameters,  and  AEjj.  is  the  rotational  energy  de¬ 
fect.  For  upward  (/'>/)  transitions,  detailed  balance 
gives 

yjj.  =  dr  Aj 

X  exp{-AEjj.lkT)exp(-Cj\AEjj.\/kT).  (20) 

This  general  functional  form  for  the  rotational  relax¬ 
ation  rate  matrix  describes  the  relaxation  patterns  of 
hydrogen  halids  fairly  well  [22  [ ,  and  has  an  informa¬ 
tion  theoretic  foundation  [23].  If  a  set  of  Cj  is  assumed, 
then  the  values  of  the  A  j  parameters  follow  from  eqs. 
(13)  and  ( 1 4). 

Values  of  Cj  >  1  correspond  to  a  diagonally  domi¬ 
nant  rate  matrix;  that  is,  one  in  which  values  of  7fJ 
close  to  the  diagonal  dominate  (A J  =  ±  2).  The  other 
limit,  Cj  <  1 ,  corresponds  to  a  situation  in  which  all 
rotational  levels  are  strongly  coupled  by  collisions, 
and  all  values  of  A  /  (even)  occur  with  comparable 
probability.  In  the  calculations  to  follow,  the  isolated 
linewidths  reported  in  [10]  were  scaled  linearly  with 
pressure,  and  shifts  were  not  included. 


4.  Experimental 

CARS  spectra  of  gaseous  nitrogen  were  obtained 
by  means  of  broadband  collinear  CARS  generation  in 
a  high  pressure  vessel.  Bandwidths  (fwhm)  were  meas¬ 
ured  with  a  high  resolution  monochromator  as  a  func¬ 


tion  of  pressure  from  one  to  one  hundred  atmosphere 
for  comparison  with  the  theoretical  predictions.  For 
these  investigations,  the  532  nm  output  from  a  Q- 
switched,  frequency-doubled  neodymium:  YAG  puls¬ 
ed  laser  (Quantu-Ray  Model  DCR)  was  combined  with 
the  broadband  output  of  a  circulating-dye  laser,  end- 
pumped  with  a  portion  of  the  doubled  Nd.YAG  out¬ 
put.  The  532  nm  radiation  was  narrowed  to  a  band- 
widthof  ~  0.4  cut-1  (as  measured  by  a  double  mono¬ 
chromator)  with  an  intracavity  etalon.  The  Stokes 
beam  from  the  dye  laser  was  wavelength  centered  at 
607  nm  by  adjusting  the  dye  concentration;  the  band¬ 
width  was  about  200  cm' 1 .  The  dye  laser  energy  was 
approximately  10  mJ,  in  a  ~8  ns  pulse.  More  detailed 
descitpttons  of  the  experimental  aspects  of  CARS  are 
given  in  [5 1  and  |8I 

Collinear  CARS  was  achieved  by  combining  the 
532  nm  pump  beam  and  the  Stokes  dye  laser  beam 
with  a  dichrotc  mirror  and  focussing  into  the  cylindri¬ 
cal  high  pressure  cell.  1  he  emerging  beams,  including 
the  nascent  CARS  beam,  were  collimated  and  spectral¬ 
ly  dispersed  with  an  extra  dense  flint  prism.  The  pump 
and  Stokes  beams  were  trapped  and  the  CARS  beam 
directed  to  the  double  monochromator  for  scans  of 
the  bandwidth.  At  the  lugher  pressures, p  >20  atm, 
the  pump  laser  beam  was  attenuated  by  a  factor  of 
100  with  neutral  density  filters.  This  attenuation  pre¬ 
vented  breakdown  and  stimulated  Raman  scattering 
(SRS)  from  taking  place.  The  CARS  signals  were  aver¬ 
aged  with  a  PAR  boxcar  averager  in  the  A/B  mode; 
the  B  channel  normalizing  signal  was  obtained  by  split¬ 
ting  off  a  portion  of  the  CARS  beam  before  entering 
the  spectrometer.  The  CARS  spectra  were  scanned  at 
the  slowest  available  rate,  1  cm' ’/min.  The  resolution 
of  the  double  monochromator  was  measured  by  scan¬ 
ning  a  low-pressure  Pen-Ray  xenon  lamp  and  found  to 
be  0.38  cm-1. 


5.  Discussion 

The  overall  instrumental  resolution  was  sufficient¬ 
ly  good  to  give  an  N2  signature  at  one  atmosphere  in 
which  individual  even-/  vibrational  Q-branch  transi¬ 
tions  could  be  observed  (fig.  1 ).  At  higher  pressures,  it 
was  no  longer  possible  to  observe  these  Q-branches, 
and  the  experimental  spectrum  tended  to  assume  a 
more  symmetrical  shape,  with  an  overall  width  slightly 
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Fig.  t.  CARS  signatures  of  300  K  N2  at  one  and  100  atm, 

- experiment, - theory.  The  isolated  linewidths 

reported  in  [10|  were  scaled  linearly  with  pressure  in  the  cal¬ 
culations. 

less  than  the  one  atmosphere  value.  Fig.  I  also  shows  a 
typical  100  atm  spectrum  with  a  width  of  about  1 .5 
cm Also  shown  in  fig.  1  are  theoretical  predictions 
based  on  eqs.  (6),  and  (19),  (20). 

In  generating  the  theoretical  curves  shown  in  fig.  1 , 


it  was  found  that  calculations  based  on  the  selection 
rule  \AJ\  =  2  resulted  in  a  predicted  width  of  about 
1  cm-1  at  100  atm,  a  value  about  337f  lower  than  that 
observed.  The  good  agreement  exhibited  at  100  atm 
resulted  from  assuming  that  the  C j  parameters  (eqs. 
(19).  (20)  all  had  a  value  of  1 .75.  In  fact,  collisional 
narrowing  was  found  to  play  a  role  in  obtaining  the 
good  agreement  shown  at  one  atmosphere.  In  the  iso¬ 
lated  line  approximation,  the  O-branches  at  one  at¬ 
mosphere  stood  out  too  strongly,  and  better  agree¬ 
ment  was  obtained  by  taking  narrowing  into  account. 
It  is  apparent  from  fig.  1  that  the  spectrum  maximum 
undergoes  a  shift  of  slightly  more  than  1  cm  1  over 
this  pressure  range.  The  theoretical  curve,  which  does 
not  include  line  shifts,  was  translated  on  the  frequen¬ 
cy  axis  until  the  peaks  coincided. 

The  very  significant  role  that  collisional  narrowing 
plays  in  the  high  pressure  experimental  signatures  can 
be  appreciated  by  considering  fig.  2.  where  theoretical 
predictions  based  on  the  isolated  line  approximation 
(dotted  line)  and  overlapping  line  theory  (solid  line) 
are  compared.  The  width  of  an  isolated  line  signature 
will  be  determined  by  the  homogeneous  linewidths  if 
these  quantities  are  larger  than  a  spectral  bandwidth 
determined  by  the  rotational  population  distribution 
and  the  spacing  between  Q-branch  transitions.  At 
100  atm  the  homogeneous  linewidths  will  have  values 
of  about  10cm~‘,and  the  predicted  isolated  line  sig¬ 
nature  has  an  overall  bandwidth  of  about  this  value. 


FREOJFNCY  _  cm  -  ' 

Fig.  2.  Comparison  cf  calculated  100  otm  signature  with  and 

without  collisional  nirttwin?, - with  narow-n*. 

without  narrowine. 
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Collisional  narrowing  caused  by  line  overlap  thus  re¬ 
sults  in  a  contraction  of  the  100  atm  N2  CARS  signa¬ 
ture  of  about  85%. 

The  pressure  dependence  of  the  CARS  signature 
bandwidth  (fwhm)  is  shown  in  fig.  3.  The  experimen¬ 
tal  points  represent  data  which  were  collected  over  a 
period  of  two  weeks:  a  particular  symbol  denotes  data 
from  one  day’s  runs.  The  scatter  in  the  data  most  like¬ 
ly  result  from  the  extreme  alignment  sensitivity  which 
arises  from  working  at  the  resolution  limit  of  the  dou¬ 
ble  monochromator,  coupled  with  the  very  long  opti¬ 
cal  lever  arms  which  direct  the  CARS  beam  to  the 
monochromator.  Additionally,  the  assignment  of 
bandwidth  to  each  spectrum  is  prone  to  subjective  er¬ 
ror  in  reading  of  peak  height  and  full  width  at  half¬ 
height  .  This  subjective  error  is  random,  and  is  estimat¬ 
ed  to  be  at  least  +0.1  cm-1.  This  seems  reasonable  in 
view  of  the  fact  that  the  spectrometer  resolution  is 
~0.4cm~1.  With  these  remarks  in  mind,  the  experi¬ 
mental  data  suggest  that  the  bandwidth  is  essentially 
constant,  after  the  narrowing  of  the  Q-branch  struc¬ 
ture  which  occurs  from  one  to  about  10  atm. 

Theoretical  calculations  for  two  values  of  the  C  pa¬ 
rameter  (suppressing  the  J  subscript)  are  shown  in 


Fig.  3.  Pressure  dependence  of  the  CARS  signature  width 
(fWhm).  The  solid  lines  give  the  theoretical  predictions,  and 
the  various  o,  a,  and  □  symbols  denote  the  results  of  different 
series  of  experimental  scans  taken  over  a  period  of  several 

weeks. 


fig.  3.  with  the  C  =  1 .75  curve  giving  slightly  better 
overall  agreement.  There  appears  to  be  a  slight  discrep¬ 
ancy  between  experiment  and  theory  for  low  pressures 
up  to  approximately  10  atm,  where  the  theory  exhib¬ 
its  a  pronounced  local  minimum.  Otherwise,  the  agree¬ 
ment  is  quite  good.  The  calculated  profiles  have  a  com¬ 
plicated  dependence  on  C.  As  has  been  stated,  use  of 
a  diagonally  dominant  (AJ  =  ±2)  rotational  relaxation 
rate  matrix  leads  to  predicted  profiles  that  are  nar¬ 
rower  than  those  observed.  As  C  is  reduced  in  magni¬ 
tude,  the  predicted  widths  increase;  for  C  =  1 ,  the  100 
atm  value  is  1 .75  cm-1.  Further  reductions  in  the  val¬ 
ue  of  C,  however,  result  in  a  contraction  of  the  band- 
widths.  For  C  =  0.75,  good  agreement  with  experi¬ 
ment  is  achieved  at  10  and  100  atm,  but  the  inter¬ 
mediate  widths  are  much  too  large.  In  the  limit  of 
strongly  coupled  rotational  states  (C-+  0),  the  predict¬ 
ed  bandwidths  are  again  much  narrower  than  those 
observed  experimentally.  Thus,  these  experimental  re¬ 
sults  are  most  consistent  with  a  relaxation  rate  matrix 
that  lies  between  the  limits  of  diagonal  dominance 
and  strong  coupling.  In  the  best  fit  calculation  AJ  = 

±  2  steps  are  predominant,  but  small  contributions 
from  AJ  =  ±4,  ±  6,...  are  important  as  well. 

In  fig.  3,  the  theoretical  width  reaches  a  shallow 
maximum  between  50  and  75  atm,  and  then  resumes 
a  continuous  narrowing  with  increasing  pressure.  Cal¬ 
culations  for  pressures  higher  than  those  shown  in 
fig.  3  show  that  the  bandwidth  should  undergo  nar¬ 
rowing  out  to  about  300  atm.  This  pressure  depen¬ 
dence  is  somewhat  different  than  that  reported  for 
the  spontaneous  Raman  spectrum  of  N2  [14,15], 
where  a  maximum  at  about  30  atm  is  followed  by  a 
continuous  narrowing  out  to  about  350  atm.  This  dif¬ 
ference  is  not  surprising  due  to  the  nonlinear  nature 
of  CARS. 

CARS  could  be  employed  for  species  concentration 
measurements  by  making  integrated  intensity  meas¬ 
urements.  At  low  pressures,  where  line  overlap  is  small, 
the  integrated  intensity  for  a  majority  species  is  pre¬ 
dicted  to  vary  linearly  with  pressure.  This  assumes,  of 
course,  that  the  pressure  is  not  so  low  that  Dicke  nar¬ 
rowing  or  Doppler  broadening  effects  are  important. 
For  higher  pressures,  the  integrated  CARS  signal  is 
predicted  to  vary  much  more  rapidly  with  pressure. 

As  the  pressure  increases  from  50  to  300  atm,  the  pre¬ 
dicted  dependence  goes  from  quadratic  to  nearly  cu¬ 
bic,  and  then  reverts  at  the  higher  range  to  approxi¬ 
mately  a  quadratic  dependence.  These  are  the  depen- 
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dences  expected  for  a  single  transition  with  a  width 
that  either  is  independent  of  pressure  (quadratic),  ot 
varies  inversely  with  pressure  (cubic).  This  is  consis¬ 
tent  with  the  picture  of  collisional  narrowing  causing 
a  coalescence  of  lines  into  a  single  effective  transition 
radiating  at  an  average  frequency.  While  this  predicted 
quadratic-cubic  pressure  dependence  has  favorable  im¬ 
plications  for  CARS  species  concentration  measure¬ 
ments,  interference  from  the  background  nonlinear 
susceptibility  will  make  this  variation  considerably 
less  dramatic  for  a  minority  species. 

Although  the  particular  model  for  the  off-diagonal 
linewidth  elements  employed  here  is  physically  plau¬ 
sible  and  is  commonly  employed  in  studies  of  rota¬ 
tional  relaxation,  it  is  possible  that  other  semi-empiri¬ 
cal  functional  forms  [24,25)  for  the  rotational  relax¬ 
ation  rates  will  also  give  good  results.  It  would  be  most 
desirable  not  to  have  to  rely  on  the  use  of  adjustable 
parameters.  Ultimately,  it  should  be  possible  to  calcu¬ 
late  all  the  linewidth  elements  for  Nj  in  terms  of  inter- 
molecular  potential  parameters  in  much  the  same  way 
that  the  isolated  linewidths  and  shifts  are  calculated 
[20,26 j .  Also,  the  role  of  line  shifts  will  have  to  be  ac¬ 
counted  for,  particularly  if  higher  density  regimes  are 
investigated. 


6.  Summary  and  conclusions 

Investigations  of  the  300  K  N2  CARS  spectrum 
over  the  range  from  one  to  100  atm  have  shown  that 
collisional  narrowing  plays  an  important  role  in  deter¬ 
mining  the  shape  of  the  profile  over  the  entire  pres¬ 
sure  range.  The  pressure  dependence  of  the  N2  CARS 
spectrum  is  somewhat  different  from  that  reported 
for  the  spontaneous  Raman  spectrum.  An  expression 
for  the  third-order,  nonlinear  electric  susceptibility 
has  been  derived  which  is  valid  for  overlapping  spectral 
lines,  and  calculations  based  on  this  generalized  sus¬ 
ceptibility  expression  are  in  good  agreement  with  the 
experimental  results. 

The  authors  wish  to  acknowledge  helpful  discus¬ 
sions  about  CARS  theory  with  Dr.  R.L.St.  Peters,  and 
the  very  capable  assistance  of  Mr.  Normand  Gantick 
with  the  experiments. 
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Abstract.  Coherent  anti-Stokes  Raman  spectroscopy  (CARS)  appeals  vet, 
promising  for  the  remote,  spatially  and  temporally  precise  probing  of  hostile 
combustion  environments  due  to  its  large  signal  conversion  efficiency  and 
coherent  signal  nature  CARS  is  a  wave  mixing  process  in  which  incident  laser 
beams  at  frequencies  u>,  and  u^.  wdh  a  frequency  difference  tuned  to  a  Raman 
resonance  in  the  molecular  species  being  probed,  interact  to  generate  a  coher 
ent  signal  at  frequency  103  =  2<u,  -luj.By  analyzing  the  spectral  distribution  of 
the  CARS  signal,  temperature  measurements  can  be  performed  Speces  con 
centration  measurements  derive  from  the  intensity  of  the  CARS  radiation  or.  in 
certain  cases,  from  its  spectral  shape  CARS  spectra  have  been  recorded  in  a 
variety  of  flames  from  the  major  flame  constituents  and  generally  show  very 
good  agreement  with  computer  synthesized  spectra  Significantly,  CARS  has 
been  successfully  demonstrated  with  both  liquid  and  gaseous  fuels  in  the 
,.,'imary  zone  and  exhaust  of  practical  combustors  Both  thermometry  and 
species  concentration  measurements  have  been  performed  High  pressure 
effects  on  CARS  spectra  have  also  been  examined 
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I.  INTRODl  CTION 

I  aver  spectroscopy  is  playing  an  increasing  role  in  the  diagnostic 
probing  of  the  hostile,  but  easily  perturbed,  env  ironments  character¬ 
istic  of  combustion  processes.  Physical  probes  can  seriously  perturb 
the  (lame  properties  they  seek  to  measure,  are  often  limited  in  their 
spatial  resolution  and  temporal  response,  and  may  not  survive  at 
high  temperatures  and  pressures.  Spatially  precise  laser  spectros¬ 
copy.  on  the  other  hand,  has  the  potential  for  the  remote,  nonper¬ 
turbing.  in  .uni  measurement  of  temperature  and  species  concentra¬ 
tions  in  combustion  processes.  Laser  techniques  are  also  capable  of 
high  temporal  resolution.  Three  approaches  which  are  spatially  and 
temporally  precise  have  received  much  attention  in  the  last  several 
years:  spontaneous  Raman  scattering,  coherent  anti-Stokes  Raman 
spectroscopy  (CARS),  and  laser-induced  fluorescence.1  '  Raman 
scattering  and  CARS  are  best  suited  to  thermometry  and  major 
species  concent  rat  ion  measurements,  while  laser-induced  fluorescence 
is  applicable  to  measurement  of  flame  radicals  at  trace  levels.  Raman 
scattering  and  laser-induced  fluorescence  are  incoherent  scattering 
processes  in  which  the  generated  signal  is  dispersed  into  4rrsr  CARS 
is  a  coherent  wave-mixing  process  in  which  the  signal  emerges  as  a 
laser-like  beam  in  a  precise  direction. 

t  |M<M  Ktetvcd  Match  2.  WK  l :  accepted  tor  publication  March  12. 1%I. 
received  b\  Managing  fr'ditor  March  19  I9KI 
c  I9KI  Sik-icH  of  Photo-Optical  Instrumentation  f-.npincers 


Spontaneous  Raman  scattering  has  been  widely  investigated  and 
is  well  understood  It  has  been  utilized  in  a  number  ol  practical 
situations,  but  its  application  is  often  restricted  to  certain  tuel- 
stoichiometries,  cycles,  and  locations,  line  to  its  inherent  weakness, 
it  isgenerallylimited  to  application  in  relatively  "clean"  flames  W  uh 
increases  in  fuel  droplet  fragment  and  soot  concentrations,  lasei- 
induced  interferences  such  as  fluorescence  and  incandescence  can 
exceed  and  mask  detection  of  the  Raman  signals  1  With  greater 
emphasis  being  directed  toward  less  clean,  alternate  fuels,  stronger 
diagnostic  techniques  are  required  which  can  operate  over  a  broad 
range  of  fuel  specifications,  stoichiometries,  and  combustion 
approaches.  Both  CARS  and  laser-induced  fluorescence  appear  to 
possess  this  capability.1  but  only  CARS  has  been  demonstrated  to 
date  in  practical  combustion  environments. ''  '*  The  two  techniques 
have  complementary  capabilities  and  have  been  under  development 
in  our  and  other  laboratories  for  several  years.  In  this  paper,  atten¬ 
tion  will  be  restricted  to  CARS.  The  next  section  of  the  paper 
summarizes  the  theory  and  application  of  CARS  for  combustion 
diagnostics.  Succeeding  sections  survey  the  use  of  CARS  for  tempera¬ 
ture  and  species  concentration  measurements  in  a  variety  of  flame 
and  combustion  systems. 

2.  CARS 

The  theory  and  application  of  CARS  are  well  explained  in  several 
reviews  which  have  appeared  recently.1  I0'l!  As  illustrated  in  Fig.  1. 
incident  laser  beams  at  frequencies  and  a*;  (termed  the  pump  and 
Stokes  respectively)  interact  through  the  third  order  nonlinear  sus¬ 
ceptibility  of  the  medium.  x|||||(— «uv  <U|.  «*,.  -  <uj)  .  to  generate 
polarization  and  coherent  radiation  at  frequency  <u,  =  2a>,  -  a.. 
When  the  frequency  difference  (ai,  —  ut^)  is  close  to  the  frequency  , 
of  a  Raman  resonance  of  a  certain  species,  the  magnitude  of  the 
radiation  at  aij  is  resonantly  enhanced  and  results  in  a  signature 
uniquely  characteristic  of  that  species.  In  isotropic  media  such  as 
gases,  the  third  order  susceptibility  is  actually  the  lowest  order  non¬ 
linearity  exhibited,  i.e..  due  to  symmetry  considerations,  second 
order  effects  are  nonexistent.  The  third  order  nonlinear  susceptibility 
tensor  is  of  fourth  rank.  The  subscripts  denote  the  polarization 
orientations  of  the  four  Fields  in  the  order  listed  parenthetically.  In 
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Fig.  1 .  Coherent  anti-StokM  Raman  apactroacopy  (CARS). 


isotropic  media,  the  tensor  must  be  invariant  to  all  spatial  symmetry 
transformations  and  the  81  tensor  elements  reduce  to  just  three 
independent  components.  In  CARS,  which  is  frequency  degenerate, 
there  are  only  two  independent  elements. 

In  all  of  the  experiments  to  be  described  in  this  paper,  vibrational 
Raman  resonances  are  probed,  and,  because  the  vibrational  modes  of 
the  molecules  are  mainly  symmetric  (small  depolarization),  only 
Q-branch  transitions  with  no  change  in  rotational  quantum  number 
are  important.  The  CARS  spectrum  is  governed  by  the  squared 
modulus  of  the  third-order  electric  susceptibility,  which  in  the  iso¬ 
lated  line  approximation  (small  overlap  between  neighboring  transi¬ 
tions)  may  be  expressed  as: 


2(<U|  —  iu 2  —  Wj)  —  i  Tj 


+  xNR 


(I) 


where  the  j  summation  is  over  all  Q-branch  transitions:  N  is  the 
number  density  of  Raman-active  molecules;  aj,  Ap:(0),  a>j,  and  Tj  are 
the  polarizability  matrix  element,  Boltzman  population  difference, 
frequency,  and  isolated  linewidth,  respectively,  of  transition  Q(j); 
and  the  xNR  term  is  an  essentially  dispersionless,  background  contri¬ 
bution  arising  from  distortion  of  electron  orbits,  and.  in  some  cases, 
from  remote  vibrational  resonances.  The  tensor  subscripts  denoting 
polarization  orientation  of  the  electric  field  have  been  dropped. 

The  susceptibility  thus  consists  of  resonant  components  from 
Raman-active  transitions  in  the  species  of  interest  and  a  nonreso¬ 
nant,  mainly  electronic,  contribution  from  all  of  the  molecular  con- 
s’ituents  present.  The  resonant  contribution  to  the  CARS  signal  is 
proportional  to  the  square  of  the  active  molecule  number  density  and 
contains  information  about  the  gas  temperature  through  the  Boltz¬ 
mann  population  factors.  For  very  low  concentrations,  the  desired 


resonant  signal  merges  into  a  baseline  level  due  to  the  nonresonant  | 
background  susceptibility.  When  the  modulation  of  this  dispersion-  * 
less  background  becomes  undetectable,  the  trace  species  is  nominally 
no  longer  measureable.  At  one  time,  this  was  perceived  to  be  a  major  r 
limitation  to  CARS  diagnostics.  However,  the  resonant  and  nonre-  f 
sonant  terms  have  different  tensor  properties.  By  proper  orientation 
of  the  laser  field  and  CARS  detection  polarizations,  the  nonresonant 
electronic  contributions  can  be  cancelled,  permitting  measurements 
of  lower  concentrations  if  the  signal  level  is  adequate.  In  certain  1 
concentration  ranges,  the  presence  of  the  nonresonant  susceptibility 
can  actually  be  used  to  advantage  As  long  as  the  background  modu¬ 
lation  is  detectable,  concentration  measurements  can  actually  be 
made  from  the  shape  of  the  CARS  spectrum,  a  unique  feature  of 
CARS  spectroscopy.  Both  of  these  aspects  will  be  subsequently 
demonstrated. 

Measurements  of  medium  properties  are  performed  from  the 
shape  of  the  spectral  signature  and  or  intensity  of  the  CARS  radia¬ 
tion  The  CARS  spectrum  can  be  generated  in  either  of  two  ways,  as 
seen  in  Fig.  !.  The  conventional  approach  is  to  employ  a  narrow  band 
Stokes  source  which  is  scanned  to  generate  the  CARS  spectrum 
sequentially.  For  nonstationary  and  turbulent  combustion  diagnos¬ 
tics.  this  approach  is  not  appropriate  due  to  the  nonlinear  depen¬ 
dence  of  CARS  on  temperature  and  density  Generating  the  spectrum 
piecewise  in  the  presence  of  large  density  and  temperature  fluctua¬ 
tions  leads  to  signatures  weighted  toward  the  high  density,  low 
temperature  excursions,  and  true  medium  averages  cannot  be  \ 
obtained.  The  alternate  approach"  is  to  employ  a  broadband  Stokes 
source.  This  leads  to  weaker  signals,  but  the  entire  CARS  spectrum  is 
generated  with  each  pulse,  permitting  instantaneous  measurements 
of  medium  properties.  For  efficient  signal  generation,  the  incident  ; 
beams  must  be  so  aligned  that  the  three-wave  mixing  process  is 
properly  phased.  In  gases,  phase  matching  occurs  when  the  input  laser 
beams  are  essentially  aligned  parallel  to  or  collinear  with  each  other  j 

In  many  diagnostic  circumstances,  however,  this  collinear  phase  I 

matching  leads  to  poor  spatial  resolution  because  the  CARS  radia¬ 
tion  undergoes  an  integrative  growth  process  This  difficulty  can  be 
circumvented  by  employing  crossed  beam  phase  matching,  such  as  | 
BOXCARS.14  or  a  variation  thereof. "  "  In  these  approaches.  Fig.  1. 
the  pump  beam  and  Stokes  beam  are  crossed  at  a  point  to  generate  the 
CARS  signal.  CARS  generation  occurs  only  where  the  beams  inter¬ 
sect  and  very  high  spatial  precision  is  possible.  The  intersecting  laser 
beams  need  not  be  coplanar1*  and  experimental  simplifications  can  be 
gained  using  three-dimensional  phase  matching 

CARS  spectra  are  more  complicated  than  spontaneous  Raman 
spectra,  which  are  simply  an  incoherent  addition  of  transitions 
CARS  spectra  can  exhibit  constructive  and  destructive  interference 
effects.  For  most  molecules  of  combustion  interest,  these  effects  can 
be  readily  handled  numerically  since  the  theory  of  CARS  generation 
is  well  understood.  At  United  Technologies  Research  Center 
(UTRC),  CARS  computer  codes  have  been  developed  and  validated 
experimentally  for  the  diatomic  molecules  N2.  H,.  CO.  and  Oj1*  and 
the  triatomics  H2OJ0  and  C02.  Computer  codes  are  extremely  useful 
for  studying  the  parametric  behavior  of  CARS  spectra  and  for  data 
reduction. 

Although  CARS  has  no  threshold  per  se  and  can  be  generated 
with  cw  laser  sources,  high  intensity  pulsed  laser  sources  are  required 
for  most  gas  phase  and  flame  diagnostics  to  generate  acceptable 
signal-to-interference  ratios,  particularly  with  broadband  generation  | 
and  detection. 1  In  the  CARS  work  to  be  reported  here,  a  frequency-  - 
doubled  neodymium  laser  at  5320  A  provides  the  pump  beam  and 
also  drives  the  broadband  Stokes  dye  laser.  Crossed-beam  phase  ■ 
matching  (BOXCARS),  either  planar  or  three  dimensional,  is  used  | 
to  ensure  good  spatial  precision.  The  CARS  signatures  are  detected  ■ 
with  an  optical  multichannel  analyzer  (OM  A)  which  permits  capture 
oftheentireCARSspectruminasinglepulse.  In  laminar  flames  and  f 
situations  where  fluctuation  magnitudes  are  small,  the  CARS  spec-  | 
trum  can  be  averaged  on  the  OMA  or  scanned  with  a  mono¬ 
chromator  using  a  BOXCAR  averager.  Greater  detail  about  the 
apparatus  and  experimental  procedures  may  be  found  in  Refs.  2. 9.  * 

and  18.  I 
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3.  THERMOMETRY 

Temperature  measurements  derive  from  the  spectral  shape  of  the 
CARS  signature.  Thermometry  is  therefore  more  easily  performed 
than  concentration  measurements,  which  generally,  but  not  always, 
require  determination  of  absolute  signal  intensity  levels.  Tempera¬ 
ture  information  can  be  extracted  from  any  of  a  number  of  molecular 
constituents.  In  this  section.  CARS  thermometry  utilizing  several 
molecules  in  a  variety  of  measurement  situations  will  be  illustrated. 

3.1.  Nitrogen 

N it  rogen  is  the  dominant  constituent  in  airfed  combustion  processes 
and  is  everywhere  present  in  large  concentrations.  Performing 
temperature  measurements  from  N2  provides  information  on  the 
location  of  the  combustion  heat  release  and.  to  some  degree,  the 
extent  of  chemical  reaction.  In  Fig.  2  is  displayed  the  calculated 
temperature  dependence  of  the  Nj  CARS  spectrum  for  parameters 
corresponding  to  the  experimental  approach  employed  here.  A  de- 


Fig  2.  Temperature  variation  of  tha  CARS  ipactrum  of  N2  for  a  pump 
linawidth  of  0  8  cm  1 .»  Stoke*  bandwidth  of  ISO  cm-1  andal  cm-'  (lit 
width. 


scription  of  the  computer  code  is  summarized  in  Refs.  I9and  21.  In 
the  code,  the  dependence  of  the  Raman  linewidth  on  temperature  and 
rotational  quantum  number  is  taken  into  account.22  At  low  tempera¬ 
tures.  one  sees  the  v  =  0  —  I  band  containing  low  J  value  Q  branch 
transitions,  i.e.,  AJ  =  0  where  J  is  the  rotational  quantum  number. 
As  is  apparent,  the  low  J  Q  branches  are  unresolved.  As  the  tempera¬ 
ture  increases,  the  band  broadens  as  the  rotational  population  distri¬ 
bution  shifts  to  higher  J  values.  At  very  high  temperatures,  individual 
even  J  Q  branches  ranging  fromQ(20)  to  Q(40)  can  be  resolved.  The 
odd  Q  branches,  which  have  a  nuclear  spin  weighting  equal  to  half  of 
the  even  number  branches,  are  reduced  in  intensity  by  about  a  factor 
of  four  and  do  not  stand  out.  At  intermediate  temperatures,  fewer  Q 
branch  transitions  are  resolvable.  For  Q  branches  beyond  Q(40), 
overlap  with  the  v  =  I  —  2  band  transitions  occurs,  giving  rise  to  two 
prominent  peaks  in  the  “hot”  band.  At  lower  spectral  resolution,  e.g., 
~  3  cm-1,  the  fine  structure  shown  in  Fig.  2  is  lost,  but  the  spectra 
still  exhibit  good  temperature  sensitivity. 

The  accuracy  of  CARS  N2  thermometry  has  been  examined  in 
premixed  flat  flames  by  comparison  with  radiation-corrected, 
coated,  fine  wire  thermocouples.2'  The  radiation  corrections  were 
experimentally  calibrated  at  different  flame  temperatures  by  sodium 
line  reversal.  In  Fig.  3,  the  CARS  spectrum  scanned  at  a  resolution  of 
I  cm"'  in  a  21 10  K  flame  is  displayed.  The  dotted  curve  is  a  theoreti¬ 
cal  least  square  fit  to  the  experimental  signature.  The  best  fit  temper¬ 
ature  is  2 104  K,  with  a  standard  deviation  of  9  K.  Similar  agreement 
was  obtained  at  lower  flame  temperatures. 


Frequency  -  cm" 1 


Fig.  3.  Comparison  of  experimental  (solid  line)  and  theoretical  least 
square  fit  (dashed)  N2  CARS  spectra  at  a  measured  temperature  of 
2110K. 


(a)  Single  pulse  (b)  Averaged 


Fig.  4.  Comparison  of  single-pulse  and  averaged  CARS  spectra  from 
flame  N2  in  a  sooting  propans  diffusion  flame. 


The  capability  of  CARS  for  use  in  highly  sooting  flames  has  been 
demonstrated.22  25  Although  incoherent  and  coherent  interferences 
from  C2  can  be  encountered  in  the  N2  CARS  bands  from  a  5320  A 
pump  laser,  they  are  for  the  most  part  suppressible  by  appropriate 
Stokes  laser  bandwidth  selection  and  use  of  a  polarization  filter.  The 
C2  is  produced  by  the  laser  vaporization  of  soot  which  occurs  even 
for  0  ( 1 0— 8  sec)  duration  laser  pulses.  In  Fig.  4  both  a  single  pulse 
( 1 0~ 8  sec)  and  an  averaged  N2  CARS  spectrum  from  a  highly  soot¬ 
ing.  laminar  propane  diffusion  flame  are  displayed  for  a  spatial 
resolution  of  0.3  X  I  mm.  The  spectra  are  interference  free  and  of 
high  quality.  The  single  pulse  CARS  spectrum  is  possible  at  laser 
energies  an  order  of  magnitude  lower  than  those  typically  employed 
for  single  pulse  spontaneous  Raman  scattering.  Furthermore,  in  this 
sooting  flame,  laser  modulated  particulate  incandescence1  would 
exceed  the  Raman  signal  by  several  orders  of  magnitude,  precluding 
successful  measurements.  In  a  turbulent  flame,  these  single  pulse 
temperature  measurements  would  lead  to  determination  of  the 
temperature  probability  distribution  function  (PDF).  Data  such  as 
those  displayed  in  Fig.  4  have  been  employed  to  perform  detailed 
axial  and  radial  temperature  surveys  in  laminar  propane  diffusion 
flames.21  Figure  5  displays  the  radial  variation  of  temperature  in  the 
sooting  flame  with  height  above  the  burner. 

Recently,  the  feasibility  of  CARS  for  measurement  in  practical 
combustion  systems  has  been  demonstrated.  In  tests  at  our  labo¬ 
ratory,9  BOXCARS  thermometry  was  performed  in  two  different, 
liquid-fueled  combustors  housed  in  a  50  cm  dia  combustion  tunnel. 
Delicate  instrumentation  was  housed  in  a  control  room  adjacent  to 
the  burner  test  cell  and  the  CARS  signals  were  piped  out  employing 
20  m  long,  60  m  dia  fiber  optic  guides.  In  Fig.  6  are  shown  CARS 
signatures  of  N2  averaged  for  10  sec  at  two  different  axial  locations  in 
the  primary  zone  of  a  Jet  A  fueled  swirl  burner.  At  the  upstream,  x 
6  cm.  location.  CARS  measurements  were  made  through  the  fuel 
spray  and  the  temperature  was  found  to  be  about  900  K  at  an  overall 
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Fig  5  Radial  temperature  profile*  determined  by  CARS  in  a  laminar 
propane  diffusion  flame. 


Fig.  6  Spatial  variation  of  temperature  from  averaged  CARS  spectra  of 
Nj  in  swirl  burner  with  Jet  A  fuel. 


Fig.  7.  Comparison  of  averaged  and  single  pulse  N2  CARS  spectra  in 
swirl  burner  with  refractory  back  wall  fueled  with  Jet  A  at  an  overall 
equivalenca  ratio  of  0.8. 


equivalence  ratio  of  0.8  At  the  downstream  location,  where  the 
flame  was  very  luminous,  the  temperature  increased  to  1500  K  In 
Fig  7  is  shown  a  comparison  of  a  single  pulse  (10  *  sec)  and  an 
averaged  CARS  spectrum  ( 10  sec)  in  the  Jet  A  fueled  swirl  burner 
fitted  with  a  refractory  back  wall  to  simulate  a  furnace  more  closely 
The  single  pulse,  although  of  slightly  lower  quality,  is  fairly  good  and 
would  permit  credible  measurements.  Measurements  were  also  suc¬ 
cessfully  performed  in  the  exhaust  of  a  JT- 1 2  combustor  burning  Jet 
A.9  Single  pulse  spectra  were  virtually  indiscernible  from  the  aver¬ 
aged  data  CARS  measurements  downstream  in  the  exhaust  were  in 
very  good  agreement,  i.e  .  10  to  50  K,  with  temperatures  determined 
by  aspirating  thermocouple  probes. 

Most  practical  combustor  devices,  of  course,  operate  at  pressures 
considerably  in  excess  of  atmospheric.  It  is  known  from  spontaneous 
Raman  studies,  however,  that  remarkable  changes  in  spectra  can 


occur  w  hen  the  coiiisi.ut-indiiced  Raman  line  b'oadening  becomes 
comparable  to  the  separation  between  the  individual  vibration- 
rotation  transitions  At  this  point,  the  isolated  Raman  linewidths 
continue  to  broaden  linearly  with  pressure,  but  collisional  energy 
transfer  between  adjacent  J-states  causes  adiacent  transitions  m 
coalesce  or  collapse  toward  a  frequency  center  of  grav  ity  1  he  simple 
representation  of  x*’1  (Fq.  ( I))  as  a  sum  of  independently  broadened 
transitions  is  no  longer  valid  in  regions  of  strong  line  overlap  This  is 
referred  to  as  collisional  narrow  ing  and  its  eflect  must  be  included  in 
l  he  diagnostic  interpretation  of  high  pressure  CARS  spectra  In  most 
diatomic  molecules,  collisional  narrowing  becomes  important  above 
a  few  atmospheres  of  pressure,  with  the  effect  becoming  less  pro¬ 
nounced  as  temperature  is  increased 

For  recently  developed  theories  account  tor  the  phenomenon 
quite  well  at  300  K,’J  and  the  agreement  w  it h  data  at  higher  tempera¬ 
tures  also  appears  to  be  good  As  an  illustration.  Fig  X  displavs 
theory  experiment  comparisons  ol  V  CARS  spectra  at  -  I  MX)  K 
and  pressures  of  7  8  and  28.2  atmospheres.  At  7  K  atmospheres  line 


Fig.  8.  Comparison  of  observed  and  predicted  pressure  dependence  of 
Nz  CARS  spectra  at  1630K. 


structure  is  still  apparent,  and  collisional  narrowing  plays  a  very 
small  role.  At  28.2  atmospheres,  the  isolated  Q-branch  linewidths 
have  increased  to  the  point  where  individual  Q-branches  can  no 
longer  be  resolved,  and  collisional  narrowing  is  play  ing  a  significant 
role.  In  the  isolated  line  approximation,  the  width  of  the  fundamen¬ 
tal  band  and  the  hot  band  peak  height  are  both  about  JO  percent 
larger  at  28.2  atm  than  the  values  observed  in  the  experimental 
signature.  With  narrowing  included  in  the  calculation,  the  theory  is 
in  very  good  agreement  with  experiments,  as  can  be  seen  in  Fig.  8 

3.2.  Hydrogen 

H-,.  when  sufficiently  abundant,  is  ideal  for  combustion  thermome¬ 
try  because  of  the  simplicity  of  its  spectrum,  as  seen  in  Fig.  9.:<  The 
line  spectra  are  components  of  the  v  ibrational  Q  branch,  v  =  0  —  I. 
■VJ  =  0.  The  adjacent  rotational  components  are  well  separated  due 
to  the  large  H,  vibration-rotation  interaction,  making  the  spectrum 
particularly  simple  to  interpret.  The  intensity  alteration  between 
even  and  odd  J  lines  again  is  due  to  nuclear  spin  statistics  in  which  the 
odd  lines  have  three  times  the  statistical  weight  of  the  even-numbered 
transitions.  H,  BOXCARS  signatures  such  as  those  shown  in  Fig.  9 
have  been  used  to  measure  temperatures  at  various  locations  in  a  flat 
H,-air  diffusion  flame.-’'  Approximate  temperatures  were  deduced 
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Fig.  9  CARS  spactra  of  H2  in  H2-air  diffusion  flame  at  temperatures 
determined  from  relative  strengths  of  indicated  Q-branch  transitions. 
Frequency  scale  corresponds  to  0  60  cm  1  per  dot. 


3.3.  \\  atcr  vapor 

Water  vapor  is  the  major  product  of  airfed  hydrogen  combustion 
and  often  the  dominant  product  species  of  hydrocarbon-fueled  com¬ 
bustion.  Its  measurement  is  an  important  gauge  of  the  extent  ot 
chemical  reaction  and  of  overall  combustion  efficiency,  f  igure  1 1 
displays  an  H,0  CARS  spectrum  at  I  cm-1  resolution  in  the  posi- 
flame  region  of  an  atmospheric  pressure,  premixed  CH4  air  flame  at 
1700  K.-’°  Alsoshown  in  Fig.  II  is  the  CARS  computer  code  predic- 


Experiment  Theory 

ISO  pulse  average  Slit  =  1.8  cm"1  I' =  0.2  cm"1 


from  the  ratios  of  the  CARS  intensities  of  the  Q(l).  Q(3).  and  Q(5) 
transitions.  Agreement  among  the  temperatures  deduced  from  the 
three  intensity  ratios  was  good.  i.e..  the  standard  deviation  varied 
from  2  to  8  percent  over  the  range  of  900  K  to  2100  K.  Better  accuracy 
would  probably  result  from  using  ratios  of  integrated  line  strengths 
instead  of  the  peak  intensities  The  results  of  temperature  profiling 
the  flat  diffusion  flame  are  shown  in  Fig.  10.  where  temperatures 
deduced  from  Hi  and  0:  CARS  spectra  are  compared  to  measure- 


ince  across  burner,  y  (mm) 


Fig.  10.  Temperature  measurements  in  s  flat  H2-air  diffusion  flame. 
Symbols:  circles,  radiation-corrected  thermocouple:  open  triangles,  H2 
CARS:  end  solid  triangles,  02  CARS.  Dotted  curve  it  locus  of  maximum 
temperature. 


ments  with  a  radiation  corrected  thermocouple.  The  CARS  spec¬ 
trum  of  Oj.  omitting  spectroscopic  details,  is  qualitatively  quite 
similar  to  that  of  N2.  The  temperatures  agree  quite  well  in  cooler 
regions  of  the  flame.  Larger  discrepancies  occur  at  the  higher 
temperatures  where  the  concentrations  are  low  and  signal-to-noise 
ratio  is  poorer.  By  spectrally  integrating  the  spectra  and  compar¬ 
ing  the  integrated  intensity  to  that  generated  simultaneously  from  a 
high  pressure  gas  cell  in  a  parallel  reference  leg,  density  measure¬ 
ments  have  been  made  in  the  flat  diffusion  flame.  These  were  in  fair 
comparison,  i.e..  ~50  percent,  with  H2  concentrations  measured  by 
spontaneous  Raman  scattering. 


Frequency  - 


Fig.  11.  Comparison  of  experimental  and  theoretical  CARS  spectra  of 
water  vapor  in  a  flame  at  1 700  K.  A  pressure-broadened  linewidth  r  -  0.2 
cm  1  was  assumed  for  all  transitions  in  the  calculation. 


tion  for  the  experimental  conditions.  While  only  the  v,  symmetrical 
mode  of  HjO  is  Raman  active,  the  involved  rotational  energy  level 
structure  of  this  asymmetric  rotor  gives  rise  to  a  complex  CARS 
signature.  In  the  calculations,  nearly  one  thousand  vibrational- 
rotational  Q  branch  transitions  are  included  at  flame  temperatures 
The  peaks  seen  in  the  spectrum  arise  from  spectral  overlapping  of  the 
various  transitions. 

Good  agreement  between  experiment  and  theory  has  also  been 
achieved  in  test  cells  of  pure  H20  for  temperatures  in  the  range  300  K 
to  800  K.w  The  predicted  and  observed  spectra  exhibit  a  pronounced 
sensitivity  to  temperature,  making  the  H20  molecule  potentially 
attractive  for  CARS  thermometry.  However,  theory-experiment 
agreement  is  obtained  only  for  relatively  large  Raman  linewidths.10 
and  H2OCARS  will  be  a  reliable  diagnostic  tool  when  more  informa¬ 
tion  about  the  temperature-,  concentration-,  and  Q-branch  depen¬ 
dences  of  these  linewidths  is  obtained.  Work  is  currently  underway  at 
our  laboratory  to  resolve  these  questions. 

3.4.  Carbon  dioxide 

The  other  major  product  of  hydrocarbon  combustion  besides  water 
vapor  is.  of  course,  carbon  dioxide.  In  Fig.  12  is  shown  the  CARS 
spectrum  of  C02  in  the  postflame  rone  of  a  CO  air  flame  at  I520K 
The  experimental  signature  displays  a  fundamental  band  at  a  Raman 
shift  of  about  1388  cm-1,  corresponding  nominally  to  the  transition 
between  the  ground  vibrational  state  and  the  first  symmetric  stretch 
state.  A  number  of  “hot  bands,"  corresponding  to  vibrationally 
excited  initial  states,  also  appear  at  larger  shifts.  Studies  of  CO, 
CARS  spectra  in  different  flames  at  varying  temperatures  show  that 
the  strength  of  the  hot  bands  is  strongly  dependent  on  temperature, 
making  these  features  useful  for  thermometry.  Unlike  the  other 
spectra  presented  earlier,  there  is  little  rotational  smearing  of  the 
vibrational  bands,  a  consequence  of  the  extremely  small  vibration- 
rotation  interaction  in  C02.  Each  vibrational  band  represents  the 
contributions  of  approximately  100  Q-branch  transitions.  Not 
shown  in  Fig.  12  is  a  slightly  weaker  sequence  of  vibrational  bands 
whose  fundamental  occurs  at  a  Raman  shift  of  1 285  cm- 1 .  and  w  hose 
hot  bands  occur  at  smaller  shifts.  The  fact  that  two  distinct  vibra¬ 
tional  sequences  occur  is  a  consequence  of  "Fermi  resonance"; 
because  the  symmetric  stretch  frequency  nearly  coincides  with  twice 
that  of  the  bending  mode,  anharmonic  terms  in  the  intramolecular 
potential  can  cause  a  large  perturbation  of  the  corresponding  eigen- 
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(a)  Experiment 


(b)  Theory 


Frequency  —  cm  1 


Fig  12.  Comparieon  of  experimental  and  theoretical  CAR  Saignaturea  of 
C02  in  a  llama  at  1 520  K. 


states,  and  a  mixing  of  the  symmetric  stretch  and  bending  modes 
occurs. 

Also  shown  in  Fig.  12  is  a  theoretical  calculation  for  T  =  1520  K. 
As  can  be  seen,  the  agreement  with  experiment  is  fairly  good,  with  a 
few  minor  discrepancies.  Similar  agreement  is  obtained  in  other 
flames  at  different  temperatures.  The  calculation  of  vibrational 
energy  levels,  rotational  constants,  and  spontaneous  Raman  cross 
sections  basically  follows  the  treatment  of  Courtoy  2*  It  has  been 
found,  however,  that  better  agreement  is  obtained  if  tabulated 
values27  for  the  rotational  constants  are  employed.  Further  investiga¬ 
tions  relating  to  spectral  constants,  linewidths,  and  possibly  colli- 
sional  narrowing  will  likely  result  in  even  better  agreement  between 
experiment  and  theory. 

4.  CONCENTRATION  MEASUREMENTS 

In  general,  concentration  information  derives  from  the  intensity  of  the 
spectrally  integrated  CARS  signal,2'  assuming  that  the  nonresonant 
susceptibility  contribution  is  negligible  or  has  been  suppressed  in  a 
manner  to  be  detailed  shortly.  As  alluded  to  earlier,  species  detection 
sensitivity  is  limited  for  conventional  CARS  approaches  with  aligned 
polarizations  due  to  the  presence  of  the  background  nonresonant 


Frequency  -cm'1 


Fig.  13  Computer  CARS  spectra  of  CO  at  various  concentrations  at 
1800K 

susceptibility.  This  situation  is  quite  evident  in  the  computer  calcula¬ 
tions  for  CO  at  1 800K  shown  in  Fig.  13.  A  similar  situation  exists  with 
O,.25  As  seen,  the  CARS  resonant  spectrum  at  20  percent  concentration 
rises  well  out  of  the  nonresonant  background,  which  is  only  slightly 
evident,  and  the  hot  band,  v  =  I— 2.  is  quite  prominent  As  the 
concentration  decreases,  the  modulation  and  the  nonresonant  back¬ 
ground  become  more  significant;  at  still  lower  concentrations,  the 
degree  of  modulation  decreases  until  the  signal,  i.e..  the  modulation,  is 
lost  in  the  nonresonant  background  and  the  species  is  no  longer 
detectable.  To  the  extent  that  the  nonresonant  susceptibility  is  known 
or  can  be  approximated,  the  concentration  of  a  particular  species  can 
be  obtained  from  the  shape  of  the  CARS  spectrum,  obviating  the 
requirement  for  absolute  intensity  measurements.  The  concentration 
range  appropriate  to  spectral  curve  fitting  for  concentration  measure¬ 
ments  will  vary  from  molecule  to  molecule  and  with  gas  composition 
and  temperature.  For  molecules  with  closely  spaced  vibrational- 
rotational  transitions  such  as  CO.  O,.  and  N,  at  flame  temperatures,  the 
range  is  approximately  0. 1  to  20  percent"  This  approach  has  been 
quantitatively  verified  in  the  case  of  CO.  as  seen  in  Fig.  14  There,  an 
experimental  CARS  signature  of  CO  in  the  postflame  region  of  a 
premixed  CO-O-,  flat  flame  is  displayed  The  thin  solid  line  is  the 
theoretical  CARS  overlay  at  a  CO  concentration  of  3.6  percent  corre¬ 
sponding  to  the  level  experimentally  measured  with  a  quartz  micro¬ 
probe  sampling  system  and  a  nondispersive  infrared  analyzer.  As  can 
be  seen,  the  agreement  between  the  CARS  and  probe  mea-'iremenis  is 
quite  good.  This  approach  was  also  verified  by  CARS  measurements 


Fig.  14.  Comparieon  of  experimental  and  computed  (eolid  lino)  CARS 
tpodra  In  •  flame  at  a  concentration  of  3.6%.  T  -  210OK. 
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in  cold,  calibrated  gas  mixtures.  The  spectral  fitting  approach  has  been 
employed  to  perform  02  measurements  as  well2’  and  to  map  02  decay 
through  the  flat  H2-air  diffusion  flame  previously  mentioned. 

For  concentration  measurements  at  lower  levels  than  permitted  by 
spectral  fitting,  the  background  nonresonant  susceptibility  can  be 
suppressed  by  appropriately  orienting  the  laser  field  and  CARS- 
detection  polarizations.-"'  w  Number  densities  are  then  obtained  from 
the  strength  of  the  spectrally  integrated  signal  after  appropriate  cali¬ 
bration.  e  g.,  using  an  external  reference  cell.  In  Fig.  15,  background 
susceptibility  suppression  is  demonstrated  for  CO  in  a  flat  CO-air 
diffusion  flame.  Folded,  nonplanar  BOXCARS,1*  was  employed 
together  with  broadband  Stokes  generation  and  optical  multichannel 
spectral  detection.  With  aligned  polarizations,  the  characteristic  modu¬ 
lated  CARS  spectrum  (Fig.  13)  is  obtained.  By  orienting  the  w,  pump 
field  polarizations  and  CARS  analyzer  relative  to  the  horizontally 
polarized  at,  Stokes  field,  as  shown,  the  nonresonant  susceptibility  was 
reduced  by  over  two  orders  of  magnitude  and  the  background  free 
CARS  spectrum  of  CO  was  obtained. 

Unfortunately,  the  polarization  approach  leads  to  a  considerable 
loss  in  signal  intensity  The  resonant  mode  signal  reduction  is  at  least 
a  factor  of  sixteen.  In  addition,  there  is  the  signal  loss  associated  with 
the  suppression  of  the  nonresonant  susceptibility  itself.  In  a  recent 
experimental  study,11  the  CO  detectivity  of  broadband  CARS  with 
and  without  background  suppression  was  compared.  With  crossed- 
beam  phase  matching,  better  concentration  sensitivity,  by  about  a 
factor  of  five,  was  found  by  not  suppressing  the  background.  With 
the  background  suppressed,  the  CARS  signal  level,  i.e.,  number  of 
photons,  limited  the  detection  sensitivity.  With  the  background  pres¬ 
ent,  detectivity  was  limited  by  the  loss  of  modulation  in  the  nonres¬ 
onant  susceptibility  profile.  In  those  circumstances  where  detectivity 
is  not  photon  but  background  limited,  e  g.,  low  spatial  resolution. 

high  pressure . concentration  sensitivity  will  be  enhanced  by 

employing  nonresonant  susceptibility  suppression. 

With  the  background  suppressed,  density  measurements  are 
based  on  the  strength  of  the  spectrally  integrated  CARS  signature. 
This  usually  requires  the  use  of  a  reference  cell  to  normalize  pulse-to- 
pulse  laser  power  fluctuations,  to  serve  as  a  calibrating  standard,  and 
to  account  for  small  optical  misalignments.  Reference  cells  may  not 
be  practical  when  probing  turbulent,  soot-  and  droplet-laden  gases 
because  they  cannot  account  for  effects  such  as  beam  steering,  defocus- 
ing.  and  attenuation.  It  seems  desirable,  then,  to  perform  concentra¬ 
tion  measurements  from  spectral  shapes  whenever  possible.  In 
essence  one  utilizes  the  nonresonant  susceptibility  at  the  measure¬ 
ment  location  as  an  in  situ  reference.  For  unmodulated  spectra,  i.e.. 
when  the  nonresonant  susceptibility  is  suppressed  or  unimportant, 
there  are  other  approaches  to  in  situ  referencing.  These  approaches 
are  receiving  considerable  attention  for  turbulent  combustion  diag¬ 
nostics  because  of  their  potential  for  circumventing  the  foregoing 
problems  and  are  reviewed  in  Refs.  31  and  32. 

Hydrocarbon  fuels  are  Raman  active,  often  with  large  Raman 
sections,  and  CARS  could  in  principle  be  used  to  monitor  fuel 
mixing  and  pyrolysis  processes.  Other  than  observing  CARS  spectra 
of  hydrocarbons  in  flames  or  cells,  little  work  of  a  quantitative  nature 
has  been  reported,  although  this  should  be  a  fruitful  area  for  future 
investigations. 

5.  SUMMARY 

C  A  RS  is  a  laser  spectroscopic  technique  appropriate  to  spatially  and 
temporally  resolved  measurement.!  of  temperature  and  major  species 
concentrations  in  combustion.  It  has  been  generated  from  all  of  the 
dominant  constituents  in  air-fed,  hydrogen-,  and  hydrocarbon- 
fueled  combustion.  Quite  importantly,  CARS  has  been  demonstrat¬ 
ed  to  be  applicable  to  practical  combustion  systems. 
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CARS  thermometry  in  reacting  systems* 

J.  F.  Verdieck,  J.  A.  Shirley,  R.  J.  Hall,  and  A.  C.  Eckbreth 

United  Technologies  Research  Center,  East  Hartford.  Connecticut  06108 

CARS  (Coherent  anti-Stokes  Raman  Spectroscopy)  diagnostic  techniques  provide  accurate  temperature 
measurements  in  environments  ranging  from  controlled  laboratory  flames  to  hostile  systems  such  as  gas 
turbine  combustors  and  internal  combustion  engines.  CARS  is  an  optical  wave  mixing  process  wherein 
incident  laser  beams  at  frequencies  to ,  and  to  „  with  a  frequency  difference  appropriate  to  the  molecular 
species  being  probed,  interact  to  generate  the  CARS  signal  at  <a,  =  Zto,  —  to,.  Through  analysis  of  the 
frequency  distribution  of  the  CARS  spectrum,  temperature  measurements  can  be  obtained.  Concentration 
measurements  derive  from  the  intensity  of  the  CARS  radiation  or,  in  certain  cases,  from  its  spectral  shape. 
CARS  is  a  remote,  nonperturbing  technique  which  determines  temperature  in  a  matter  of  seconds  with  good, 
spatial  resolution  (of  the  order  of  millimeters).  Under  some  circumstances,  CARS  signal  levels  are  intense 
enough  to  permit  "instantaneous,"  single  pulse  temperature  measurements,  i.e„  — 10~'  seconds.  Because 
CARS  is  a  coherent,  nonlinear  optical  technique,  it  offers  decided  advantages  over  conventional  incoherent 
spectroscopic  methods.  These  advantages  derive  both  from  the  strong  CARS  signal  and  its  laser-like  beam 
character  which  furnish  excellent  discrimination  against  interferences  such  as  combustor  luminosity  and 
laser-induced  interferences,  which  are  often  quite  severe.  In  this  paper,  the  fundamentals  of  CARS 
thermometry  from  both  theoretical  and  experimental  points  of  view,  with  emphasis  on  the  nitrogen  molecule, 
will  be  discussed  in  detail.  Areas  covered  will  include  the  accuracy  of  CARS  thermometry,  computer¬ 
generated  CARS  spectra,  linewidth  considerations  and  high  pressure  effects.  Additionally,  the  use  of  other 
molecular  thermometric  probes  for  CARS  diagnostics,  H„  O,,  H,Op  CO,  and  CO„  will  be  shown.  Two 
important  experimental  techniques,  BOXCARS,  a  particular  laser  beam  configuration  which  achieves  good 
spatial  resolution,  and  single-shot  thermometry,  wherein  the  CARS  spectrum  from  a  single  pulse  is  recorded, 
will  be  treated.  The  utilization  of  CARS  thermometry  will  be  illustrated  by  its  application  in  several  different 
flame  systems  and  in  practical  combustor  devices. 


INTRODUCTION 

Temperature  measurements  utilizing  spectroscopic 
techniques  such  as  absorption  and  emission  spectro¬ 
scopy  are  well  established  methods  which  have  served 
well  over  many  decades.  Mere  recently,  Raman 
scattering  and  fluorescence,  using  laser  sources,  have 
been  applied  with  considerable  success.  These  spec¬ 
troscopic  techniques  offer  several  advantages  over 
probe  thermometry  such  as  remote  operation,  non- 
perturbation  of  the  sample,  and  extremely  high 
temperature  capability.  The  major  disadvantage  of 
spectroscopic  thermometry  is  that  the  system  under 
investigation  must  be  "seen",  i.e.,  be  visible 
through  some  type  of  optical  ports. 

This  paper  describes  a  more  recent  thermometric 
method  based  upon  a  coherent,  nonlinear  optical 
technique  named  CARS  (for  coherent  anti-Stokes  Raman 
spectroscopy) .  In  addition  to  the  advantages  listed 
above,  the  CARS  method  offers  several  others;  good 
accuracy  (better  than  IX  demonstrated  with  nitrogen 
CARS  thermometry) ,  rapid  measurement  and  processing 
(typically  less  than  ten  seconds),  and  most  impor¬ 
tantly,  the  capability  for  a  high  degree  of  spatial 
resolution  (better  than  one  millimeter).  Additionally, 
the  exciting  potential  for  performing  virtually 
Instantaneous  (in  10"®  seconds)  temperature  measure¬ 
ments  is  available  with  single  pulse  CARS.  In  this 
paper  the  term  reacting  system  shall  be  loosely 
Interpreted  to  mean  combustion  systems.  The 
extension  of  CARS  measurement  to  other  types  of 
chemically  reacting  systems  is  easily  made.  It 
should  be  emphasized,  that  in  combustion,  the  single 
most  Important  parameter  to  measure  is  the  temperature. 
Given  the  initial  conditions,  and  the  temperature  at 
a  point  in  a  combusting  system,  the  products,  heat 
release,  and  other  parameters  can  be  estimated  quite' 
well.  For  this  reason,  highly  accurate,  spatially 
precise  measurements  of  temperature  are  a  necessity 
in  understanding  and  controlling  combustion  processes. 


CARS  thermometry  has  been  applied  to  many  diverse 
systems  ranging  from  carefully  controlled  ’aboratory 
flames  to  such  difficult  cases  as  burning  propellent 
strands.  Internal  combustion  engines,  gas  turbine 
combustor  cans,  furnaces,  shock  tubes,  and  plasmas 
are  among  the  several  devices  which  have  been  probed 
by  the  CARS  method.  Some  of  these  applications  will 
be  covered  in  detail,  following  a  discussion  of  the 
theoretical  and  experimental  aspects  of  CARS  measure¬ 
ments.  A  comparison  of  the  CARS  spectra  of  several 
molecular  species,  of  particular  interest  in  com¬ 
bustion,  which  have  been  used  for  thermometry  is  3lso 
presented.  The  effect  of  high  pressure  on  the  CARS 
nitrogen  spectrum  is  discussed  to  emphasize  the 
dramatic  spectral  changes  which  occur  with  pressure 
increase , 

CARS  FUNDAMENTALS-  THEORETICAL  CONSIDERATIONS 

Temperature  measurements  derived  from  spectro¬ 
scopic  techniques  ultimately  depend  upon  the  Boltzmann 
equation,  which,  of  course,  determines  the  energy 
level  population  distribution  with  temperature.  This 
fundamental  dependence  is  also  true  of  CARS  tem¬ 
perature  measurements;  however,  because  of  the  non¬ 
linear  nature  of  the  CARS  effect,  the  derivation  of 
temperature  is  somewhat  more  complicated.  This 
complication  makes  the  fitting  of  computer-modelled 
CARS  spectra  with  experimental  CARS  Bpectra  mandatory. 
For  this  reason  an  analytical  modelling  capability  is 
of  equal  Importance  to  the  experimental  program  in 
any  CARS  Installation. 

CARS  is  a  coherent,  nonlinear  process  which 
involves  the  mixing  of  four  optical  fields  in  a 
material  medium  (liquid,  gas,  or  solid).  In  the 
usual  CARS  experiment,  two  of  the  fields  are  derived 
from  the  same  laser  beam,  called  the  pump  beam,  of 
frequency  The  other  beam,  the  Stokes  beam,  has 

frequency  which  is  less  than  u, ,  The  top  portion 

of  Fig,  1  sketches  one  experimental  means  of  mixing 
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“1  and  “2  t0  Senerate  u>3,  the  CARS  frequency.  The 
middle  section  of  Fig,  1  demonstrates  the  photon 
energy  conservation  requirement  for  CARSj  namely, 
that 

“3  -  -  “2  U> 

Mote  also  that  -  Wj  =  u»  .  When  corresponds  to 
a  Raman-active  vibrational-rotational  transition  in 
the  molecule  of  interest,  the  CARS  effect  is  strong, 
because  of  a  quantum-mechanical  resonance  which  will 
not  be  discussed  here.  Because  of  this  resonance, 
the  CARS  effect  is  usually  much  stronger  than  the 
conventional  spontaneous  Raman  effect.  Equally 
important  is  the  fact  that  the  CARS  radiation  is 
generated  in  a  laser-like  beam  nearly  identical  to 
the  input  beams.  This  means  that,  unlike  Raman 
scattering,  the  entire  CARS  beam  can  be  collected. 
Furthermore,  this  provides  good  discrimination 
against  background,  if  aperturlng  is  employed. 
Obviously,  this  provides  a  tremendous  advantage 
when  making  measurements  in  highly  luminous  com¬ 
bustion  devices.  This  point  will  be  illustrated 
in  more  detail  in  the  section  on  CARS  applications. 

The  bottom  portion  of  Fig,  1  depicts  two 
different  ways  to  generate  an  entire  CARS  spectrum. 
In  the  first  of  these,  the  pump  frequency  u>j  is 
combined  with  a  narrowband  Stokes  frequency,  (02, 
which  is  tuned,  generating  the  CARS  spectrum  step¬ 
wise.  The  second  method  shown  employs  a  broadband 
Stokes  laser,  u>2,  which  in  combination  with  , 
generates  the  entire  CARS  spectrum,  <*>3,  from  each 
laser  pulse.  This  latter  method  offers  the 
potential  for  making  temperature  measurements  from 
a  single,  10~8  second,  laser  pulse,  as  will  be  shown 
later. 


•  Approach 


The  thereotlcal  basis  for  CARS  will  be  briefly 
sketched)  for  a  more  detailed,  yet  readable,  account, 
the  reader  may  refer  to  Ref.  1,  and  to  the  references 
contained  therein.  The  physical  origin  of  the  CARS 
effect  arises  from  the  nonlinear  response  of 
molecules  (both  the  electronic  and  nuclear  portions) 
to  the  incident  electric  fields.  This  response  is 
a  macroscopic  polarization  (a  charge  distortion) 
that  is  proportional  to  the  product  of  the  three 
electric  field  amplitudes.  This  time  dependent 
polarization  then  acts  as  an  oscillating  dipole 
source  to  generate  the  CARS  radiation.  The  material 
parameter  which  relates  the  macroscopic  polarization 
to  the  product  of  the  electric  fields  is  named  the 
third  order  nonlinear  susceptibility,  It  *s 

a  complex  quantity  and  depends  upon  all  molecules 
present  in  the  sample.  The  electronic  portion  of  x^ 
is  nearly  frequency  Independent  (unless  an  electronic 
transition  is  present)  and  is  termed  the  nonresonant 
susceptibility,  xnr«  In  contrast,  the  nuclear  response 
of  ^O)  £s  strongly  frequency  dependent  and  complex, 
and  is  called  the  resonant  susceptibility.  The 
total  susceptibility  Is  written  as 

^  nr 

X1  -  £  (x’  ^  iX  ) ,  +  x  .  (2) 

j 

where  the  bracketed  term,  the  resonant  portion,  is 
summed  over  all  transitions,  j.  In  terms  of 
and  the  incident  laser  intensities,  the  CARS  intensity 
is  given  by  (for  the  case  of  plane  wave  excitation, 
which  is  sufficient  for  our  use  here) 


where  n^  is  the  refractive  index  at  frequency  0)3;  Ij 
and  I2  are  the  laser  intensities  at  frequencies 
and  u>2  respectively.  L  is  an  interaction  length  over 
which  the  input  laser  beams  are  strongly  focused,  and 
phase  matched.  K  Is  a  constant.  Phase  matching  is 
achieved  when  the  vector  equation 


is  satisfied.  The  wave  vector  k  is  defined  by 

1*1 1  *  n^u>i/c,  n^  is  the  refractive  Index  at  frequency 

Equation  (3)  assumes  the  form  shown  when  phase 
matching  is  satisfied.  Phase  matching  is  mainly  a 
problem  of  directing  the  laser  beams  appropriately 
and  will  be  discussed  in  detail  In  the  experimental 
section. 

The  third-order  susceptibility  can  be  derived 
quantum-mechanical ly  and  expressed  in  spectroscopic 
terms  as 


•  Spactrum 

Scanned 


“2  *>i  *3 


Fig.  1.  Coherent  anti-Stokes  Raman  spectroscopy  (CARS). 


»  j  2(u1--2-«J)-ir  *  (5) 

where  xnr  has  been  separated  from  the  reeonent 
susceptibility,  the  summed  term.  The  sum  is  often 
restricted  to  so-cslled  Q-branch  transitions,  with 
&J  •  0  (no  change  in  rotational  quantum  number) 
because  the  Q-branch  transitions  often  dominate  the 
CARS  spectrum.  In  equation  (5) ,  N  is  the  number 


density  of  the  molecular  species  of  interest;  Oj , 
Uj,  and  Tj  ere  the  polarleablllty  matrix  element 
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(related  to  the  Rattan  cross  section),  the  transition 
frequency,  and  the  Isolated  Raman  llnevidth 
respectively,  of  transition  Q(j).  The  term  is 
the  Boltzmann  population  difference  and  is 
responsible  for  the  temperature  dependence  of 
Extracting  the  temperature  from  a  CARS  signal,  which 
depends  upon  the  squared  modulus  of  x^  •  is  not  a 
simple  task  and  the  assistance  of  computed  model 
spectra  is  a  necessity,  as  previously  stated. 

Several  examples  of  computer -mod el ,  predicted  CARS 
spectra  which  demonstrate  considerable  sensitivity 
to  temperature  will  be  shown  in  a  later  section. 

In  simple  terms,  the  temperature  is  derived  from  the 
spectral  distribution,  i,e.#  the  "shape"  of  the  CARS 
signal.  What  factors  determine  the  shape  requires 
closer  examination.  Based  upon  consideration  of 
equation  (2),  there  arc  two  limiting  cases  which 
determine  the  shape  of  a  CARS  spectrum  as  shown 
in  Fig.  2,  For  the  "strong  signal"  case,  it  is 
assumed  that  the  resonant  susceptibility  is  large 
compared  to  xnr*  This  leads  to  the  resonance  type 
of  line  shape  indicated.  An  example  of  the  strong 
case  would  be  that  of  a  majority  species  such  as 
nitrogen  in  an  air-fed  combustion  system.  CARS 
spectra  of  the  strong  line  case  provide  a  basis 
for  accurate  temperature  measurement.  The  other 
extreme,  the  "weak  line"  case,  occurs  when  xnr  is 
large  compared  to  the  resonant  susceptibility.  The 
desired  resonant  CARS  signal  then  appears  as  a 
modulation  or  dispersion-type  spectrum  on  the  ynr 
background.  This  type  of  signal  can  be  extremely 
useful  for  determining  the  concentration  of  a 
species  by  fitting  the  shape  of  the  spectrum,  without 
making  an  absolute  signal  intensity  measurement. 

This  topic  will  be  treated  briefly  in  a  later  section 
of  this  paper,  A  more  complete,  detailed  description 
of  the  modelling  of  CARS  spectra  can  be  found  in 
Refs.  2,  3.  Those  papers  present  the  calculation  of 
X<3>  from  spectroscopic  constants  and  the  Boltzmann 
distribution  of  energy  levels.  Instrumental  factors 
such  as  laser  llnevldths  and  spectrometer  resolution 
are  taken  into  consideration. 


STftONO  LINE  WEAK  LINE 


Fig.  2.  Illustration  of  strong  and  weak  line  CARS  in¬ 
tensity  profiles. 

CARS  FUNDAMENTALS i  EXPERIMENTAL  CONSIDERATIONS 

CARS  is  a  coherent,  nonlinear  optical  process, 
as  stated  previously.  Laser  sources  are  required  in 
order  to  achieve  signal  coherence  (obtain  phase 
’Hatching),  and  nonlinearity  of  response  (requires 
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Fig,  3.  Illustration  of  CARS  phase-matching  approaches. 
Subscripts  denote  beams:  1,  pump;  2,  Stokes;  3,  anti- 
Stokes. 


high  intensities).  Although  CARS  has  been  observed 
with  continuous  lasers  (Ref.  4),  high  peak  power, 
pulsed  lasers  are  almost  universally  employed  for  the 
generation  of  CARS,  particularly  for  gas  phase 
diagnostics.  This  is  the  type  of  system  which  will 
be  described  in  this  section  after  a  short  description 
of  phase  matching. 

In  order  to  obtain  the  maximum  CARS  signal ,  che 
incident  laser  beams  must  be  so  aligned  that  the 
wave  mixing  process  is  properly  phased.  The  CARS 
signal  generated  at  a  certain  point  will  be  in  phase 
with  CARS  generated  at  other  points  in  space  leading 
to  a  constructive  buildup  of  signal  (rather  than 
destructive  interference,  if  out  of  phase).  In 
Fig.  3a  is  shown  the  pictorial  representation  of 
equation  (4),  the  general  condition  for  phase 
matching.  Underneath  the  vector  diagrams  are  shown 
the  experimental  arrangement  of  the  beams.  In  a 
gas,  the  refractive  index  changes  little  with  wave¬ 
length,  and  the  collinear  arrangement  of  Fig.  3b 
satisfies  the  phase-matching  condition.  Although 
very  easy  to  implement,  the  collinear  configuration 
has  very  poor  spatial  resolution,  as  demonstrated  by 
Eckbreth  (Ref.  5).  The  length  of  sample  interrogated 
by  the  CARS  effect  depends  upon  the  focal  length  of 
field  lens  and  can  be  quite  large.  For  example, 
a  10  cm  lens  yields  a  sample  length  of  only  %  0.07  cm, 
but  a  100  cm  lens  (which  may  be  required  in  some 
Instances)  generates  a  sample  length  of  ^7  cm. 
Clearly,  using  a  long  focal  length  lens,  a  collinear 
CARS  measurement  over  a  small  scale  burner  would 
sample  both  hot  and  cold  regions  and,  hence,  would 
be  considerably  misleading.  Recognizing  this, 

Eckbreth  (Ref.  5)  devised  a  crossed-beam  phase- 
matching  scheme  Illustrated  in  Fig.  3c,  called 
BOXCARS  because  of  the  shape  of  the  diagram.  In 
this  case,  CARS  is  generated  only  at  the  point  where 
the  three  beams  cross,  which  can  have  a  dimension 
along  the  lens  axis  of  ^  1  mm,  and  less  in  the 
transverse  direction.  A  further  advantage  of  BOXCARS 
is  that  the  CARS  beam  is  almost  completely  spatially 
(hence  spectrally)  separated  from  the  input  laser 
beams,  and  is  easier  to  isolate  than  the  collinear 
case.  Complete  separation  of  the  CARS  beam  is 
achieved  if  the  BOXCARS  diagram  is  folded  (hence, 
folded  BOXCARS)  about  the  axis  as  shown  in  Fig,  4b, 

The  arrangement  of  the  Incident  and  exiting  beams, 
also  shown,  clearly  demonstrates  how  the  beam 
is  completely  separated  from  the  <4}  and  <*2  beams. 

This  spatial,  and  consequent  spectral,  separation 
is  particularly  advantageous  for  rotational  CARS 
experiments  (Ref .  6) ,  where  the  spectral  separation 


607 


Verdieck  etui. 


from  input  frequencies  is  small.  Rotational  CARS 
spectra  obtained  by  use  of  folded  BOXCARS  are 
discussed  in  a  later  section.  There  are  several 
other  types  of  CARS  phase-matching  schemes  which 
will  not  be  discussed  here.  These  methods  are 
described  and  referenced  in  Ref.  1. 


•)  Plant! 


Tig.  4.  CARS  crossed-beam  phase-matching  approaches 
showing  tiie  phase-matching  diagram  and  actual  geometry 
of  the  optical  beams  for  (a)  planar  and  (b)  folded 
BOXCARS.  Subscripts  denote  beams;  1,  pump,  2,  Stokes, 
),  anti -Stokes. 

A  typical  broadband  CARS  arrangement  for 
combustion  diagnostics  is  schematically  illustrated 
in  Fig.  5  (Refs.  1,  7,  8).  The  neodymium: YAG  laser 
provides  two  2nd  harmonic  beams  at  532  nm  by 
frequency-doubling  the  fundamental  1.06  pm  output 
from  the  laser,  then  frequency-doubling  the  residual 
1.06  u  emerging  from  the  first  doubler.  The  primary 
2nd  harmonic,  i.e.,  after  the  first  doubler,  is 
typically  200  mJ,  while  the  secondary,  i.e.,  from 
the  second  doubler,  is  about  an  order  of  magnitude 
less.  The  secondary  passes  through  a  slab  of  KG3 
Schott  glass  placed  at  the  Brewster  angle  to  absorb 
any  remaining  1.06  pra.  The  secondary  is  then  directed 
by  mirrors  and  focussed  to  pump,  slight!  v  off  axis, 
the  Stokes  dye  cell  oscillator.  The  output  from  the 
oscillator  is  amplified  in  a  second  dye  cell,  pumped 
by  splitting  a  portion  of  the  primary  (typically  33%) 
with  a  beamsplitter  as  shown.  The  Stokes  laser  and 
532  nm  pump  component  each  pass  through  Gian  laser 
polarizers  to  ensure  polarization  purity,  which  is 
important  when  performing  polarization  sensitive 
CARS  experiments.  Each  beam  passes  through  expansion 
or  contraction  telescopes  whose  function  will  be 
detailed  shortly.  The  Stokes  beam  passes  through  a 
rotatable  optical  flat  and  then,  in  folded  BOXCARS 
as  shown,  directly  to  the  focussing  lens.  In  planar 
BOXCARS  (Ref.  5)  the  Stokes  beam  would  pass  through 
a  dichroic  mirror  used  to  reflect  one  of  the  primary 
pump  components  as  depicted  in  Fig.  5,  One  advantage 
of  folded  BOXCARS  is  the  elimination  of  this  dichroic 
element.  Another  advantage  of  folded  over  planar 
BOXCARS  is  the  complete  angular  and  spatial  separation 


Fig.  5.  Schematic  of  CARS  r  imeiit  .irr  in.-.o’u  ni  . 


of  the  CARS  beam  which  eliminates  the  need  to  disperse 
or  filter  the  CARS  beam  from  one  of  the  pump  com¬ 
ponents.  The  primary  pump  beam  is  split  at  a  second 
beamsplitter  (50%)  to  produce  the  two  components  for 
BOXCARS.  These  beams  bass  through  low  order  half 
wave  plates  which  control  the  polarization  orientation 
of  the  pump  components.  The  pump  components  are 
directed  by  mirrors  to  the  focussing/crossing  lens. 

If  the  three  beams  are  aligned  parallel  to  each  other, 
they  will,  by  definition,  cross  at  the  focal  point  of 
the  lens.  However,  due  to  wavefront  sphericity  and 
chromatic  aberatlon,  they  may  not  necessarily  waist 
at  the  crossing  point.  The  adjustable  telescopes 
are  used  to  position  the  beam  waists  at  the  crossing 
point.  At  high  laser  intenstit ies ,  the  beams  are 
readily  visualized  near  the  focal  region  from 
the  room  air  Rayleigh  scattering.  The  primary 
component  waists  will  coincide  only  if  the  distances 
travelled  by  each  pump  component  from  the  second 
beamsplitter  to  the  focussing  lens  are  equal.  The 
telescopes,  depending  on  their  magnification,  also 
permit  the  focal  diameters  to  be  varied.  Rotation 
of  the  optical  flat  in  the  Stokes  beam  permits  its 
displacement  on  the  focussing/crossing  lens, 
permitting  the  phase-matching  angle  6  to  be  varied. 
After  passing  through  the  crossing  point,  the  four 
beams  i.e.,  CARS  at  ,  u>j  ,  wj  t  ,  are  recollimated 
by  a  second  lens  of  generally  the  same  focal  length 
as  the  focussing  lens.  The  unwanted  frequencies  are 
trapped  and  the  spatially  separated  CARS  beam  is  split 
into  two  components,  not  necessarily  equivalent,  by 
a  beamsplitter.  The  component  split  off  at  right 
angles  passes  through  a  Gian  Thomson  polarization 
analyzer  set  normal  to  the  polarization  of  the  resonant 
mode  CARS  signal.  Thus,  only  the  nonresonant  CARS 
signal,  acutally  some  fraction  thereof,  is  transmitted 
and  is  monitored  after  spectral  filtering  by  a  photo¬ 
multiplier.  This  signal  is  used  for  ln-situ 
referencing  which  will  not  be  treated  here.  The 
portion  of  the  CARS  radiation  transmitted  by  the  beam¬ 
splitter  passes  through  a  polarization  analyzer  and  is 
directed  to  a  spectrograph  fitted  with  an  optical 
multichannel  detector  (0MD),  or  a  monochromator 
equipped  with  a  photomultiplier  for  spectral  scanning. 
The  second  polarization  analyzer  may  allow  any  con¬ 
tributions  to  the  CARS  radiation,  i.e.,  isotropic 
and  anisotropic  resonant  modes,  nonresonant  signals, 
to  pass.  Or,  It  may  be  oriented  to  suppress  detection 
ol  one  or  two  of  the  above.  A  small  fraction  of 
the  CARS  beam  is  split  off  at  a  glass  slide  before 
the  spectrograph  and  sent  to  a  photomultiplier  tube 
fitted  with  an  appropriate  narrowband  interference 
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filter.  This  signal  is  averaged  on  a  BOXCAR 
averager  and  used  to  monitor  and  "tweak"  the 
alignment  of  the  CARS  system.  It  is  relatively 
immune  to  angular  and  translational  motions  of  the 
CARS  beam  which  occur  during  peaking  up  the  alignment. 
The  alignment  of  a  CARS  system  should  never  be 
monitored  through  a  narrow  angular  and  spatial 
acceptance,  such  as  a  slit,  where  signal  loss  due 
to  slight  steering  can  be  misinterpreted  for 
alignment  detuning.  The  spatial  resolution  of  a 
given  configuration  is  readily  checked  by  generating 
CARS  from  within  a  translatable,  thin  microscope 
slide  cover  with  the  laser  beams  suitably  attenuated 
to  prevent  optical  damage  to  the  glass.  Microscope 
slide  covers  are  generally  less  than  200  thick  and 
allow  the  pointwise  CARS  signal  contribution  to  be 
ascertained.  Upon  integration  of  the  pointwise 
contributions,  the  spatial  resolution  can  be  found. 

In  general,  it  is  far  better  practice  to  measure 
the  spatial  resolution  than  to  calculate  it. 

One  advantage  of  broadband  CARS  is  the  simplicity 
of  the  Stokes  dye  laser  which  requires  no  elements 
for  tuning  and  spectral  condensation.  A  flowing 
dye  cell,  oriented  at  Brewster's  angle,  resides 
within  a  planar  Fabry-Perot  oscillator  cavity. 

Because  the  dyes  amenable  to  532  am  pumping  typicall> 
exhibit  high  gain,  slighly  off-axis  pumping  works 
quite  well  and  leads  to  very  good  beam  quality,  i.e., 
low  angular  beam  divergence,  in  contrast  to  that 
often  obtained  with  transverse  pumping.  The  dye 
spectrum  is  centered  at  the  desired  wavelength  by 
selecting  the  dye  appropriately  and  by  adjusting  its 
concentration.  Solvent  tuning  can  also  be  employed. 

With  the  planar  Fabry-Perot  dye  oscillator 
arrangement,  bandwidths  vary  from  100  to  200  cm”1, 
depending  on  the  pump  energy  and  dye  employed.  Binary 
dye  mixtures  are  often  used  to  improve  dye  conversion 
efficiency  in  the  desired  wavelength  region.  For 
single  pulse  CARS  diagnostics  or  for  spectral 
scanning  of  laminar  flames,  it  is  important  that 
the  dye  spectrum  be  smooth  and  reproducible  from 
pulse  to  pulse.  The  dye  cavity  is  purposely 
designed  to  have  a  high  Fresnel  number  and  to 
accommodate  as  many  modes  as  possible  to  "fill  in" 
the  spectrum.  The  pulse  to  pulse  spectral  stability 
of  the  dye  laser  has  been  examined  on  the  0MD  and 
found  to  be  fairly  good.  Single  pulse  dye  spectra 
display  an  irregular  fine  structure  with  an  amplitude 
variation  of  +  5%,  presumably  due  to  spatial  "hole 
burning"  in  the  dye  laser.  This  fine  structure 
averages  out  in  time  to  produce  spec t rally-smooth 
profiles.  High  quality,  averaged  CARS  spectra  are 
readily  achieved,  but  single  pulse  CARS  spectra  may 
exhibit  some  spurious  structure,  particularly  at 
moderate  spectral  resolution,  <  1  cm"*.  At  lower 
resolutions,  >.  1  cm”1,  this  fine  structure  is 
■spectrally  smeared  and  is  not  a  problem.  In  broad¬ 
band  CARS,  either  the  linevldth  of  the  pump  laser  or 
the  resolution  of  the  monochromator/spectrograph 
determines  the  ultimate  resolution  of  the  spectrum 
1-m  monochromators  typically  have  a  limiting 
resolution  of  about  0.5  cm"1  in  the  visible.  2xNd 
lasers  with  intracavity  etalons  have  llnewidths  in 
the  0.1  to  0.4  cm”1  range,  and  0.8  cm"1  without. 

CARS  spectra  so  obtained  thus  have  a  limiting 
resolution  between  0.5  and  1  cm”1.  This  moderate 
resolution  is  generally  more  than  sufficient  for 
diagnostics.  Furthermore,  it  is  Important  to  note 
that  the  large  pump  laser  linevidth  is  not  detri¬ 
mental  in  regard  to  the  strength  of  the  CARS 
radiation.  For  broadband  CARS  generation,  it  is 


easy  to  show  that  the  spectrally  Integrated  intensity 
of  a  CARS  transition  is  independent  of  the  pump  laser 
linewidth. 

Narrowband  CARS  systems  for  combustion 
diagnostics  (Ref.  9)  are  conceptually  similar  to  that 
described  above  with  the  obvious  exception  of  the 
tunable  Stokes  source.  For  narrowband  Stokes  work, 
a  variety  of  dye  laser  configurations  exist,  are 
commercially  available  or  readily  assembled.  The 
oscillator  sections  are  generally  transversely-pumped 
to  obtain  gains  high  enough  to  overcome  the  insertion 
losses  engendered  by  the  spectral  condensation  scheme 
employed.  The  Hansch  design  (Ref.  10)  of  a  circular 
telescope,  large  two  dimensional  grating  has  generally 
given  way  to  one  dimensional  expansion  schemes 
employing  a  grazing  incidence  grating  (Ref.  11)  or 
multiple-prism  beam  expanders  (Ref.  12).  The  one- 
dimensional  expanders  require  a  grating  large  in  one 
dimension  only  and  greatly  reduce  cost. 

The  output  from  the  nochromater  can  be  detected 
by  a  photomultiplier  tube  or  with  an  optical  multi¬ 
channel  analyzer  OMA  (Ref.  13).  Photomultiplier 
signals  are  usually  averaged  with  a  boxcar  integrator. 
Detection  with  an  OMA  is  a  better  method,  because  the 
entire  spectrum  from  each  CARS  pulse,  is  captured  and 
many  spectra  may  be  summed  over  a  desired  number  of 
pulses.  Given  sufficient  signal  strength,  a  single 
CARS  pulse  (typically  10  nanoseconds)  can  provide 
enough  information  to  extract  the  temperature,  as 
stated  previously.  Moreover,  if  the  OMD  signal  is 
processed  with  an  on-line  laboratory  minicomputer 
which  can  store  model  predicted  CARS  spectra,  then 
temperature  measurements  can  be  made  in  a  matter  if 
seconds  (10  seconds  corresponds  to  100  laser  pulses) . 
This  capability  is  presently  available  at  UTRC. 

CARS  THERMOMETRY 

Temperature  can  be  derived  from  a  CARS  spectrum 
because  the  vibrational  -  rotational  quantum  state 
populations  depend  upon  temperature.  Statistical 
mechanics  tells  us  that  the  probability,  P  <v,J)  that 
a  state  with  vibrational  and  rotational  quantum 
numbers,  v  and  J,  is  populated  is  given  by  (to 
first  order  of  approximation) 

-E  /RT  -B  J(J+1)/RT 
gj(2  J+l)  e  V  e  V 

P(V’J)  '  ^  <6> 

where  gj  is  the  nuclear  spin  degeneracy;  and  QR(v) 
are  respectively  the  vibrational  and  rotational 
partition  functions  for  level  v;  Ey  li  vibrational 
energy  and  By  is  the  rotational  constant  in  state  v; 
k  is  the  Boltzmann  constant  and  T  the  equilibrium 
temperature.  For  most  molecules  of  interest  in  com¬ 
bustion  the  vibrational  CARS  spectrum  is  dominated 
by  Q-b ranch  transitions,  whose  frequencies  are  given 
by  (for  a  diatomic  molecule) 

WVJ  *  “e  *  2  "eXe  (V+1)  -  ae  J(J+1)  (7) 

Thus,  the  two  equetions  (6)  and  (7),  through  the  basic 
equations  for  the  CARS  effect  discussed  In  the 
theoretical  section,  deteralne  the  shape  of  the 
resonant  CARS  spectrun.  Experience  In  this  laboratory 
has  demonstrated  that  temperature  determination  from 
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experimental  CARS  spectra  requires  fitting  to 
computer-generated  spectra  (Refs.  2,  3,  14).  The 
possible  exception  might  be  for  hydrogen  CARS  spectra, 
where,  because  of  the  large  vibration-rotation 
interaction,  ae  -  2.933  cm"1,  the  Q-branch  transitions 
are  well  separated  and  temperature  can  be  found  from 
relative  peak  heights.  However,  hydrogen  is  present 
in  abundance  in  combustion  systems  only  for  special 
cases.  A  more  useful  CARS  probe,  present  in  large 
quantity  in  all  air-supported  combustion,  is  the 
nitrogen  molecule. 

CARS  NITROGEN  THERMOMETRY 

In  addition  to  the  nearly  ubiquitious  presence 
of  nitrogen,  the  spectroscopic  parameters,  including 
the  Raman  linewidth,  needed  for  accurate  CARS  spectral 
modelling,  are  quite  well  known.  That  the  N2 
vibrational  CARS  spectrum  provides  an  excellent  basis 
for  thermometry  can  be  appreciated  through  examination 
of  the  computer  calculated  temperature  dependence  of 
the  N2  spectrum  shown  in  Fig.  6.  These  calculations, 
carried  out  by  Hall  (Ref.  3),  correspond  to  multiplex 
(broadband  Stokes)  CARS  employing  a  0.8  cm-1  band¬ 
width  u>2  ,  a  130  cm-1  bandwidth  Stokes  and  a  1  cm 
instrumental  slit  function.  Clearly,  these  spectra 
exhibit  pronounced  temperature  sensitivity,  partic¬ 
ularly  at  higher  temperatures,  and  show  much  spectral 
detail  even  though  the  instrumental  resolution  is  not 
particularly  high.  These  spectra  all  correspond  to 
the  "strong  line"  case.  It  is  assumed  that  nitrogen 
is  present  at  70%,  hence,  there  is  little  contribution 
from  the  nonresonant  susceptibility  background 
although  its  inclusion  in  the  calculations  is 
important  for  accurate  thermometry.  Some  notable 
features  should  be  pointed  out.  At  room  temperature 
a  single  narrow  peak,  without  structure  at  this 
resolution,  is  observed.  This  peak  broadens  with 
increasing  temperature  and  shows  resolved  Q-branch 
transitions  corresponding  to  Q(J  =  20)  to  Q(J  =  40) 
at  about  1000  K.  The  peaks  observed  correspond  to 
even-J  Q-branch  transitions;  the  odd-J  transitions, 
are  reduced  in  intensity  by  one-fourth  (nuclear 
statistical  factors  are  squared  in  CARS,  because 
X  is  squared) •  A  second  vibrational  band,  a  "hot 
band,"  v*»l  to  v=2  ,  becomes  noticable  above  2000  K, 
but  is  not  shown  in  Fig.  6.  For  the  moderate 
resolution  conditions  stated  (and  readily  achieved) 
for  Fig,  6,  vibrational  CARS  thermometry  would  be 
most  accurate  for  measurements  above  ^  1000  K, 

Given  better  instrumental  resolution,  and  lasers 
of  smaller  bandwidth  (both  requiring  greater  cost 
and  effort),  more  spectral  detail  would  be  observ¬ 
able  at  lower  temperatures,  which  would  increase 
the  accuracy  of  CARS  measurement  at  lower  temper¬ 
atures.  In  contrast,  under  conditions  of  poorer 
spectral  resolution,  when  the  fine  structure  is 
lost,  temperature  still  can  be  found  from  the 
width  of  the  fundamental  and/or  the  relative 
strength  of  the  hot  band.  From  the  preceding 
discussion,  it  is  seen  that  vibrational  CARS 
thermometry  can  be  applied  over  a  wide  range  of 
temperature  with  a  varying  degree  of  accuracy. 

As  an  example  of  the  type  of  accuracy  that  can 
be  achieved  with  CARS  N2  thermometry,  experiments 
were  performed  at  UTRC  in  premixed  flat  flame 
burners  (Refs.  14,15).  Comparison  of  the  CARS  determined 
temperatures  was  made  with  radiation-corrected  thermo¬ 
couples  (in  turn  calibrated  by  sodium  line  reversal) . 
Figure  7  shows  a  scanned,  broadband  Nj  CARS  spectrum 
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Fig.  6.  Temperature  variation  of  the  CARS  spectrum  of 
N2  for  a  pump  linewidth  of  0.8  cm"1,  a  Stokes  bandwidth 
of  150  cm"1  and  a  1  cm"1  slit  width. 


taken  in  a  clean  (non-sooting)  methane/air  flame  at 
a  point  where  the  temperature,  measured  by  the 
thermocouple,  was  determined  to  be  2110  K.  The 
dotted  curve  is  the  computer  model,  least  mean 
squares  fit  to  the  experimental  spectrum;  the  best 
fit  temperature  is  2104  K  +  9  K.  This  corresponds 
to  an  accuracy  of  better  than  1.0%.  In  order  to 
achieve  this  type  of  accuracy  for  the  theoretical 
fit,  very  accurate  spectral  constants  for  the  N2 
vibrational-rotational  energy  levels  are  required. 

In  particular,  the  Raman  linewidths,  which  enter 
into  the  expression  for  the  resonant  susceptibility 
equation  (5) ,  must  be  known  as  a  function  of  tem¬ 
perature  and  J-state.  Fortunately,  accurate 
measurements  of  the  Raman  linewidths  of  N2  have 
been  obtained  by  Rahn  and  co-workers  (Ref.  16)  at 
300  K  and  1700  K  using  high  resolution  stimulated 
Raman  gain  spectroscopy.  Equally  important,  a 
theoretical  model  for  Raman  linewidths,  developed  by 
Bonamy  and  colleagues  (Ref.  17),  is  in  good  agreement 
with  the  experimental  values.  This  model,  anchored 
by  the  experimental  points,  permits  extrapolation 
to  higher  temperatures  for  estimates  of  the  Raman 
linewidthB,  with  an  error  probably  less  than  20k, 
Moreover,  the  T  and  J  dependence  of  the  linewidth 
is  weak  at  higher  temperatures.  The  sensitivity  of 
the  CARS  determined  temperature  to  the  Raman  line- 
width  has  been  examined  by  Eckbreth  and  Hall  (Ref.  1  ). 
At  2100  K,  it  is  estimated  that  the  temperature  error 
is  less  than  25  K  (hence,  ^  1%)  for  a  20%  error  in  the 
Raman  linewidths.  Thus,  the  sensitivity  to  linewidth 
is  not  large,  and  extrapolated  values  should  be  quite 
adequate. 
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Fig.  9.  Radial  temperature  profiles  determined  by  CARS 
in  a  laminar  nropane  diffusion  flame. 


Fig.  7.  Comparison  of  experimental  (solid  line)  and 
theoretical  least  squares  fit  (dashed)  N-,  CARS  spectra 
at  a  measured  temperature  of  2110°K. 


To  illustrate  the  application  of  CARS  thermometry 
under  adverse  conditions,  measurements  in  a  highly 
sooting  flame  will  be  described  (Ref.  3).  BOXCARS, 
with  a  spatial  resolution  of  about  0.3  by  1.0  mm, 
was  employed  in  a  laminar  propane/air  diffusion  flame. 
The  axial  and  radial  variation  of  temperature  in  this 
flame  are  shown  in  Figs.  8,  9,  respectively.  The 
quality  of  the  sooting  flame  CARS  spectra  is 
equivalent  to  that  obtained  from  clean  flames.  The 
significance  of  CARS  temperature  measurements  in 
the  highly  luminous,  sooting  flame  Is  that  the  CARS 
method  succeeds  where  conventional  Raman  methods 
would  fail  because  laser-induced  particulate 
incandescence  (Ref.  18)  would  swamp  the  weak  Raman 
signal.  In  contrast,  the  CARS  signal,  which  is 
completely  collected,  is  sufficiently  intense  to 
determine  temperature  from  a  single,  10“^  second 
pulse.  An  example  of  single  shot  CARS  spectra  will 
be  presented  under  the  section  dealing  with 
applications. 


Fig.  8.  Axial  temperature  profile  determined  in  a 
laminar  propane  diffusion  flame. 


PRESSURE  EFFECTS  IN  THERMOMETRY 

To  achieve  maximum  usefulness  as  a  diagnostic 
tool,  CARS  must  be  extended  to  pressures  much  higher 
than  one  atmosphere.  Combustion  at  high  pressure  is 
important  in  such  devices  as  gas  turbines,  internal 
combustion  engines  (especially  diesels),  and  rockets. 
The  pressure  encountered  in  such  systems  ranges  from 
a  few  tens  of  atmospheres  to  well  over  100  atmospheres. 
If  the  theoretical  CARS  analysis  used  to  explain  one 
atmosphere  CARS  flame  spectra  were  employed  at  high 
pressure,  it  would  predict  that  the  overall  band¬ 
width  of  the  specturm  would  increase  linearly  with 
pressure.  The  unmodified  theory  would  also  predict 
that  the  1  -*  2  "hot  band"  would  merge  with  the 
fundamental  band,  and  therefore  would  be  lost  as  a 
sensitive  indicator  of  temperature. 

The  effect  of  pressure  on  CARS  spectra  has 
been  investigated  experimentally  and  theoretically 
over  the  pressure  range  of  one  to  100  atmospheres 
(Ref.  19).  A  similar  study  has  been  performed  by 
Roland  and  Steele  (Ref.  20)  for  pressure  up  to  30 
atmospheres.  The  room  temperature  studies  of  Hall, 
Eckbreth,  and  Verdieck  show  experimentally  that, 
after  an  initial  narrowing  of  the  Q-branch  band 
over  the  first  five  to  ten  atmospheres,  the  width 
of  the  band  remains  constant  with  Increasing  pressure 
to  ICO  atmospheres.  Figure  10  shows  collinear  CARS 
Q-branch  spectra  of  nitrogen  at  one  am*  100 
atmospheres  scanned  with  a  monochromator  spectral 
width  of  'v  0.4  cra-^.  The  dotted  line,  the  theoretical 
fit,  demonstrates  the  excellent  agreement  of  the 
computer  model  modified  to  include  "collision-induced 
band  narrowing".  In  simple  terms,  collisional 
narrowing  of  the  band  occurs  when  the  collision 
frequency  (a  function  of  pressure)  becomes  comparable 
to  the  spacing  between  individual  Q-branch  transitions. 
When  this  happens,  the  molecule  moves  rapidly  through 
many  J  states  and  appears  to  spend,  on  the  average, 
most  of  the  time  in  a  most  probable  J-state,  defined 
by  the  Boltzmann  distribution.  The  phenomenon  has 
also  been  observed  in  NMR  (Ref.  21),  and  Raman 
spectroscopy  (Ref.  22).  For  greater  detail  on  the 
quantum  mechanical  treatment  of  collisional  narrowing 
in  CARS  spectra,  the  Interested  reader  may  consult 
(Refs.  1  and  19).  Figure  11  Illustrates  more 
dramatically  the  Importance  of  collisional  narrowing 
by  a  comparison  of  a  collision-narrowed  line  with  a 
line  calculated  not  including  narrowing.  Note  that 
for  the  latter  case  a  very  broad  line,  a*  10  cm-*, 
results. 


Verdieck  et  sf 


601 


Normalized  intensity 


T-300K 


21122  21112  21122  21112 
Frsqusncy  -  cm- 1 


P  ■  7.0  aim 


Experiment 


Thaory 


Fig.  10.  CARS  signatures  of  300 ®K  N2  at  one  and  one 
hundred  atmospheres  pressure.  The  solid  line  is  the 
experimentally-scanned  spectrum  in  each  case,  the 
dashed  line  the  theoretical  prediction. 


Fig.  LI.  Comparison  of  calculated  100  atm  signature 
with  amt  without  collision.il  narrowing  -  with  narrow¬ 
ing,,  -  without  narrowing. 


The  effect  of  high  pressure  on  high  temperature 
CARS  spectra  has  been  initiated  at  UTRC;  some  results 
are  shown  in  Fig.  12  (Ref.  2  3).  {Experimental  results 
are  shown  at  •<  1600  K  for  pressures  of  •  8  and  28 
atmospheres.  The  model  predicted  spectra  are  presented 
for  comparison.  It  is  noted  that  although  the 
rotational  structure  collapses  over  this  pressure 
change,  the  hot  band  persists;  hence,  the  potential 
for  temperature  measurements  at  high  pressure  still 
exists.  Obviously  the  area  of  high  pressure,  high 
temperature  CARS  spectra  must  be  investigated  more 
carefully  and  thoroughly  in  order  to  assess  CARS 
diagnostics  for  the  several  Important  high  pressure 
applications  listed  earlier. 

Pl'RK  ROTATIONAL  CARS  THKRMOMF.TRY 

Pure  rotational  CARS  may  offer  advantages 
over  vibrational-rotational  CARS  at  high  pressures, 
and  also  at  low  temperatures  (under  1000  K,  at  low 
pressures).  The  reason  for  this  is  that  the  pure 
rotational  lines  have  better  spectral  separation,  4B; 
where  B  is  the  rotational  constant.  In  N2  ,  this 
amounts  to  8  cm* ' .  Pure  rotational  CARS  presents  the 


P«28.2  atm 


Fig.  12.  High  pressure  CARS  spectra  for  S,  it  lnjncK, 

same  classic  problem  that  pure  rotational  :.aman 
spectroscopy  does;  the  frequency  shifts  from  :  lie 
incident  radiation  are  small,  and  the  CARS  radiation 
is  difficult  to  separate  spectral l v  from  the  input 
beams.  One  method,  employed  in  conventional  Ranan 
spectroscopy  is  to  employ  a  double  tor  triple) 
monochromator.  In  this  laboratory  the  net  «.  ssa rv 
spectral  separation  of  the  CARS  beam  is  achieved  In 
the  complete  spatial  separation  of  tlu  beams,  without 
a  dispersing  optic,  through  use  oj  the  three- 
dimensional  phase-matching  technique  called  "folded 
BOXCARS"  (Ref.  6),  J 1  scussed  earlier  in  tin  <  :-.peri- 
nental  section.  Referring  back  to  lie.  ,  it  an  be 
seen  that  the  CARS  bean,  .5,  emerge^  well  separated 
angularly  from  the  input  beams.  In  this  manner,  tin- 
pure  rotational  CARS  spectrur*.  of  room  air  (one 
atmosphere,  300  °K)  was  recorded  bv  spectrally  ;c  anil  in 
the  isolated  beam  (Fig.  I  n.  This  experimental 
spectrum  is  in  good  agreement  with  the  theoretical 
model  (Kef.  ft). 


Foldad  BOXCARS.  Aligned  Polarizations 


Wavalangth  (nm) 


Fig.  1J,  Pure  rotation.il  CARS  spfetra  uJ  air. 
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CARS  THERMOMETRY  WITH  OTHER  MOLECULES 

At  UTRC,  several  other  molecules  of  interest  to 
combustion  science  have  been  examined  for  CARS 
thermometry  and  concentration  measurements.  CARS 
spectra  have  been  obtained  for  hydrogen,  oxygen, 
carbon  monoxide,  carbon  dioxide,  water,  and  methane. 
Space  does  not  permit  a  discussion  of  all  of  these; 
hydrogen,  water,  and  carbon  dioxide  CARS  spectra  are 
selected  for  discussion  with  application  to  temperature 
measurement . 

Computed  CARS  spectra  of  hydrogen  at  four 
different  temperatures  are  illustrated  in  Fig.  14. 

The  vibrational-rotational  interaction,  a,  is  much 
larger  in  hydrogen  (2.993)  than  in  nitrogen  (0.0187); 
therefore,  the  spectrum  is  simple  (few  lines),  and 
quite  spread  out  (over  nearly  200  cm”^) .  The 
alternation  in  intensity  due  to  nuclear  spin 
statistics  is  quite  apparent.  Recall  that  ortho/ 
para  forms  of  hydrogen  bear  statistical  weights  of  3:1. 
However,  in  CARS  spectroscopy,  population  differences 
are  squared  (because  x  Is  squared),  therefore,  the 
intensity  contrast  between  adjacent  lines  is  quite 
large,  and  hence  useful  for  temperature  measurements. 
For  comparison,  experimental  CARS  spectra  for  H0 
recorded  on  an  OMA  and  taken  in  a  H2/air  diffusion 
flame,  are  presented  in  Fig.  15.  The  temperatures 
listed  were  obtained  from  comparison  of  the  peak 
height  ratios  of  Q(l),  Q(3),  and  Q(5)  transitions. 

As  a  practical  application,  temperature  profiles  of 
a  slot-shaped  H^/air  diffusion  flame  were  determined 
from  CARS  spectra,  and  also  from  O2  CARS  spectra 
(02  CARS  spectral  temperatures  are  obtained  in  manner 
analogous  to  that  for  N«)*  an< i  are  displayed  in 
Fig.  16.  The  CARS  derived  temperatures  are  in  fairly 
good  agreement  with  a  radiation-corrected  thermo¬ 
couple.  A  detailed  discussion  of  all  aspects  of  the 
hydrogen  CARS  studies  is  found  in  Ref.  24. 


Fig.  14,  Theoretical  CARS  spectra  of  H2  over  the  tem¬ 
perature  range  1000°K  <  T  1  2000°K. 


Because  water  is  a  major  product  of  air  supported 
combustion  of  hydrogen-containing  fuels,  it  can 
serve  as  a  thermometr ic  probe  molecule,  if  present  in 
sufficient  abundance.  Moreover,  the  measurement  of 


Fig.  15.  CARS  spectra  of  H2  in  H2~air  diffusion  flame 
at  temperatures  determined  from  relative  strengths  of 
indicated  Q-branch  transitions.  Frequency  scale 
corresponds  to  0.60  cm~^  per  dot. 


A h2cars 

A 02 CARS 
o  Thermocouple 


Fig.  16.  Temperature  measurements  in  a  flat  H2-air 
diffusion  flame.  Symbols:  circles,  radiation-corrected 
thermocouple;  open  triangles,  H2  CARS;  and  solid  tri¬ 
angles,  02  CARS.  Dotted  curve  is  locus  of  maximum 
temperature. 

water  concentration  as  a  function  of  distance  in  a 
combustion  system  could  yield  chemical  kinetic 
information  and  a  better  understanding  of  the 
combustion  process. 

The  water  molecule  is  an  asymmetric  top  with 
three  vibrational  modes,  of  which  the  symmetric 
stretch  is  strongly  Raman  allowed,  with  a  Raman  shift 
of  3657  cm-*.  Because  it  is  an  asymmetric  top,  the 
rotational  structure  is  complex  and  the  rotational 
energy  levels  cannot  be  expressed  in  closed  form. 
However,  the  structure  of  the  Q-branch  of  the  mode 
is  tractable,  and  Hall,  et  al.  (Ref.  25)  have  modelled 
the  CARS  spectrum.  Figure  17  displays  the  comparison 
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Fig.  17.  Comparison  of  r.ilcul at  ul  CARS  ‘pec t rum  of 
iure  H90  with  spectrum  measured  n  a  heated  cell  at 
/73°K. 

between  computed  and  experimental  CARS  spectra  of 
water  in  a  static,  heated  cell  at  771  K.  Because 
Raman  linewidth  data  were  not  available  for  this 
calculation,  a  value  of  0.*)  cm  independent  of  .1 
was  selected  on  the  basis  of  tho  good  agreement 
achieved  between  the  calculated  and  experimental 
aspects.  A  similar  comparison  of  theory  and 
experiment  is  shown  in  Fig.  1-S  for  in  a  1700  K 

premixed  methane/air  flame.  Because  of  the  higher 


temperature,  the  assumed  Raman  linewidth  was 
estimated  to  be  0.2  cm"^.  In  these  calculations, 
nearly  one  thousand  vibrational-rotational  Q-branch 
transitions  were  included  in  the  computed  spectrum. 
Because  the  shape  of  the  water  spectrum  is  quite 
sensitive  to  temperature,  the  water  molecule  may  offer 
the  potential  for  thermometry  over  a  wide  range  of 
temperature . 

The  other  major  product  of  hydrocarbon  combustion 
is  carbon  dioxide.  CO2 ,  a  linear  triatomic  molecule, 
has  three  normal  vibrational  modes,  a  symmetric  stretch 
(v.),  an  antisymmetric  stretch  (vj),  and  doubly- 
degenerate  bend,  (v'2).  Moreover,  an  accidental 
degeneracy  of  and  2v2,  the  prototvpe  case  of  Fermi 
resonance,  mixes  these  two  states  and  makes  both 
transitions  from  the  ground  state  Raman  active.  Not 
only  these  two  states  are  in  Fermi  resonance,  but 
several  other  combinations,  occurring  up  the  vibrational 
;tate  ladder,  are  also  in  Fermi  resonance  which  makes 
the  pr^olem  of  calculating  C09  vibrational  energy 
levels  difficult.  In  spite  of  this  complexity, 
theoretical  CARS  spectra  for  Cf>2  have  been  calculated 
by  Hall  (Ref.  26). 

BOXCARS  spectra  of  < 09  have  been  obtained  by 
Eckbreth  (Ref.  1)  in  several  propane-air  and  CO-air 
flames,  as  shown  in  Fig.  19  where  a  measured  thermo¬ 
couple  temperature  is  given  for  each  spectrum.  The 
experimental  signature  displays  a  fundamental  band  at 
a  Raman  shift  of  about  1388  cnfA,  assigned  to  the 
transition  between  the  ground  vibrational  state  and 
the  first  symmetric  stretch  state.  A  number  of  hot 
bands,  originating  in  vibrational ly  excited  initial 
states,  appear  at  larger  shifts.  The  relative  strength 
of  the  hot  bands  is  seen  to  be  moderately  sensitive  to 
temperature,  making,  the  molecule  potentially  attractive 
for  the rmome t ry . 

Figure  20  compares  the  experimental  CC^  CARS 
spectrum  at  1520  K  with  a  theoretical  spectrum  based 
on  the  isolated  line  approximation  (Ref.  26).  As  can 
be  seen,  the  agreement  is  fairly  good,  except  for  the 
iiighest  hot  bands.  In  the  theoretical  calculation, 
vibrational  energy  levels,  rotational  constants,  and 
polairzabil ity  matrix  elements  were  computed  following 
the  treatment  of  Courtoy  (Ref.  27).  This  basically 
involves  diagonalizing  the  system  Hamiltonian  for  each 
set  of  near  resonance  states.  It  was  found,  however. 


Frequency 


Fig.  18.  Comparison  of  experimental  and  theoretical 
CARS  spectra  ol  water  vapor  in  .1  flame  at  1700<>K.  A 
best-fit  pressure-broadened  linewidth  Y  *  0.19  cm  * 
was  inferred  for  all  transitions  in  the  calculation. 
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•)  1520  K 


Channel  number  (0.34  cm  ~  1/chan) 

Fig.  19.  Experimental  carbon  dioxide  spectra  at 
selected  temperatures. 

cna'  .etter  agreement  with  experiment  resulted  from 
the  use  cf  tabulated  values  for  the  rotational 
constants  (Ref.  28).  Similar  agreement  was  obtained 
with  the  other  experimental  spectra  of  Fig.  19. 
Further  investigations,  mainly  theoretical,  will  be 
needed  before  CO2  CARS  can  be  regarded  as  a  reliable 
diagnostic  tool.  For  example,  the  calcul¬ 
ations  display  some  sensitivity  to  assumed  Raman 
llnewidth,  and  more  information  concerning  the 
magnitude  of  these  quantities  is  needed.  Also,  the 
high  degree  of  line  overlap  within  each  band  may 
make  the  isolated  line  approximation  questionable 
even  at  one  atmosphere  pressure. 

CARS  CONCENTRATION  MEASUREMENTS 

although  this  review  is  intended  to  deal 
prlrarly  with  CARS  temperature  measurements,  species 
concentration  measurements  from  CARS  will  be  included 
for  completeness.  CARS  measurements  of  species 
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FiR.  21.  Computed  CARS  spectral  behavior  of  CO  at 
various  concentrations  at  1800°K. 

concentration  in  combustion  are  made  in  either  of 
two  quite  different  ways;  from  the  absolute  intensity 
of  the  CARS  signal,  or  in  certain  concentration  ranges, 
from  the  shape  of  the  spectrum.  The  latter  case, 
which  will  be  presented  here,  results  when  the 
resonant  and  nonresonant  contributions  are  comparable 
and  the  weak  signal  limit  applies,  shown  in  Fig.  2. 

This  behavior  can  be  better  understood  by  reference 
to  Fig.  21,  which  displays  a  set  of  model  predicted 
spectra  of  carbon  monoxide  at  1800  K.  The  con¬ 
centration  values  range  from  0.5  percent  to  100 
percent;  however,  the  useful  range  extends  to  ^  30 
percent  at  the  high  end.  Beyond  this  value  some 
sort  of  reference  is  required  to  scale  the  intensity. 
Below  ^0.5  percent  level,  the  resonant  signal 
disappers  into  the  baseline. 

Spectral  shape  fitting  for  concentration 
measurements  has  been  demonstrated  experimentally 
in  CO,  as  illustrated  by  the  two  spectra  in  Fig.  22. 

The  upper  portion  is  the  experimental  CARS  spectrum 
(dotted  line)  of  a  2.1  percent  CO  (in  argon) 
calibrated  gas  mixture  at  300  K.  The  thin  solid 
line  is  the  computer  generated  CARS  spectrum  cal¬ 
culated  for  the  same  concentration,  and  is  in 
excellent  agreement.  The  bottom  part  of  Fig.  22 
displays  a  similar  comparison  for  CO  in  a  CH^A^/Ar 
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Fig.  20,  Comparison  of  experimental  and  calculated 
CARS  carbon  dioxide  spectra. 
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flame  doped  with  CO.  The  flame  temperature  was 
2100  K.  The  theoretical  spectrum  was  calculated  for 
3.6  percent  CO,  which  was  determined  with  a  quartz 
microprobe  and  NDIR  CO  detector.  Again  the  agreement 
is  good.  For  more  detailed  information  on  CO  con¬ 
centration  measurements,  including  a  discussion  of 
other  CARS  methods,  the  reader  may  consult  Ref.  29. 
Concentration  measurements  by  the  CARS  spectral 
shape  fitting  technique  have  been  applied  to 
oxygen  (Ref.  30)  and  nitrogen  (Ref.  J).  Spectral 
fitting  could  be  used  for  measuring  water  vapor 
concentrations  as  well. 

Before  closing  the  discussion  of  concentration 
measurement  by  spectral  shape,  it  must  be  emphasized 
that  both  the  temperature,  and  the  nonresonant 
susceptibility,  must  be  known  in  order  to  determine 
the  concentration  of  the  desired  species.  Because 
the  nonresonant  susceptibility  changes  little  (10- 
15%)  through  the  air  fed  combustion  reaction  going 
from  reactants  to  products,  intelligent  estimates 
can  be  made  for  xnr ,  if  the  location,  and  temperature 
at  that  point,  are  known.  From  the  experiments  shown 
above,  it  is  clear  that  concentration  measurement 
by  spectral  shape  fitting  is  a  viable  technique, 
easily  applied,  where  appropriate. 

PRACTICAL  APPLICATIONS  OF  CARS  THERMOMETRY 

In  this  sect  ion  some  illustrations  of  CARS 
measurements,  mainly  temperature  measurements,  in 
industrial  scale  combustion  devices  are  listed.  A 
detailed  description  of  CARS  measurements  performed 
at  UTRC  on  a  gas  turbine  test  combustor  concludes 
this  final  section. 

Apparently  the  first  application  of  CARS 
measurements  on  a  practical,  large-scale  combustion 
device  was  performed  by  Taran  and  co-workers  at  ONERa 


T  =  300  °K  2.1%  CO 


Frequency 


Fig,  22,  Comparison  of  experimental  and  computed 
(solid  line)  CARS  spectra  at  300°K  (2,1%  CO)  and 
2100°K  (3.6%  CO). 


in  1978  (Ref.3i).  Temperature  and  concentration 
measurements  of  N2  and  CO2  were  made  in  a  simulated 
turbomachine  combustor  burning  kerosene.  More 
recently  this  group  has  made  measurements  in  a 
commercial  gas  turbine  combustion  can. 

Similar  measurements  at  Wright-Patterson  AFB  have 
been  carried  out  in  a  large  scale  diffusion  flame 
which  was  stabilized  with  a  bluff  body  (Ref.  32). 

CARS  measurements  in  an  internal  combustion  engine 
have  been  made  by  Stenhouse  et  al,  (Ref,  33)  The  CARS 
input  laser  pulses  were  synchronized  with  the  engine 
cycle  to  generate  the  CARS  spectrum  in  time  steps,  and 
could  be  positioned  at  known  points  on  the  engine  cycle. 
More  recently,  workers  at  Ford  Motor  Company 
Research  have  measured  temperature  and  carbon 
monoxide  concentration  in  a  research  scale,  single 
cylinder  engine  (Ref.  34).  A  special  type  of  non- 
collinear  phae  matching  was  employed  which  achieved 
a  spatial  resolution  of  'v-  2  mm  along  the  beam  and 
100  um  transverse  to  the  beam.  Single  pulse,  10"® 
second,  CARS  spectra  were  obtained,  from  which 
temperature  was  determined.  CARS  measurements  also 
have  been  made  in  a  commercial  diesel  engine  (Ref.  35) 
at  Komatsu  Corporation.  Further  measurements  in 
internal  combustion  engines  using  CARS  are  scheduled 
at  many  laboratories  throughout  the  world. 

In  experiments  at  UTRC,  Eckbreth  (Ref.  7) 
made  CARS  measurements  in  a  0.5  m  diameter  combustion 
tunnel  fitted  with  a  variety  of  burners,  one  of  which 
was  a  gas  turbine  (JT-12)  combustor.  CARS  temperature 
measurements  (using  BOXCARS)  from  nitrogen  were  made 
in  the  primary  zones  of  flames  and  in  the  exhaust. 

Both  gaseous  (propane)  and  liquid  (Jet  A)  fuels  were 
used. 


Fig.  23.  Cross-section  through  combustor  test  tunnel 
indicating  the  CARS  transmitter  on  the  left  and  CARS 
receiving  optics  on  the  right.  Abbreviations  are 
defined  in  the  text. 

The  experimental  scheme  Is  shown  in  cross-section 
in  Fig.  23.  The  pump  and  Stokes  lasers,  along  with 
the  requisite  optics,  are  contained  on  a  one  by  two 
meter  optical  pallet  placed  near  the  test  tunnel. 
Lenses  focus  and  collect  the  BOXCARS  beams  through 
the  tunnel  windows.  It  should  be  noted  that  the  flame 
in  the  test  combustor  is  so  luminous  that  most  con¬ 
ventional  optical  methods  would  fail.  The  emergent 
laser  beams  are  manipulated  In  a  receiver  which 
spectrally  disperses  the  beams,  and  traps  the 


806 


CARS  thermometry  in  reacting  systems 


unwanted  and  ^  f requnci ■ - .  A  reference  PMT 
detects  a  small  fraction  of  the  CARS  signal  and 
serves  to  monitor  the  optical  alignment  of  the  input 
optics.  The  CARS  beam  Is  focussed  into  a  fiber 
optic  (20  m  long)  and  piped  out  to  a  control  room. 

The  control  room  furnishes  a  much  quieter  environment 
for  delicate  instrumentation,  such  as  the  spectrometer, 
and  the  optical  multichannel  analyzer  which  is  subject 
to  microphonics. 

As  examples  of  CARS  measurements  from  this 
experiment.  Fig.  24  displays  averaged  N2  CARS  system 
at  two  different  locations  downstream  of  the  burner 
exit  face.  These  CARS  spectra  correspond  to  100-150 
laser  pulses,  or  to  a  10-15  second  average.  The 
spectrum  at  x  ■  b  cm  was  made  in  the  fuel  spray, 
hence,  the  relatively  cool  temperature.  The  second 
location,  at  39  cm,  is  in  the  flame,  and  is  much 
hotter.  Figure  2.5  demonstrates  that  temperature 
can  be  determined  in  a  highly  luminous,  noisy 
combustion  system,  from  a  single,  10“®  second,  pulse 
spectrum.  Comparison  is  made  with  a  130  pulse 
average  spectrum.  The  single  pulse  spectrum  displays 
photon  statistical  noise  (shot  noise);  however,  it 
is  of  sufficiently  good  quality  to  allow  instanta¬ 
neous  temperature  measurement.  The  single  shot 
measurement  capability,  together  with  the  high 
spatial  resolution  of  BOXCARS,  fulfills  an  obvious 
need  for  studies  in  highly  turbulent,  inhomogeneous 
combustion  systems. 


I’TRO  CARS  Laboratory,  in  East  Hartford, 
a  multipurpose  portable  CARS  apparatus  is  being 
assembled  for  use  at  the  Government  Products  Division 
(CPD)  of  Pratt  &  Whitney  Aircraft.  When  complete, 
measurements  will  be  made  in  the  exhaust  of  a 
production  run  1130  (a  modified  F-100)  gas  turbine 
engine.  In  addition  to  N2  temperature  measurements, 
the  concentration  of  water,  carbon  dioxide,  and  oxygen 
will  be  determined,  as  well  as  some  measure  o)  total 
unburned  hydrocarbons.  These  measurements  are 
scheduled  to  take  place  before  the  end  of  1981. 

The  application  of  CARS  diagnostics  to  both 
research  scale  and  production  line  combustion  devices 
is  increasing  rapidly.  Academic,  government,  and 
industrial  research  laboratories  are,  or  soon  will 
be,  performing  measurements  in  internal  combustion 
engines,  gas  turbine  combustors,  furnaces,  chemical 
process  streams,  and  propellant  burning  strands. 

The  equipment  necessary  for  CARS  is  commercially 
available,  and  special  methods  for  employment  near 
noisy  environments  can  be  developed,  as  illustrated 
by  the  previous  discussion.  The  capabilities  of  the 
CARS  technique  will  be  further  improved  by  experi¬ 
mental  modification,  and  more  importantly,  from 
advances  in  computer  data  processing  (including 
better  modelling  of  spectra).  Considering  the  great 
amount  of  knowledge  that  conventional  diagnostic 
methods  have  provided,  CARS  should  increase  our 
understanding  of  fundamental  and  applied  combustion 
processes  even  further. 


x  =  6cm  T  =  900K  x=39cm  T  =1500K 


Frequency - - 

0.586  cm"1  /  dot 


Fig.  24.  Spatial  variation  of  temperature  from  averaged 
CARS  spectra  of  N2  in  swirl  burner  with  Jet  A  fuel. 
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A  theoretical  model  of  CARS  intensities  based  on  Gordon’s  rotational  diffusion  model  has  been  found  to  give  good  agree¬ 
ment  with  300  K  Nj  spectra  over  the  pressure  range  one  to  100  atmospheres.  The  model  requires  only  that  the  isolated  pres¬ 
sure-broadened  linewidths  be  known  and  greatly  simplifies  the  task  of  calculating  CARS  spectra  for  overlapping  transitions. 


Coherent  anti-Stokes  Raman  spectroscopy  (CARS) 
has  come  to  prominence  in  recent  years  because  of  its 
promise  as  a  versatile  diagnostic  tool  [1-4].  The  appli¬ 
cation  of  CARS  to  combustion  research  in  particular 
is  an  area  of  great  current  interest,  with  the  technique 
proving  useful  for  thermometry  and  species  concentra¬ 
tion  measurements  [5-8],  The  ability  to  model  CARS 
spectra  for  molecules  of  interest  in  combustion  (e.g. 
N2,  CO,  C02,  H20, ...)  is  an  important  aspect  of  the 
process  of  extracting  diagnostic  information  from  ex¬ 
perimental  spectra,  because  of  unique  spectral  inter¬ 
ference  effects  that  preclude  simplified  data  reduction 
in  most  cases  [9] .  One  important  problem  in  the  the¬ 
ory  of  gas-phase  CARS  is  the  effect  of  pressure  on  pre¬ 
dicted  spectra.  Experimentally,  it  is  most  common  to 
probe  vibrational  Q-b ranch  transitions;  for  most  of  the 
molecules  of  interest,  the  splitting  of  adjacent  lines 
due  to  vibration-  rotation  interaction  is  relatively  small, 
giving  rise  to  appreciable  line  overlap  even  at  modest 
gas  pressures.  When  there  is  strong  line  overlap,  the 
overall  band  contour  is  no  longer  adequately  described 
in  terms  of  isolated  lines  undergoing  ordinary  pressure 
broadening.  It  is  then  necessary  to  account  for  the 
phenomenon  of  collisional  narrowing,  which  causes  a 
coalescence  or  collapse  of  adjacent  transitions  toward 
a  frequency  center  of  gravity  determined  by  the  most 


populous  rotational  state  [10-14],  This  effect  is  anal¬ 
ogous  to  narrow  bandwidth  detection  of  an  emitter 
which  is  rapidly  switching  between  discrete  frequen¬ 
cies  with  no  change  of  phase.  Diatomic  ^-branch  fre¬ 
quencies  are  separated  by  2 aeJ,  where  <*„  is  the  vibra¬ 
tion-rotation  interaction  constant.  If  the  mean  time 
between  inelastic  molecular  collisions  is  of  the  order 
of  (2ae/)_1 ,  the  transition  frequency  of  the  molecule 
will  appear  to  switch  rapidly  during  the  radiation  pro¬ 
cess,  and  the  detector  will  record  an  averaged  signal 
in  which  the  individual  Q-branch  frequencies  have 
merged  or  blended  together. 

The  pressure-induced  narrowing  effect  has  been 
modelled  for  CARS,  and  good  agreement  with  experi¬ 
ment  has  been  obtained  for  300  K  N2  over  the  pres¬ 
sure  range  one  to  1 00  atm  [  1 5] .  As  is  well  known 
from  theoretical  studies  of  overlapping  spectral  lines, 
it  is  necessary  to  specify  not  only  the  isolated  pres¬ 
sure-broadened  linewidths,  but  also  an  off-diagonal 
transition  rate  matrix  for  line  amplitudes  whose  ele¬ 
ments  sensibly  give  the  rates  of  collisional  energy 
transfer  between  rotational  energy  levels  [10-14]. 
Unless  these  off-diagonal  linewidth  parameters  can  be 
measured  or  calculated  from  first  principles,  however, 
unknown  parameters  will  have  to  be  introduced  into 
the  calculation.  Also,  these  theoretical  approaches  re- 
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quire  at  each  frequency  the  inversion  of  a  complex  ma¬ 
trix  (sometimes  called  the  “G-matrix”)  (13,14]  result¬ 
ing  in  cumbersome  calculation  for  conditions  requir¬ 
ing  the  inclusion  of  large  numbers  of  rotational  states. 

We  report  here  a  simpler  calculation  of  gaseous  N2 
CARS  spectra  based  on  Gordon’s  extended  7-diffusion 
model  of  rotational  motion  [16,17],  The  analysis 
draws  on  Brueck’s  [18]  application  of  Gordon’s  model 
to  the  lineshapes  of  two-photon  resonances  in  liquids. 
The  rotational  diffusion  theory  has  primarily  been  ap¬ 
plied  to  liquid  phase  infrared  and  Raman  spectra  [16, 
17];  however,  because  it  is  based  on  a  gas-like  rather 
than  a  structure-dominated  view  of  the  liquid  state, 
it  can  with  equal  validity  be  applied  to  the  gaseous 
state.  It  is  an  impact  theory  in  which  the  times  during 
which  intermolecular  torques  act  on  a  molecule  are 
assumed  to  be  short  compared  to  the  time  for  free  ro¬ 
tation  of  the  molecule,  with  free  rotation  assumed  be¬ 
tween  collisions.  Strong  collisions  are  assumed  in 
which  both  the  magnitude  and  orientation  of  the  rota¬ 
tional  angular  momentum  are  randomized  with  each 
collision.  Thus,  the  values  of  the  rotational  quantum 
numbers  after  a  collision  are  unrelated  to  those  before, 
with  the  magnitude  of  the  angular  momentum  assumed 
to  be  redistributed  in  accordance  with  the  Boltzmann 
distribution. 

Vibrational  CARS  signals  for  the  molecules  of  inter¬ 
est  are  generated  mainly  by  the  isotropic  part  a  of  the 
polarizability  tensor,  which  is  sensitive  only  to  the 
magnitude  of  the  angular  momentum  and  not  to  its 
orientation.  Brueck  has  derived  the  third-order  electric 
susceptibility  in  these  approximations  using  a  density 
matrix  formulation  which  properly  accounts  for  vibra¬ 
tion-rotation  interaction.  The  off-diagonal  vibrational 
density  matrix  element  p10  (the  temperature  is  as¬ 
sumed  to  be  so  low  that  the  influence  of  hot  bands 
can  be  neglected)  is  assumed  to  be  a  function  of  the 
rotational  angular  frequency  wr  as  well  as  time,  and 
its  equation  of  motion 

3Plo(«r-  0/9f  =  -0/*)  [H0  +  h(t),  P]  io 

~  fv  Pl°  “  Tj(Lr)  ^Pl°  “-^iXpio*}  O) 

is  solved  to  first  order  with  the  perturbation  interac¬ 
tion  energy  h(t)  =  -  \<xE2  and  a  collision  term  that 
randomizes  the  distribution  of  angular  frequency.  In 
eq.  (1)  Hq  is  the  unperturbed  molecular  hamiltonian; 


Tv  is  a  characteristic  time  for  vibrational  processes; 
r  j(o>j)  is  a  characteristic  time  between  angular  momen¬ 
tum  randomizing  collisions;  /(cor)  is  the  thermal  dis¬ 
tribution  of  angular  velocity,  and  <x>  =  /  doJr  x(cor). 

The  mixing  of  fields  with  frequencies  Wj  and  co2  by 
the  third-order  susceptibility  will  give  rise  to  a  CARS 
frequency  component  in  the  macroscopic  polarization 
p(3)  given  by 

/>(3)  =  Ara,0<p,0>£1  ,  (2) 

where  Vis  the  number  density  of  Raman-active  mo¬ 
lecules,  and  from  this  the  vibrationally  resonant  third 
order  susceptibility  may  be  defined  as 

Xof  =  (A'tt120/ft)</(wr)/(w1  -  co2  -  w01(wr)~  i/Tj(u)r)> 

X  {1  +  i<(/(cor)/rJ(cor)] /(ol>,  -  to2  -  co01(tot) 

-  i/Tj(wr)>}_1  ,  (3) 

where  to01(u>r)  is  the  Q-branch  transition  frequency; 
the  vibrational  broadening  (7”v)  term  has  been  neglect¬ 
ed,  and  the  dependence  of  the  angular  momentum  ran¬ 
domizing  rate  on  wr  (rotational  quantum  number)  has 
been  retained.  The  Boltzmann  population  factor /(cor) 
in  the  numerator  should  strictly  be  the  population  dif¬ 
ference  between  the  upper  and  lower  Q-branch  energy 
levels,  but  the  difference  between  the  two  quantities 
will  be  small  for  modest  temperatures.  If  the  tempera¬ 
ture  is  sufficiently  high  that  the  populations  of  higher 
vibrational  states  cannot  be  neglected,  then  the  reso¬ 
nant  contribution  to  x*3*  can  be  represented  as  a  sum 
of  contributions  from  all  significant  bands; 

*(3)=?e.  w 

with 

(3)  M&+1  <A/>r)AP„K)> 

xv,vl-  ft  l+i<[/(wr)/r,(wr)] /©„(«,)>’  (4) 

where  A /^(w,)  is  now  the  population  difference  fac¬ 
tor  with  /;(wr)  normalized  to  unity  for  all  vibrational 
and  rotational  mates,  and  Dv( cor)  =  Wj  -  w2 
-  wu  u+i(cjr)  -  i/Tj(«r).  The  above  formulation  is 
based  on  the  assumptions  that  vibrationally  inelastic 
energy  transfer  will  be  relatively  slow,  and  that  neither 
the  rotational  Boltzmann  factor  /nor  the  angular  mo¬ 
mentum  randomizing  frequency  (tj)-1  will  have  a  sig¬ 
nificant  dependence  on  vibrational  state.  In  the  ab- 
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sence  of  vibrational  effects,  the  angular  momentum  re¬ 
distribution  frequency  ri  (cj,)"1  is  equivalent  to  the 
pressure-broadened  halfwidth.  Similar  formulations  of 
the  Gordon  rotational  diffusion  theory  have  been  given 
by  Lynden-Bell  [  19] ,  Eagles  and  McClung  [20] ,  and 
Temkin  and  Burshstein  [21]. 

Given  x*3\  the  predicted  CARS  spectral  density 
is  given  by  the  convolution  of  |x(3*l2  with  the  spectral 
profiles  of  the  pump  and  Stokes  sources  [9,22].  Nu¬ 
merical  calculations  were  performed  using  eq.  (3)  to 
simulate  the  pressure  variation  of  the  300  K  N2  CARS 
spectrum  reported  in  [15].  For  these  calculations  the 
N2-N2  Raman  linewidths  measured  by  Rahn  and 
Owyoung  [23]  were  scaled  linearly  with  pressure  to 
give  rj(cjr)_1 .  Because  these  widths  have  a  significant 
dependence  on  rotational  quantum  number,  it  was  not 
possible  to  obtain  an  analytic  expression  for  x(3)  as 
was  done  by  Brueck,  and  the  integrations  over  to,  re¬ 
quired  in  eq.  (3)  were  performed  using  explicit  sum¬ 
mations  over  the  discrete  rotational  states  [21].  Com¬ 
parison  of  the  rotational  diffusion  model  with  the  ex¬ 
perimental  data  of  [1 5]  presents  a  good  test  for  the 
theory  because,  with  the  transition  linewidths  r  j(or)~ 1 
known,  the  problem  is  completely  specified,  and  there 
are  no  adjustable  parameters. 

A  -'omparison  of  the  signature  bandwidths  predict¬ 
ed  by  eq.  (3)  with  the  experimental  data  of  [15]  is 
shown  in  fig.  1 .  As  explained  in  [  1 5  ] ,  the  various  sym¬ 
bols  denote  the  results  of  different  series  of  broadband 
experiments  taken  over  a  period  of  several  weeks,  and 
the  assignment  of  bandwidth  to  each  signature  is  pi  one 
to  a  subjective  error  in  reading  of  peak  height  and  full 
width  at  half  maximum  that  is  estimated  to  be 
±0.1  cm-1 .  Bearing  this  in  mind,  the  comparison  be¬ 
tween  experiment  and  theory  is  seen  to  be  quite  satis¬ 
factory,  with  agreement  to  within  the  experimental 
error.  The  theory  predicts  a  slight  minimum  in  the  vi¬ 
cinity  of  2  atmospheres,  a  gradual  increase  in  the 
bandwidth  to  about  20  atmospheres,  and  then  a  steady 
contraction  out  tq  100  atmospheres.  The  theory  does 
tend  to  the  low  side  of  the  experimental  data  at  the 
higher  pressures;  this  may  indicate  the  presence  of  a 
small  vibrational  dephasing  contribution.  Also  shown 
in  the  figure  is  the  bandwidth  predicted  on  the  basis 
of  isolated  line  theory,  that  is,  neglecting  the  <  )  term 
in  the  denominator  of  eq.  (3).  At  100  atm,  the  predic¬ 
ted  isolated  line. bandwidth  exceed  10  cm-1 ,  empha¬ 
sizing  the  very  significant  role  that  collisional  narrow- 


Fig.  1 .  Bandwidth  (FWHH)  of  300  K  N2  CARS  signature  ver¬ 
sus  pressure.  Symbols  denote  experimental  points  [15);  solid 
lines  denote  theoretical  calculations  for  an  assumed  overall  re¬ 
solution  of  0.4  cm'1 . 

ing  plays  at  higher  pressures. 

Scatter  in  the  data  presented  in  fig.  1  precludes  a 
detailed  verification  of  the  predicted  bandwidth  varia¬ 
tion  at  low  pressures.  However,  the  initial  contraction 
of  the  bandwidth  at  2-4  atmospheres,  with  a  subse¬ 
quent  increase  out  to  approximately  20  atmospheres, 
does  seem  to  be  confirmed  qualitatively  by  the  experi¬ 
ments  of  Roland  and  Steele  [24] ,  in  which  the  spec¬ 
trum  of  300  K  N2  was  examing  using  tuned,  narrow- 
band  sources  for  pressures  up  to  33  atm. 

In  view  of  the  fact  that  the  theoretical  spectra  con¬ 
tain  no  adjustable  parameters,  the  overall  agreement 
exhibited  in  fig.  1  is  gratifying.  Use  of  eq.  (3)  to  model 
collisional  narrowing  effects  in  CARS  spectra  offers 
great  computational  advantages,  because  the  calcula¬ 
tion  of  X<3>  is  only  slightly  more  involved  than  in  the 
isolated  line  approximation  and  therefore  requires 
practically  no  extra  computation  expense.  It  will  be 
interesting  to  see  whether  the  rotational  diffusion  the¬ 
ory  is  also  applicable  to  the  gas-phase  CARS  spectra  of 
other  molecules  of  interest;  for  molecules  such  as  H20 
and  C02  in  particular,  whose  spectra  are  comprised  of 
very  large  numbers  of  transitions,  the  rotational  diffu¬ 
sion  theory  would  make  possible  narrowing  calcula¬ 
tions  that  would  otherwise  be  very  costly. 
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Some  Applications  of  Gas  Phase  CARS  Spectroscopy 


J.  F.  Verdieck,  R.  J.  Hall,  J.  A.  Shirley,  and  A.  C.  Eckbreth 

United  Technologies  Research  Center,  East  Hartford,  CT  06108 

CARS  (coherent  anti-Stokes  Raman  spectroscopy)  is  a 
unique  nonlinear  optical  technique  made  possible  with  high 
power  lasers.  It  is  a  spectroscopic  method  which  obeys  Raman 
selection  rules  but  is  a  much  stronger  effect  than  conventional 
Raman  spectroscopy.  Equally  important,  the  CARS  radiation 
emerges  from  the  sample  as  a  laser-like  beam  which  may  be 
completely  collected,  thus  providing  good  discrimination 
against  incoherent  interfering  radiation.  For  these  reasons, 
CARS  is  finding  numerous  applications  involving  small  vol¬ 
ume,  highly  dilute  samples  such  as  biochemical  substances 
in  solution,  or  gas  phase  samples.  The  application  of  the  CARS 
spectroscopic  technique  as  a  diagnostic  method  for  measuring 
temperature  and  species  concentrations  in  combustion  sys¬ 
tems  of  various  types  will  be  described  in  this  paper. 

Problems  in  Combustion  Diagnostics 

Combustion  diagnostics  is  a  fascinating  and  rewarding,  yet 
very  challenging,  area  for  spectroscopic  analysis.  Although 
non-spectroscopic  methods,  such  as  mass  spectrometry  and 
gas  chromatography,  have  contributed  significantly  to  our 
understanding  of  combustion  processes,  spectroscopic 
methods  offer  several  advantages  over  probe  sampling 
methods.  Some  of  the  advantages  of  spectroscopic  techniques 
are:  (1)  remote  operation,  (2)  rapid  measurement,  (3)  com¬ 
bustion  process  is  not  disturbed,  (4)  good  spatial  precision  is 
possible,  and  (5)  very  high  temperature  capability.  Each  of 
these  advantages  will  be  illustrated  by  several  CARS  appli¬ 
cations  described  below. 

The  advantages  just  listed  apply  equally  well  for  CARS 
spectroscopic  measurement  of  temperature.  It  Bhould  be 
emphasized,  that  in  combustion,  the  single  most  important 
parameter  to  measure  is  the  temperature.  Given  the  initial 
conditions  and  the  temperature  at  a  point  in  a  combusting 
system,  the  products,  heat  release,  and  other  parameters  can 
be  estimated  quite  well.  For  this  reason,  highly  accurate, 
spatially  precise  measurements  of  temperature  are  a  necessity 
in  understanding  and  controlling  combustion  processes.  It  is 
noted  that  thermocouple  probes  are  limited  to  temperatures 
below  ~2500  K.  In  contrast,  spectroscopic  methods  have  much 
higher  limits,  up  to  several  thousand  degrees  for  the  case  of 
nitrogen  vibrational  thermometry,  as  will  be  explained  in 
detail.  Moreover,  spectroscopic  temperature  measurement 
is  more  accurate  than  thermocouple  probe  thermometry  be- 
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Figure  1.  Basic  Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  Process. 


cause  of  the  corrections  which  must  be  applied  to  thermo¬ 
couples,  namely,  radiative  and  conduction  corrections. 

To  gain  an  appreciation  of  why  a  powerful,  wide-ranging 
technique  such  as  CARS  is  required  for  combustion  diag¬ 
nostics,  it  is  useful  to  list  some  examples  of  combustion  sys¬ 
tems  which  are  investigated.  Combustion  devices  can  range 
from  simple  candle  or  Bunsen  burner  type  flames  to  controlled 
detonations  which  occur  in  rockets  or  gun  barrels.  In  between 
these  extremes  are  home  heating  and  industrial  furnaces, 
gasoline  and  diesel  engines,  and  gas  turbine  engines  (both 
airborne  and  ground-based).  These  extremes  represent  very 
wide  ranges  of  temperature  and  pressure.  In  an  energy-scarce 
society,  equally  concerned  about  its  environment,  it  is  es¬ 
sential  to  design  our  combustion  systems  to  be  energy  efficient 
and  as  nearly  pollution  free  as  we  can.  Reliable  diagnostic 
methods  are  required  to  gain  the  knowledge  necessary  to 
achieve  these  goals,  and  CARS  can  help  attain  that  knowledge. 
For  more  detailed  discussions  of  these  types  of  CARS  appli¬ 
cations,  the  reader  may  check  references  (1-3). 

CARS  Essentials:  Experimental 

An  introduction  to  the  CARS  technique  is  made  by  first 
describing  the  CARS  experiment  and  the  apparatus  required. 
The  CARS  theory  will  be  described  in  the  following  section. 
CARS  is  a  coherent,  nonlinear-optical  process  in  which  light 
waves  of  different  frequencies,  <*>i  and  u>2,  are  mixed  together 
in  a  material  medium  (liquid,  gas,  or  solid)  to  generate  a  new 
frequency,  W3,  which  appears  in  a  coherent,  laser-like  beam. 
This  is  shown  by  the  simple  sketch  of  Figure  1.  When  the 
frequency  difference  between  the  two  input  frequencies  cor¬ 
responds  to  a  Raman-active  molecular  vibration  or  rotation 
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Figure  2  CARS  experimental  arrangement  lor  flame  measurements 


(or  electronic  transition),  then  the  CARS  effect  is  strong  and 
easily  detected.  Because  the  CARS  process  is  a  coherent 
nonlinear-optical  process,  laser  beams  are  required  to  provide 
the  proper  phase  of  waves  in  space  and  time,  and  the  intensity 
to  enhance  the  nonlinear  aspect  of  the  process.  Usually, 
high-peak  power  pulsed  lasers  are  employed  to  generate 
CARS,  but  lower  power  cw  lasers  have  been  used  in  a  few  cases 
(4). 

At  the  United  Technologies  Research  Center  (UTRC), 
CARS  experiments  are  based  upon  a  frequency-doubled 
needy nium:Y AO  laser,  hereafter  simply  referred  to  as  2X 
Nd.YAG.  The  frequency-doubled  output  occurs  in  the  visible 
at  532  nm  and  is  used  for  one  of  the  CARS  input  frequencies, 
uii.  The  second  CARS  frequency,  w2,  is  generated  from  a  dye 
laser  driven  bv  a  portion  of  the  2XNd:YAG  beam.  Thus,  the 
dye  laser  pulses  are  nearly  synchronous  with  the  2XNd:YAG 
pulses.  The  2XNd:YAG  pulses  are  ~10  nanoseconds  (10~8sec) 
in  duration  and  occur  at  a  repetition  rate  of  10  Hz.  The  high 
repetition  rate  is  advantageous  because  it  permits  high  data 
rate  acquisition,  and  for  adjusting  optics  for  critical  alignment 
of  the  CARS  input  beams.  The  peak  power  in  a  2XNd:YAG 
pulse  is  — 10  MW;  the  average  power  is  ~2  W. 

The  arrangement  of  the  lasers  is  shown  in  Figure  2.  A  por¬ 
tion  of  the  2XNd:YAG  radiation  (532  nm)  is  picked  off  with 
a  beamsplitter  to  pump  the  oscillator  and  amplifier  stages  of 
the  dye  laser.  The  532  nm  radiation  is  ideally  suited  to 
pumping  laser  dyes  in  the  550-700  nm  range,  particularly 
rhodamine  type  dyes  which  can  exhibit  conversion  efficiencies 
of  more  than  30%.  The  dye  lasers  are  of  quite  simple  con¬ 
struction,  utilizing  either  spectrophotometer  cells  or  home¬ 
made  dye  cells.  The  dye  solution  flows  through  the  cells  to 
prevent  heat  buildup  and  subsequent  optical  distortion  of  the 
laser  beam.  The  532  nm  is  directed  nearly  parallel  to  the  dye 
laser  optical  axis  (defined  by  the  dye  laser  mirrors  DM;  and 
DM2).  This  so-called  end  pumping  (or  axial  pumping)  yields 
good  dye  laser  beam  quality  and  good  conversion  efficiency. 
The  purpose  of  the  Galilean  telescope  (negative  lens-positive 
lens  pair,  T)  in  the  dye  laser  beam  is  to  insure  that  the  dye 
beam  will  focus  at  the  same  region  as  the  532  nm  beam. 

The  two  laser  beams  are  combined  on  a  dichroic  (two-color) 
mirror,  D|.  In  this  case  the  dye  laser  frequency,  w2,  also  termed 
the  Stokes  frequency,  passes  through  the  dichroic  mirror  while 
the  532  nm  beam,  or  pump  frequency,  oi|,  is  reflected  by  it. 
This  type  of  configuration  is  called  collinear  phase  matching 
and  will  be  compared  to  other  beam  configurations  below.  The 
collinear  beams  are  focussed  (by  lens  L)  through  the  flame  or 
combustion  system  of  interest  and  recollimated  by  a  similar 
lens.  The  newly  generated  CARS  frequency,  (  =  2wj  -  w2), 
can  be  separated  from  the  u>i  and  w2  frequencies  by  use  of 
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Figure  3.  Generation  ot  CARS  spectra 


(a)  Ganaral  (b!  Coltlnaar  (c)  BOXCARS 


Figure  A.  CARS  phase-matching  techniques. 

another  dichroic,  D2,  which  transmits  blue  light,  but  not  green 
(532  nm),  or  dye  laser  light  (yellow  to  red).  Dispersing  prisms 
can  also  be  used  to  isolate  the  CARS  beam.  The  C ARS  beam 
is  then  directed  to  a  monochromator  for  spectral  analysis,  from 
which  temperature  and  concentration  information  is  extracted 
as  described  in  the  theoretical  section  below.  A  portion  of  the 
CARS  beam  is  picked  off  by  a  beamsplitter  and  detected  bv 
a  photomultiplier-filter  combination  ( I’MT  l ).  The  detector 
serves  as  a  monitor  of  the  total  (spectrally  wideband)  CARS 
signal.  This  monitor  can  be  used  as  a  normalizing  signal  and 
is  also  used  to  adjust  optical  components  for  maximum  C  ARS 
signal. 

Spectral  analysis  of  the  CARS  signal  is  obtained  by  scan¬ 
ning  with  a  monochromator  and  processing  the  output  pho¬ 
tomultiplier  signal,  FMT-2  in  Figure  2,  with  a  laser  triggered 
electronic  gate,  a  so-called  boxcar  averager.  A  more  useful 
alternative  is  to  employ  an  optical  multichannel  analyzer 
(OMA)™1  which  can  capture  the  entire  CARS  spectrum,  with 
gain,  from  each  laser  shot .  This  information  is  sent  to  a  com¬ 
puter  for  storage,  analysis,  and  displav.  This  latter  technique 
offers  the  potential  for  measuring  temperature  and  concen¬ 
tration  from  a  single,  10  nanosecond,  laser  pulse.  Both 
methods  ol  spectral  analysis  are  in  use  at  UTRC  and  have 
been  described  in  several  references  (3, 5-6).  Several  examples 
of  CARS  spectra  are  given  below. 

It  must  be  mentioned  that  the  method  of  generating  CARS 
spectra  described  here  is  called  "broadband  CARS”  or 
sometimes  “multiplex  CARS.”  It  is  only  one  of  several  dif¬ 
ferent  techniques.  One  other  method,  used  for  high  resolution 
CARS  spectroscopy  (8, 9)  makes  use  of  a  tunable,  narrowband 
dye  laser.  The  CARS  spectrum  is  generated  by  combining  the 
u>i  and  ui2  laser  frequencies,  and  while  tuning  w2,  detect  the 
amplitude  changes  in  the  CARS  signal  (whose  frequency,  ws, 
is  changing).  These  two  methods  are  contrasted  pictorially  in 
Figure  3.  For  combustion  diagnostics,  the  broadband  C ARS 
method  appears  more  suitable  because  it  generates  a  complete 
spectrum  with  each  laser  pulse.  Single-shot  thermometry,  with 
the  aid  of  an  OM  A™  and  computer,  is  then  possible,  as  men¬ 
tioned  previously 

For  the  broadband  CARS  method,  the  dye  laser  is  operated 
wideband,  that  is,  with  appropriate  mirrors  which  generate 
laser  output  which  has  a  large  bandwidth  (full  width  at  half 
maximum)  of  MOO  cm-1,  or  greater.  The  center  of  this  spec¬ 
tral  band  should  be  located  at  <«>■>  =  u>i  -  wham,  where  wram 
corresponds  to  the  desired  Raman  frequency  shift.  The  de¬ 
sired  frequency  is  obtained  by  appropriate  dye  selection  and 
concentration  adjustment. 
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a)  Planar 


Figure  5  Illustration  ot  touted  BOXCARS  geometry. 


Figure  6  Comparison  of  nonlinear  optical  processes. 


The  CABS  effect  requires  that  the  proper  relationship  ex 
ists  between  the  input  frequencies,  u>i,  u>2,  and  the  CARS 
frequency,  a'n,  namely:  W3  =  2u>i  -  w2.  This  is  a  conservation 
of  energy  requirement.  A  similar  conservation  of  linear  mo¬ 
mentum  requirement  demands  that  the  laser  beams  travel  in 
certain  directions  with  respect  to  each  other.  Satisfying  these 
directional  requirements  is  known  as  phase-matching.  The 
relationships  which  must  be  satisfied  for  maximum  CARS 
generation  are  shown  in  Figure  4.  The  wave-vector  k,  has 
magnitude  w,n,/c  (n,  is  the  refractive  index  at  frequency  «,■), 
and  direction  given  by  the  propagation  direction  of  the  laser 
beam.  Note  that  the  phase-matching  condition  is  a  vector  sum 
as  illustrated  in  Figure  4a.  In  gases,  the  refractive  index  varies 
little  with  frequency  which  means  that  the  collinear  condition, 
Figure  4b,  satisfies  the  phase  matched  condition.  Although 
it  is  the  easiest  configuration  to  implement  and  yields  high 
intensity  CARS  signals,  collinear  operation  can  suffer  from 
poor  spatial  resolution.  The  reason  for  this  is  that  the  u>i  and 
u>2  input  beams  interact  and  generate  CARS  over  a  consider¬ 
able  sample  length,  dependent  upon  the  focal  length  of  the 
lens  used.  For  example,  a  10-cm  lens  yields  a  sample  length 
of  ~0.07  cm,  whereas  a  100-cm  lens  (which  might  be  manda¬ 
tory  for  some  cases)  yields  a  sample  length  of  7.3  cm  (3). 
Clearly  this  length  is  much  too  large  and  could  sample  both 
hot  and  cold  regions  of  a  combusting  system. 

Recognizing  this,  Eckbreth  19)  has  shown  how  a  phase- 
matched,  crossed-beam  configuration  leads  to  excellent  spatial 
resolution,  because  CARS  is  generated  only  from  the  region 
of  beam  crossing.  The  method  has  been  termed  BOXCARS, 
Figure  4c,  from  the  shape  of  the  phase  matching  diagram. 
Resolution  of  a  millimeter  or  less  is  easily  achieved,  even  for 
relatively  shallow  angles  between  input  beams.  BOXCARS 
has  been  employed  in  this  laboratory  t>  map  out  the  tem¬ 
perature  distribution,  including  both  the  axial  and  radial 
variation,  in  a  conventional  diffusion  flame  (10).  It  is  not 
mandatory  to  confine  phase  matching  to  a  plane.  One  type  of 
three-dimensional  phase  matching,  illustrated  in  Figure  5,  is 
called  "folded  BOXCARS”  for  obvious  reasons.  An  out¬ 
standing  advantage  of  folded  BOXCARS  is  that  the  013  CARS 
beam  is  completely  separated  from  the  u>i  and  w2  beams 
without  using  a  dichroic  mirror  or  a  prism.  This  allows  one  to 
examine  very  small  Haman  shifts  such  as  found  in  rotational 


CARS  spectroscopy.  Folded  BOXCARS  has  been  utilized  at 
1  TRC  for  nitrogen  and  air  rotational  CARS  spectra  in  flames 
(11).  Several  other  phase  matching  schemes  for  CARS  exist 
(12,  13)  which  are  not  discussed  here  because  of  space  limi¬ 
tations. 


CARS  Essentials:  Theoretical 

In  this  brief  section  an  elementary  presentation  of  CARS 
theory  is  given  which  provides  a  basis  for  temperature  and 
concentration  measurements.  More  complete  descriptions  of 
the  theory  may  be  found  in  references  14-16. 

Perhaps  an  appropriate  place  to  begin  is  to  draw  a  com¬ 
parison  with  a  nonlinear  optical  process  which  may  be  famil 
iar,  na  mely  second -harmonic  generation  (frequency  doubling). 
This  second-order  (because  it  is  quadratic  in  laser  intensity) 
effect  converts  laser  radiation  of  frequency  w  into  radiation 
of  frequency  2w.  One  way  of  explaining  this  effect  is  to  state 
that  two  photons  are  absorbed  simultaneously  and  a  single 
photon  of  twice  the  initial  energy  is  emitted.  The  absorption 
and  emission  occurs  through  "virtual  energy  states;”  real 
energy  states  are  not  directly  involved.  Figure  6a  presents  this 
explanation  pictoriaily.  Figure  6b  shows  the  more  complicated 
third-order  process  of  CARS.  CARS  involves  the  mixing  of 
three  light  waves  of  frequency  wlt  oj],  and  w2  to  generate  a 
fourth  beam  at  frequency  o> 3.  The  intensity  of  the  CARS  beam 
is  given  by 


/3  =  Ar/,2/2|x‘:'’l2^2  (■) 

The  third-order  dependence  of  laser  intensity  is  readily  seen 
(/1,  /2,  and  13  are  the  intensities  at  frequencies,  coa,  w2,  and  W3, 
respectively).  x(3>  is  the  third-order  nonlinear  susceptibility, 
a  macroscopic  parameter  which  relates  the  electric  fields  of 
the  incoming  light  waves  to  the  polarizations  (charge  distor¬ 
tion)  in  the  macroscopic  medium.  x(:u  depends  on  both  nu¬ 
clear  and  electronic  motions  of  the  molecules  which  compose 
the  medium  (liquid,  solid,  or  gas).  The  term  z  is  an  interaction 
length  over  which  the  nonlinear  process  is  strongly  generated. 
K  is  a  constant. 

X(3>  can  be  separated  into  the  terms  shown  in  eqn.  2. 

X13’  “  x'  +  ix"  +  xnr  (2) 


The  first  two  terms  account  for  the  resonant  portion  of  the 
susceptibility  which  relates  to  nuclear  motion;  hence,  is  res¬ 
onant  with  certain  (Raman  active)  vibrational-rotational 
frequencies.  xnr  accounts  for  the  electronic  contribution  from 
all  molecular  species  present. 

The  resonant  susceptibility  depends  upon  individual  mo¬ 
lecular  transitions  in  a  complicated  fashion  given  bv  eqn. 
(3), 


X‘  +  <x'  =  ^4IVA/|t£i 


U>2 


1 


dilli  l2Aw,  -  if, 


(3) 


K2  is  a  constant,  oi2  is  the  Stokes  frequency,  N  the  molecular 
number  density,  A;  is  a  population  difference  factor  (essen¬ 
tially  the  Boltzmann  factor  from  which  temperature  is  ob¬ 
tained),  (dir/dibj  is  the  differential  Raman  scattering  cross- 
section,  Au>  is  the  detuning  from  resonance  (i.e.,  wj  -  u >■>  - 
u>R*m  =  0  at  resonance),  and  finally,  V,  is  the  Raman  linewidth 
of  the  y  th  transition.  Several  important  points  should  be  made 
about  eqn.  3.  As  mentioned  parenthetically,  temperature  in¬ 
formation  is  contained  in  the  A,  term;  concentration  from  the 
number  density  N.  In  order  to  extract  this  information,  it  is 
mandatory  to  interpret  experimental  CARS  spectra  by  direct 
comparison  with  computer  generated  model  spectra.  This 
further  implies  that  all  the  molecular  data,  Raman  cross- 
sections,  linewidths  (including  temperature  and  quantum 
number  dependences)  should  be  known.  Note  that  w...  n  Aw, 
=  0,  the  resonant  portion  of  x<3’  is  measured;  hence,  one  is 
dealing  with  ail  the  terms  of  eqn.  (4) 

lx3|J  »  lx’)2  +  2x'x,r  4  <\"rR  +  <\'l2  HI 
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Figure  7  Limiting  cases  ot  CARS  spectra 


This  equation  is  best  understood  by  referring  to  Figure  7 
where  the  strong  line  and  weak  line  limits  lor  CARS  spectra 
are  depicted.  In  the  strong  line  limit,  \  is  negligible  com 
pared  to  the  resonant  portion.  \  and  the  spectrum  reflects 
the  character  of  l  y '•■>)-'  In  the  weak  limit  y  cannot  be  tie 
fleeted  and  contributes  a  background  baseline  upon  which 
the  resonant  CARS  contribution  appears  as  a  modulation.  The 
limit  where  the  modulation  can  just  be  observed  gives  the  limit 
of  detection  of  the  desired  species  in  the  particular  environ¬ 
ment  which  yields  that  non-resonant  background  (contributed 
by  all  species  present  ). 

The  non-resonant  background  is  generally  considered  a 
nuisance,  and  techniques  have  been  designed  to  suppress  or 
eliminate  it  (discussed  below  I.  However,  in  many  cases  of  in¬ 
terest  in  combustion,  the  shape  of  the  CARS  signal  against  the 
Xnr  background  provides  for  easily  interpreted  concentration 
measurements,  without  the  need  of  intensity  measurements. 
Examples  of  this  method  will  be  shown  later. 

CARS  Spectra:  Temperature  Measurements 

In  this  section.  CARS  spectra,  both  experimental  and 
computer  model  predicted,  will  he  presented  for  several  dif¬ 
ferent  molecules.  The  emphasis  is  given  to  nitrogen  CARS 
spectra  because  nitrogen  is  present,  in  quantity,  in  air-sup¬ 
ported  combustion.  Moreover,  the  necessary  spectroscopic 
constants  are  known  quite  accurately  for  nitrogen.  The  change 
in  the  spectral  shape  with  temperature  will  be  demonstrated 
for  nitrogen  and  other  molecules  as  well.  The  effect  of  high 
pressure  on  CARS  spectra  will  be  treated  briefly  and  illus¬ 
trated  with  experimental  results.  Applications  of  CARS 
temperature  measurements  in  real  combustion  systems  will 
be  demonstrated  by  an  example  from  this  laboratory. 

Figure  8  shows  a  set  of  computer-generated  spectra  of  ni¬ 
trogen  for  typical  experimental  conditions  over  a  temperature 
range  of  300  to  2100  K,  taken  in  300  K  increments  (17).  The 
300  K  room  temperature  spectrum  is  a  narrow  lin  with  peak 
at  ~21126cm~‘,  exhibiting  no  rotational  structure  at  mod¬ 
erate  resolution.  As  the  temperature  increases,  rotational 
structure  begins  to  appear.  The  C ARS  spectrum  shown  here 
is  for  the  Q-branch  =  0).  r  =  0  to  e  -  1  transition.  At 
about  1200  K  the  hot  band  (r  =  1  to  r  =  2)  transition  begins 
to  appear.  At  1500°.  the  hot  band  has  a  peak  height  of  about 
0.1  of  the  fundamental,  and  reaches  about  0.4  of  the  funda¬ 
mental  at  2100°.  Clearly,  there  is  considerable  change  in  the 
CARS  spectrum  with  temperature,  and  from  this  change  the 
temperature  in  a  combusting  system  could  be  "eyeballed"  to 
about  ±100°.  As  Figure  9  demonstrates,  if  a  least  square 
computer  fit  to  the  experimental  data  is  made,  temperature 
can  be  measured  to  better  than  l,rr  precision.  What  is  shown 
in  Figure  9  is  the  computer  least  squares  fit  (dotted  line)  su¬ 
perimposed  on  an  experimental  CARS  spectrum  (solid  line) 
from  a  premixed  CHj/Air  flame,  whose  temperature  was  in¬ 
dependently  measured  bv  sodium  line  reversal  </8l  The 
agreement  is  excellent. 

CARS  has  been  used  to  determine  temperatures  in  highly- 
sooting  flames  (It/).  These  measurements  are  significant. 


Figure  8  Variation  of  nitrogen  CARS  s  pec  fra  wit*  temperature 


Figure  9  Comparison  ot  experiment  and  theory  for  N?  at  2  t00  k 


Figure  i0  CARS  temperature  measurement  along  mt  o»  flame 
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Figure  11.  CARS  spectra  of  nitrogen  at  1  and  100  atm  :  solid  line;  experimental, 
dotted  line;  theoretical 


Figure  12.  Model  illustration  of  CARS  linewidth  with  and  without  collisions! 
narrowing. 


Figure  13.  Pressure  variation  of  hi(pi  temperature  CARS  spectra. 


Figure  14.  M>»ogsn  CARS  apecva  as  a  function  of  temper ahrs:  modal  pre- 
dkllona. 


because  practical  combusting  systems  often  operate  under 
-  •nfimr  conditions  which  preclude  conventional  optical 
techniques,  such  as  spontaneous  Raman  scattering.  Figure  10 
shows  the  axial  variation  of  temperature  in  a  sooty  propane 
laminar  diffusion  flame  as  found  by  the  BOXCARS  method. 
The  spatial  resolution  in  these  experiments  was  about  0,;t  mm 
by  1.0  mm,  which,  as  the  temperature  profile  demonstrates, 
is  necessary  for  exploring  flame  structure. 

Pressure  Effects  in  Nitrogen  CARS  Thermometry 

In  order  to  achieve  maximum  utility  as  a  diagnostic  tool,  it 
is  necessary  to  extend  CARS  measurements  into  the  high 
pressure  regime.  Kxamples  of  important  practical  combustion 
systems  which  operate  at  elevated  pressures  are  internal 
combustion  engines  (both  gasoline  and  diesel),  gas  turbine 
engines  and  rockets.  The  pressure  encountered  in  such  sys¬ 
tems  ranges  from  a  few  atmospheres  to  well  over  100  atm.  If 
the  theoretical  CARS  analysis  which  has  been  used  to  suc¬ 
cessfully  explain  one  atmosphere  CARS  (lame  spectra  were 
employed  at  very  high  pressure  it  would  predict  that  the 
overall  bandwidth  of  the  spectrum  would  increase  linearly 
with  pressure.  It  would  also  predict  that  the  1  *  2  "Hot"  hand 
would  merge  with  the  fundamental  hand  and  therefore  would 
be  lost  as  a  sensitive  indicator  of  temperature. 

The  effect  of  pressure  on  N_>  CARS  spectra  has  been  in¬ 
vestigated  experimentally  and  theoretically  over  the  pressure 
range  of  one  to  100  atm.  (6).  A  similar  study  has  been  per¬ 
formed  by  Roland  and  Steele  (19)  for  pressures  up  to  00  atm. 
The  experimental  studies  of  Hall,  Verdieck,  and  Eckbreth 
show  that  after  an  initial  narrowing  of  the  Q-branch  band  over 
the  first  five  to  ten  atmospheres,  the  width  of  the  band  re¬ 
mains  constant  with  increasing  pressure  to  100  atm,  the  width 
of  the  band  remains  constant  with  increasing  pressure  to  100 
atm.  Figure  1 1  shows  scanned  CARS  Q-branch  spectra  at  one 
and  100  atm.  The  solid  line  in  Figure  1 1  is  the  experimental 
spectrum  obtained  from  collinear  CARS,  scanned  with  a 
monochromator  with  spectral  width  of  ~0.4  cm-1.  The  dotted 
line  represents  the  theoretical  fit  of  the  computer  model.  The 
collapse  of  the  resolved  rotational  structure  to  the  symmetric 
band  is  termed  “collision-induced  band  narrowing”  or  "mo¬ 
tional  narrowing.”  It  is  a  phenomenon  observed  in  NMR  (20) 
and  in  Raman  spectroscopy  (21).  Figure  12  illustrates  the 
importance  of  collisional  narrowing  through  theoretical 
comparison  of  a  collision-narrowed  line  with  one  calculated 
without  narrowing.  For  the  latter  case  a  very  broad  line,  ~  10 
cm-1,  results.  The  effect  of  pressure  on  high  temperature  has 
also  been  investigated,  some  results  of  which  are  shown  in 
Figure  13  (2).  Experimental  results  are  shown  for  ~1600  K  for 
pressures  of  ~8  and  ~28  atm.  along  with  model  spectra.  It  is 
noted  that  although  the  high  J  rotational  structure  collapses 
between  these  pressures,  the  “hot”  band  persists;  hence,  the 
potential  for  temperature  measurement  still  exists.  High 
pressure,  high  temperature  CARS  spectra  must  be  examined 
more  thoroughly  in  order  to  assess  CARS  diagnostics  for  high 
pressure  applications. 

CARS  Spectra  of  Other  Molecules 

Several  other  molecules  of  interest  to  combustion  studies 
have  been  examined  in  this  laboratory;  CARS  spectra  have 
been  recorded  for  hydrogen,  oxygen,  carbon  monoxide,  water, 
carbon  dioxide,  and  methane  (22-23).  Of  these,  hydrogen  and 
water  CARS  spectra  will  be  discussed  with  regard  to  tem¬ 
perature  measurement. 

Computed  CARS  spectra  of  hydrogen  are  shown  in  Figure 
14  for  four  temperatures.  Because  the  vibrational-rotational 
constant,  a,  is  so  large  in  hydrogen  (a  ”  2.993,  compared  to 
a  “  0.0187  for  nitrogen),  the  spectrum  is  simple  (few  lines), 
and  quite  spread  out  (over  nearly  200  cm-1).  The  alternation 
in  intensity  due  to  nuclear  spin  statistics  is  quite  apparent. 
Recall  that  the  ortho/para  forms  of  hydrogen  bear  statistical 
weights  of  3:1.  Because  population  differences  get  squared  in 
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Figure  1 7  CARS  water  spectm  at  7  73  K 


Figure  16  CARS  mapping  of  temperature  in  a  H2/02  flame. 


a  CARS  spectrum  (as  \  is  squared),  the  difference  in  adjacent 
lines  is  quite  dramatic  in  CARS  spectra.  For  this  reason,  H-. 
should  serve  as  a  good  thermometric  probe  molecule.  Exper¬ 
imental  CARS  spectra  for  Hj  are  shown  in  Figure  15.  as  re¬ 
corded  with  an  OMA'".  The  spectra  were  recorded  from  an 
Hj/Air  diffusion  flame,  with  temperatures  found  from  com¬ 
parison  of  peak  height  ratios  of  Q(l),  Q(3).  and  Q15).  Tem¬ 
perature  profiles  of  a  slot  shaped  H^/Air  diffusion  flame  were 
obtained  from  H2  and  CT  CARS  spectra  (Oj  CARS  spectral 
temperatures  are  obtained  in  a  fashion  similar  to  that  for  N_>) 
and  are  displayed  in  Figure  16.  The  CARS-derived  tempera¬ 
tures  are  compared  with  a  radiation-corrected  thermocouple. 
The  agreement  is  quite  good.  A  detailed  discussion  of  the 
hydrogen  CARS  studies  is  found  in  reference  22. 

Water  is  a  major  product  of  air-supported  combustion  of 
hydrogen-containing  fuels,  and  thus  is  available  as  a  ther¬ 
mometric  probe  molecule,  if  present  in  sufficient  quantity. 
Measurement  of  water  concentration  as  a  function  of  distance 
could  yield  chemical  kinetic  information  and  a  better  under¬ 
standing  of  the  combustion  process. 

H2O  is  an  asymmetric  top  molecule  with  three  vibrational 
modes,  of  which  the  iq  symmetric  stretch  is  strongly  Raman 
allowed,  with  a  Raman  shift  of- 3657  cm'1.  Because  it  is  an 
asaymmetric  top,  the  rotational  structure  is  complex  and  the 
rotational  energy  levels  cannot  be  expressed  in  closed  form. 
However,  the  structure  of  the  Q-branch  of  the  tq  is  not  too 
complex  and  the  H:0  CARS  spectrum  has  been  modeled  by 
Hall  et  al.,  (23).  Figure  17  shows  the  comparison  between 
computed  and  experimental  CARS  spectra  of  water  in  a 
heated  cell  at  500  C.  Raman  linewidth  data  were  not  available 
for  this  calculation.  An  aseumed  value  of  0.5  cm-1,  indepen¬ 
dent  of  J,  produced  the  good  agreement  shown.  A  similar 
comparison  for  H2O  in  a  1700  K  premixed  methane/air  flame 
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Figure  18  CARS  water  spectra  at  1700  K 


Folded  BOXCARS.  Aligned  Polarizations 


Figure  19.  Rotational  CARS  spectrum  of  air  using  folded  BOXCARS 


is  shown  in  Figure  18.  Here  the  Raman  linewidth  was  assumed 
to  be  0.2  cm-1.  In  these  calculations,  nearly  one  thousand  vi¬ 
brational-rotational  Q-branch  transitions  were  included  to 
produce  the  computed  spectrum  shown.  Because  the  shape 
of  the  CARS  water  spectrum  is  quite  sensitive  to  temperature, 
the  water  molecule  offers  the  potential  for  combustion  ther¬ 
mometry  over  a  wide  range  of  temperature. 

Pur*  Rotational  CARS  Spoclroacopy 

Pure  rotational  CARS  may  offer  advantages  over  vibra¬ 
tional-rotational  CARS  at  high  pressure  because  adjacent 
rotational  lines  are  better  separated  (by  8  cm-1  in  N2).  than 
in  the  Q-branch  spectrum.  For  the  same  reason — well-resolved 
spectral  lines — rotational  CARS  may  offer  advantages  at  low 
temperatures,  say,  below  800  K. 

Pure  rotational  CARS  presents  the  classic  problem  that 
pure  rotational  Raman  spectroscopy  does;  the  frequency  shifts 
are  small,  and  the  CARS  radiation  is  difficult  to  separate 
spectrally  from  the  input  beams.  One  method  used  in  con¬ 
ventional  Raman  is  to  employ  a  double  (or  triple!  mono¬ 
chromator.  At  UTRC,  the  spectral  separation  of  the  CARS 
beam  is  achieved  by  complete  spatial  separation  of  the  beams 
through  use  of  the  three-dimensional  phase-matching  method. 
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Figure  20  Concentration  variation  of  CO  CARS  spectra  at  1800  K.  model  cal¬ 
culation 
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Figure  2 1  Theory/experiment  comparison  for  CO  CARS  spectra  at  two  dtftarent 
temperatures  and  different  concentrations. 


"folded  BOXCARS”  (12),  described  earlier.  Referring  back 
to  Figure  5,  it  can  be  seen  that  the  CARS  beam,  003  emerges 
well  separated  angularly  from  the  input  laser  beams  and  is 
isolated  easily  from  them  by  simply  masking  and  trapping. 
In  this  manner,  the  CARS  pure  rotational  spectrum  of  room 
air  (300  K,  and  1  atm.)  was  obtained  by  scanning  the  w3  ra¬ 
diation  (Fig.  19).  This  experimental  spectrum  is  in  good 
agreement  with  model  calculations  (72). 

Species  Concentration  Measurements  by  Can 

CARS  measurements  of  species  concentration  in  com¬ 
busting  systems  are  made  by  either  of  two  quite  different 
methods;  from  the  absolute  intensity  of  the  spectrally-inte¬ 
grated  resonant  CARS  signal,  or  in  certain  concentration 
ranges,  from  the  shape  of  the  spectrum.  The  latter  case,  which 


will  be  discussed  here  in  detail,  results  when  the  resonant 
signal  interferes  with  the  background.  At  very  low  concen¬ 
trations,  each  method  fails  because  the  resonant  OARS  signal 
is  completely  dominated  by  the  nonresonant  signal.  However, 
there  are  CARS  techniques  which  have  demonstrably  reduced 
or  suppressed  the  nonresonant  background.  One  of  the  most 
useful  of  these  methods  is  polarization-sensitive  OARS,  which 
yields  nearly  complete  elimination  of  the  background.  But  a 
price  must  be  paid  for  the  background  suppression,  one  which 
can  reduce  the  resonant  signal  considerably.  The  background 
suppression  techniques  usually  employ  integrated  intensity 
measurements  to  make  concentration  determinations.  Be¬ 
cause  the  absolute  strength  of  the  signal  is  difficult  to  measure, 
reference  cells  (which  furnish  a  calibration)  must  lx*  employed. 
This  adds  considerable  experimental  difficulty.  For  these 
reasons,  the  methixl  of  concentration  measurement  based 
upon  spectral  shapes  is  emphasized  in  this  paper.  Both 
methods  have  been  evaluated  critically  in  a  recent  paper  liv 
Kckbreth  and  Hall  (21). 

Concentrations  from  Spectral  Shapes 

CARS  offers  the  unique  capability  of  measuring  species 
concentration,  over  a  limited  range,  solely  from  the  shape  ol 
the  CARS  spectrum.  This  remarkable  situation  applies  in  the 
weak  signal  limit  (Fig.  7)  where  the  resonant  and  nonresonant 
contributions  are  comparable.  This  behavior  is  best  under¬ 
stood  bv  reference  to  Figure  20.  which  displays  a  computer 
generated  set  of  carbon  monoxide  CARS  spectra  at  IHOO  K. 
The  concentration  ranges  from  0.5%  to  100*7;  however,  the 
useful  range  over  which  measurements  can  be  made  is  prob¬ 
ably  only  to  30%  at  the  high  end.  Beyond  this  value  some  sort 
of  reference  is  required  to  scale  the  intensity.  Below  0.5*'}  t  he 
resonant  signal  disappears  into  the  baseline. 

As  an  example  of  spectral  shape  fitting  for  concentration 
measurement,  the  experimental  CARS  spect  rum  of  2. 1  %  CO. 
in  argon  at  300  K,  is  shown  as  the  dotted  line  in  Figure  21.  The 
thin  solid  line  is  the  computer-generated  CARS  spectrum  for 
the  same  concentration.  The  agreement  is  quite  gixxl.  A 
similar  comparison  is  shown  also  in  Figure  21  for  the  case  of 
CO  in  a  CHiAVAr  flame  doped  with  CO.  The  flame  tern 
perature  was  2100  K.  The  experimental  CARS  spectrum  is 
compared  with  a  theoretical  spectrum  for  3.6*7  CO.  a  value 
determined  by  a  quartz  microprobe  and  NDIR  CO  detector 
124).  Again  the  agreement  is  good.  Similar  measurements  by 
the  spectral  shape  method  have  been  made  for  measuring  low 
concentrations  of  water  vapor  as  well.  Before  closing  the  dis¬ 
cussion  of  concentration  measurement  by  spectral  shape,  it 
must  be  pointed  out  that  the  temperature,  and  the  nonreso¬ 
nant  susceptibility  must  be  known  for  a  determination  of 
concentration  of  the  desired  species.  Because  the  nonresonant 
susceptibility  changes  little  (10-15%)  from  reactants  to 
products,  intelligent  estimates  can  be  made  for  xnr  in  a 
combustor,  knowing  the  location  and  temperature  at  that 
location.  From  the  experiments  shown,  it  is  clear  that  con¬ 
centration  measurement  by  spectral  shape  fitting  is  a  viable 
technique,  and,  where  applicable,  should  enjoy  numerous 
applications. 

Background  Suppression 

In  the  most  general  case,  CARS  is  the  mixing  of  three  input 
beams  which  can  have  three  different  frequencies  and  three 
different  polarizations.  Usually,  two  input  frequencies  are 
chosen  to  be  the  same  frequency  w\.  Often  the  polarizations 
are  selected  to  be  parallel  as  well.  However,  because  the  res¬ 
onant  and  nonresonant  susceptibilities  have  different  re¬ 
sponses  to  certain  choices  of  input  polarizations,  it  is  possible 
to  dramatically  reduce  the  nonresonant  CARS  signal  by  use 
of  a  polarization  analyzer  in  the  CARS  beam.  Figure  22  il¬ 
lustrates  this  technique  applied  to  CO  in  a  CO/air-Ar  flame. 
On  the  left  is  shown  the  BOXCARS  signal  using  aligned  po¬ 
larizations  (polarizations  of  all  three  input  beams  parallel). 
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Figure  22.  Polarization  orientation  suppression  of  background  susceptibility. 
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Figure  24  CARS  temperature  measurements  in  a  swirl  burner  with  Jet  A 
fuel. 


BOXCARS  Setup  in  Jet  Burner  Stand 


Figva  23  Experimental  scheme  tor  CAPS  measurements  In  a  combustor. 


A  signal  of  the  “weak  signal”  type  is  obtained  as  discussed  in 
the  previous  section.  In  contrast,  if  the  polarizations  of  the 
input  beams  are  arranged  in  the  pattern  shown,  the  resultant 
signal  is  of  the  “strong  signal”  limit.  The  nonresonant  con¬ 
tribution  is  now  small  compared  to  the  resonant  contribution. 
Earlier  it  was  stated  that  a  price  must  be  paid  for  this  back¬ 
ground  suppression.  Because  of  this  signal  loss,  polarization 
suppression  of  the  nonresonant  background  will  be  less  useful 
than  the  spectral  shape  method  when  the  signal  is  weak,  i.e., 
photon  limited.  For  strong  CARS  signals,  polarization  back¬ 
ground  suppression  is  a  useful  method  and  has  been  applied 
successfully  by  several  workers  (25,26). 

Practical  Applications 

Several  applications  of  CARS  in  practical  combustion  de¬ 
vices,  mainly  CARS  temperature  measurements,  are  listed. 
A  detailed  description  of  CARS  measurements  at  UTRC  on 
a  gas  turbine  test  combustor  is  illustrated  in  some  detail. 

Apparently,  the  first  application  of  CARS  measurements 
on  a  combustion  device  was  performed  by  Taran  and  co¬ 
workers  at  ONERA  in  1978  (27).  Temperature  and  concen¬ 
tration  measurements  of  N2  and  C02  were  made  in  a  simu¬ 
lated  turbo  machine  combustor  burning  kerosene.  More  re¬ 
cently,  this  group  has  made  measurements  in  a  gas  turbine 
combustion  can  (28).  Similar  measurements  at  Wright-Pat- 
terson  AFB  have  been  performed  in  a  bluff-body  stabilized 
diffusion  flame  (29). 

In  experiments  at  UTRC,  Eckbreth  (30)  used  BOXCARS 
to  make  measurements  in  a  50-cm  diameter  combustion  test 
tunnel  fitted  with  a  variety  of  burners,  one  of  which  was  a  gas 
turbine  (JT-12)  combustor.  CARS  temperature  measure¬ 
ments  from  N2  were  made  in  the  primary  zones  of  flames  and 
in  the  exhaust.  Both  gaseous  (propane)  and  liquid  fuels  (Jet 
A)  were  used. 


Tim*  average,  T  =  1450*K  Single  pulse 


— -  Light  wavelength 

Figure  25.  Comparison  o'  single  pulse  CARS  N2  spectrum  with  averaged 
spectrum  in  the  swirl  burner. 


The  experimental  layout  is  shown  in  cross-section  in  Figure 
23.  The  CARS  lasers,  along  with  the  necessary  optics,  are 
contained  in  a  box  which  is  placed  near  the  test  combustor 
tunnel.  Lenses  focus  and  collect  the  BOXCAR  beams  through 
the  tunnel  windows.  (It  should  be  noted  in  passing  that  the 
flame  in  the  test  tunnel  is  extremely  luminous,  so  luminous 
that  most  conventional  optical  methods  would  fail.)  The  laser 
beams  are  manipulated  in  a  receiver  which  spectrally  disperses 
the  beams  and  isolates  the  CARS  signal.  A  reference  PMT 
observes  a  small  fraction  of  the  CARS  signal  and  serves  to 
monitor  optimal  alignment  of  optics.  The  CARS  signal  is 
piped  out  to  a  nearby  control  room  through  a  20-m  long  fiber 
optic.  The  control  room  furnishes  a  much  quieter  environment 
for  delicate  instrumentation,  such  as  the  spectrometer  and  the 
optical  multichannel  analyzer  (which  can  be  very  micro- 
phonic). 

As  examples  of  CARS  measurements.  Figure  24  displays 
averaged  CARS  N2  spectra  at  two  locations  downstream  of 
the  burner  exit  face.  These  CARS  spectra  were  averaged  on 
the  OMA™  for  10-15  sec  which  corresponds  to  100-150  laser 
pulses.  The  spectrum  at  X  *  6  cm  was  made  through  the  fuel 
spray;  hence,  the  relatively  cool  temperature,  the  second  lo¬ 
cation  is  much  hotter.  Figure  25  demonstrates  that  tempera¬ 
ture  can  be  determined  from  a  10~*  sec,  single  pulse  spectrum. 
Comparison  is  made  with  a  130  pulse  averaged  spectrum. 
Although  the  single  pulse  spectrum  displays  photon  statistics 
noise  (shot  noise),  it  is  of  sufficiently  good  quality  to  allow 
instantaneous  temperature  measurement. 

In  the  UTRC  CARS  Laboratory,  in  East  Hartford,  a  mul¬ 
tipurpose  portable  CARS  apparatus  is  being  assembled  for 
use  at  the  Government  Products  Division  (GPD)  of  Pratt  & 
Whitney  Aircraft.  When  complete,  measurements  will  be 
made  in  the  exhaust  of  a  production  run  1 130  (a  modified 
F-100)  gas  turbine  engine.  In  addition  to  N2  temperature 
measurements,  the  concentration  of  water,  carbon  dioxide, 
and  oxygen  will  be  determined,  as  well  as  some  measure  of 
total  unburned  hydrocarbons.  These  measurements  are 
scheduled  to  take  place  before  the  end  of  1981. 

CARS  measurements  in  an  internal  combustion  engine  have 
been  made  by  Stenhouae  et  al.  (32).  The  CARS  input  laser 
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pulses  were  synchronized  with  the  engine  cycle  to  generate  the 
CARS  spectrum  in  time  steps,  hut  at  a  known  point  relative 
to  the  engine  cycle  timing.  Both  nitrogen  and  propane  CARS 
spectra  were  obtained  and  the  temperature  determined. 

More  recently,  workers  at  Ford  Motor  Company  Research 
have  measured  temperature  and  carbon  monoxide  concen¬ 
tration  in  a  research  scale,  single-cylinder  engine  (33).  A 
special  type  of  noncollinear  phase  matching  was  employed 
which  achieved  a  spatial  resolution  of  ~2  mm  along  the  beam 
and  100  fim  transverse  to  the  beam.  Single  pulse,  10-8  second, 
CARS  spectra  were  obtained,  from  which  temperature  was 
determined.  CARS  measurements  also  have  been  made  in  a 
commercial  diesel  engine  (34)  at  Komatsu  Corporation. 
Further  measurements  in  internal  combustion  engines  using 
CARS  are  scheduled  at  many  laboratories  throughout  the 
world. 

Interest  in  CARS  measurements  in  practical  combustion 
devices  is  increasing  rapidly  as  illustrated  by  the  many  labo¬ 
ratories  throughout  the  world  applying  CARS  to  internal 
combustion  engines,  furnaces,  gas  turbine  combustors,  and 
propellant  burning  for  propulsion  and  ballistics.  The  equip¬ 
ment  necessary  for  CARS  is  commercially  available,  and  the 
technique  can  be  engineered  for  application  in  noisy  envi¬ 
ronments,  as  illustrated  by  the  previous  discussion.  The 
capabilities  of  the  CARS  technique  will  be  improved  by  ex¬ 
perimental  modifications  and  from  increased  computer  data 
processing  development  (including  better  modeling  of  spec¬ 
tra).  Based  upon  comparison  with  conventional  diagnostic 
methods,  CARS  should  greatly  increase  our  understanding 
of  fundamental  and  applied  combustion  processes. 
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One  important  problem  in  vibrational  CARS  spectroscopy  concerns  the 
pressure  dependence  of  signatures  which  consist  of  closely  spaced  Q- 
branch  transitions.  As  the  pressure  and  degree  of  line  overlap  in¬ 
crease,  an  interference  effect  termed  colllsional  narrowing  arises  which 
manifests  itself  as  a  blending  or  coalescence  of  the  various  transitions 
into  a  band  whose  width  is  substantially  less  than  one  would  expect  on 
the  basis  of  isolated,  pressure-broadened  lines  [1], 

Colllsional  narrowing  has  its  origins  in  inelastic  energy  transfer  be¬ 
tween  rotational  sublevels,  and  there  is  one  model  of  colllsional  inter¬ 
actions  which  does  not  involve  adjustable  parameters  and  for  which  the 
calculation  of  polarized  Q-branch  spectra  is  particularly  simple;  namely, 
the  Gordon  "extended  J-diffusion"  model  [2].  In  this  picture  of  molecu¬ 
lar  interactions,  collisions  are  presumed  to  be  so  strong  that  a  molecule 
loses  memory  of  its  initial  rotational  state,  with  the  collisions  leading 
to  thermal  randomization  of  the  magnitude  and  orientation  of  the  rota¬ 
tional  angular  momentum. 

Brueck  [3]  has  applied  the  Gordon  model  to  the  problem  of  polarized 
two-photon  resonances  in  fluids  through  an  analysis  that  properly 
accounts  for  the  effect  of  vibration-rotation  interaction.  By  solving 
the  equation  of  motion  for  the  density  operator  with  a  collision  term 
reflecting  the  basic  assumption  about  angular  momentum  thermalization, 
the  vibrationally  resonant  CARS  susceptibility  can  be  shown  for  one 
vibrational  band  to  be  given  by:  13] 


_ ffj _ 

+  tj) 


where  N  is  the  number  density  of  Raman  active  molecules;  a.n  the  polar¬ 
izability  matrix  element;  Afj  the  population  difference  factor  for  the 
transition  (Efj»l);  Ty  a  characteristic  time  for  vibrational  dephasing 
and  relaxation  processes;  wj  end  u>2  the  pump  and  Stokes  frequencies, 
respectively;  end  uj  end  Tj  the  frequency  end  angular  momentum  random- 


ization  time  respectively  for  transition  Q(J).  The  Tv  term  has  been 
neglected  in  these  calculations. 

The  dependence  of  N2  and  CO2  vibrational  CARS  spectra  on  pressure  has 
been  investigated  experimentally  by  measurements  in  a  high  pressure  test 
cell;  details  of  the  experimental  apparatus  are  given  in  [4],  Figure  1 
presents  a  comparison  of  the  measured  N2  signature  bandwidths  at  300°K 
with  those  calculated  from  the  rotational  diffusion  theory  (Eq.  1)  for 
pressures  up  to  100  atm.  Also  shown  is  a  calculation  based  upon  isolated 
line  theory,  that  is,  setting  the  denominator  in  Equation  1  equal  to 
unity.  Good  agreement  has  also  been  obtained  at  higher  temperatures 
('1600°K,  '30  atm.). 
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Figure  1.  Bandwidth  of  300°K  N2  Figure  2.  Bandwidth  of  300°K  CO2 

CARS  Signature  vs.  Pressure.  Reso-  CARS  Signature  vs.  Pressure, 

lution  .4  cm~l.  From  (5]  With  the  Resolution  1.7  cm-*, 

permission  of  North  Holland  Publ.  Co. 

The  CARS  spectrum  of  the  1388  cm'^  v^  mode  of  CO2  has  also  been  inves¬ 
tigated  experimentally  in  the  test  cell  at  300°K.  The  measured  CO2  band- 
widths  show  a  slight  increase  with  pressure  out  to  the  highest  pressure 
recorded,  ~60  atm  (Figure  2),  with  the  predicted  bandwidth  from  the 
rotational  diffusion  theory  in  good  agreement.  The  significant  role 
played  by  collisional  narrowing  can  be  seen  by  Considering  the  width  of 
the  spectrum  that  would  be  expected  on  the  basis  of  isolated  lines. 

*The  experimental  portions  of  the  research  were  sponsored  by  the  Army 
Research  Office  under  Contract  DAAG29-79-C-0008. 
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